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ABSTRACT 



The Similarity Renormalization Group (SRG) is investigated as a powerful yet 
practical method to modify nuclear potentials so as to reduce computational require- 
ments for calculations of observables. The SRG proves to be versatile and robust in 
its treatment of these interactions and opens the door to a deeper understanding of 
the renormalization process. 

Chiral Effective Field Theory (%EFT) provides a consistent and rigorous parametriza- 
tion of the inter-nucleon interaction. While already softer than other available poten- 
tials, transformation via the Similarity Renormalization Group (SRG) brings numer- 
ous computational benefits. The hierarchy of many-body forces inherent in %EFT's 
are also treated consistently by the SRG's simple formalism. The SRG is a natural 
partner to this modern program of formulating the nuclear interaction. 

The key feature of SRG transformations that leads to computational benefits is 
the decoupling of low-energy nuclear physics from high-energy details of the inter- 
nucleon interaction. We examine decoupling quantitatively for two-body observables 
and few-body binding energies. The universal nature of this decoupling is illustrated 
and errors from suppressing high-momentum modes above the decoupling scale are 
shown to be perturbatively small. 

As implemented here, the SRG provides freedom to choose the form of its trans- 
formations and can be tailored to a given application. We explore the impacts of 

ii 



various choices and their decoupling properties, specifically a block-type transfor- 
mation inspired by previous renormalization group techniques. Sharp and smooth 
block-diagonal forms of phase-shift equivalent nucleon-nucleon potentials in momen- 
tum space are generated as examples and compared to analogous low-momentum 
interactions ( u V\ ow ) . 

To explore the SRG evolution of many-body forces, we use as a laboratory a 
one-dimensional system of bosons with short-range repulsion and mid -ransre attrac- 
tion, which emulates realistic nuclear forces. The free-space SRG is implemented for 
few-body systems in a symmetrized harmonic oscillator basis using a recursive con- 
struction analogous to no-core shell model implementations. This approach is fully 
unitary up to induced A-body forces when applied with an A-particle basis (e.g., A- 
body bound-state energies are exactly preserved). The oscillator matrix elements for 
a given A can then be used in larger systems. Errors from omitted induced many-body 
forces show a hierarchy of decreasing contribution to binding energies. An analysis of 
individual contributions to the growth of many-body forces demonstrates such a hi- 
erarchy and provides an understanding of its origins. Several other important sample 
calculations are explored in this model for future use in realistic systems. 

Building on one-dimensional results we performed the first practical evolution 
of three-dimensional many-body forces within the No-Core Shell Model basis. Re- 
sults for the 3 H binding energy are consistent with previous calculations involving 
momentum-space evolution of only two-body forces, and validate expectations from 
calculations in the one-dimensional oscillator basis. When applied to 4 He calcula- 
tions, the two- and three-body oscillator matrix elements yield rapid convergence of 
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the ground-state energy with a small net contribution of the induced four-body force. 
The radius of light nuclei is also explored in the three-dimensional basis. 
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3.9 Calculations of the 4 He ground-state energy using the NCSM. On the 
left is the energy obtained from the NCSM for potentials evolved to 
several different A values as a function of the cut (regulator) momentum 
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3.21 Evolved SRG potentials starting from Argonne Vi$ in the 1 So (left) and 
1 Pi (right) partial waves to A = 1 fm" 1 using a bizarre choice for G s 
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4.1 A snapshot of the embedding process involving the initial (unevolved) 
two-body potential. The far left is a two-body, one-dimensional, po- 
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4.4 Even part of the SRG potential \Vs 2 \p,p') + vj 2 \p, —p')\ in dimen- 
sionless units as a function of p and p' for A = 00, 5, 3, and 2 (where 



4.5 The lowest bound-state energy E 3 for a three-particle system as a func- 
tion of A with the initial two-body-only potential V a . The (red) curves 
with squares include the full evolution of the Hamiltonian while the 
(black) curves with circles use the two-body potential evolved in the 
two-particle system. The right frame shows two additional results from 



4.6 Same as Fig. 14.51 but with initial potential Vr. The right frame shows 



4.7 The lowest bound-state energy E4 for a four-particle system as a func- 
tion of A with the initial two-body potential V a and different initial 



4.8 The lowest bound-state energy E 5 for a 5-particle system as a function 
of A with an initial two-body-only V a potential for several values of 



4.9 Differences of two-body-only and two-plus-three-body A = 3 ground- 
state energies as a function of A. Each of the parameters of the poten- 
tial V a are varied in each plot as the other parameters are held constant. 
The upper panels vary the ranges while the lower vary the strengths; 



shown. Again, 7V max =12 [77 



(A = 1/s 1 / 4 ). The left plot is V a and the right plot is Vp (see Table E]). [80 



A = 1/s 1 / 4 ). The initial potential is V a from Table fl~Tl [83 



varying <j\ and V% from the values in Table 14.11 [84 



two additional results from varying a-i and V2 from the values in Table I4TT1 [86 



three-body force strengths (eg = ±0.05) [87 



N m „ 89 



the left vary the attractive part and the right vary the repulsive part. [90 



4.10 Same as Fig. U but for Vp [91 



4.11 A diagrammatic decomposition of the SRG Eq. (14.61) . A circle at a 

vertex denotes a commutator with T TC \ [92 
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4.12 Contributions from individual terms to the A = 3 ground-state ex- 
pectation value d(V x - 3 ')/d\ for several different two- and three-body 
potentials, as indicated in the plots. We emphasize that A < 2 is 
very small, comparable to A < 1.5 fm" 1 for NN forces in analogous 
calculations with the NCSM [35] 

4.13 Contributions from individual terms to the A = 4 ground-state expec- 
tation value d(V^)/d\ for several different initial three-body poten- 
tials, as indicated in the plots. We emphasize that A < 2 is very small, 
comparable to A < 1.5fm _1 for NN forces in analogous calculations 
with the NCSM [35] 

4.14 Ground-state energy of the A = 4 system with NN interaction evolved 
and then a three-body interaction term (with an exponential regulator) 
fit to the A = 3 ground-state energy. The calculation is repeated with 
a variation of the range parameter, A 

4.15 Same as Fig. I4.14l but at one value of A and a range over the regulator's 
sharpness parameter, n exp 

4.16 Same as Fig. 14. 141 but now the A = 3 ground-state energy and a second 
piece of data, the first excited state for A = 3, are fit to two three-body 
interaction terms. The first term is the same as in Fig. 14.141 and the 
second includes a gradient correction 

4.17 A sequence of plots showing the evolving components of the 3-body 
wavefunction for increasing N max . Each of these plots uses hu = 5. . 

4.18 A sequence of plots showing the evolving components of the 3-body 
wavefunction for fixed JV max and varied hu 

4.19 A sequence of plots showing the evolving components of the 3-body 
wavefunction for fixed hw and varied A^ max 

4.20 A sequence of plots at high A^ max showing the evolving components of 
the many-body wavefunction for A = 2,3, and 4 
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21 Decoupling in the three-particle system using the choice G s = T Te \. The 
initial V a potential is evolved to each A shown in a basis with iV max = 
40. On the left only the two-body potential is evolved and embedded 
while the right involves the full unitarily transformed potential. Matrix 
elements of the potential are set to zero if one or both states have 
N > N cut and the resulting Hamiltonian is diagonalized to obtain the 
ground-state energies plotted 

22 Decoupling in the four-particle system using G s = T rel . The initial V a 
potential is evolved to each A shown in a basis with iV max = 40. On 
the left only the two-body potential is evolved and embedded while the 
right involves the full unitarily transformed potential. Matrix elements 
of the potential are set to zero if one or both states have N > N cut and 
the resulting Hamiltonian is diagonalized to obtain the ground-state 
energies plotted 

23 Decoupling in the three-particle system with G s = if ho- The initial V a 
potential is evolved to each A shown in a basis with iV max = 40. On 
the left only the two-body potential is evolved and embedded while the 
right involves the full unitarily transformed potential. Matrix elements 
of the potential are set to zero if one or both states have N > N cut and 
the resulting Hamiltonian is diagonalized to obtain the ground-state 
energies plotted 

24 Same as for Fig. 14.231 but for the A = 4 binding energy. On the left 
only the two+three-body potential is evolved and embedded while the 
right involves the full unitarily transformed potential 

1 Ground-state energy of 3 H function of the SRG evolution param- 
eter, A. See Table |5] for the nomenclature of the curves 

2 Ground-state energy of 4 He as a function of A. Two hus (left and right) 
and two iV max 's (overlayed) are shown for comparison 

3 Binding energy of the alpha particle vs. the binding energy of the tri- 
ton. The Tjon line from phenomenological NN potentials (dotted) is 
compared with the trajectory of SRG energies when only the NN inter- 
action is kept (circles). When the initial and induced NNN interactions 
are included, the trajectory lies close to experiment for A > 1.7 fm -1 
(see inset) 
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5.4 Ground-state energy of 3 H as a function of the basis size iV max for an 
N 3 LO NN interaction [H] evolved with the SRG using G s = T rci to 
selected A's. The same calculation is shown at two values of fko. ... 

5.5 Ground-state energy of 3 H as a function of the basis size iV max for 
an N 3 LO NN interaction [H] plus an initial NNN interaction [701 [74] 
evolved with the SRG using G s = T rc \ to selected As 

5.6 Ground-state energy of 3 H as a function of the basis size iV max for an 
N 3 LO NN interaction [12] with and without an initial NNN interac- 
tion [TQl [71]. Unevolved ("bare") and Lee-Suzuki (L-S) results with 
Ml = 28 MeV are compared with SRG at ftw = 20 MeV evolved to 
A = 2.0 far 1 

5.7 Ground-state energy of 4 He as a function of the basis size iV max for 
an N 3 LO NN interaction [H] with an initial NNN interaction [701174]. 
Unevolved (bare) results are compared with Lee-Suzuki (L-S) and SRG 
evolved to A = 2.0 far 1 at fku = 28 and 36 MeV 

5.8 Ground-state energy of 3 H as a function of the basis size iV max for an 
evolving N 3 LO NN interaction [4S] without an initial NNN interaction. 
Three values of oscillator parameter, fiw = 28, 20, and 12 MeV are 
shown for comparison 

5.9 Same as in Figure I5TS1 but now including the initial NNN force [701 174] . 
Note the qualitative form of convergence is not affected 

5.10 The 3 H radius as a function of basis size iV max for two choices of Hoo, 
20, and 28 MeV, corresponding to the optimal values for the bare and 
some evolved Hamiltonians. The straight dashed line indicates the 
best converged value at the optimal hu (28 MeV) for the LS effective 
potential at N max =36 

5.11 The 4 He radius as a function of basis size iV max for three choices of 
hu), 20, 28, and 36 MeV (corresponding to the optimal values for the 
bare and some evolved Hamiltonians). The solid lines are calculations 
that include induced NNN forces. The dashed lines are NN-only cal- 
culations. The straight dashed line indicates the best converged value 
at the optimal hw (36 MeV) for the LS effective potential at iV max = 20 
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3 and 4 nucleon amplitudes that occur at v — 2. These terms will 
cancel each other out in an energy-independent formulation. Figure 
courtesy of Ref . [80] 
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A. 2 The leading order 3NFs at v = 3. Figure courtesy of Ref. [5U] .... Il52 



A. 3 Sample 3NF corrections and leading 4NFs at v — 4. A more complete 



A. 4 Leading 4NFs at v — 4. Permutations of vertex ordering and nucleon 



A. 5 Leading NN two pion exchange graphs with explicit A (double solid 



A. 6 Leading order 3NFs with explicit A arising at v = 3, an enhancement 



B.l On the left are shown the harmonic oscillator functions. On the right 
is the truncated momentum delta function 5(p — p')Q{p max — p). In 



B.2 A two-body potential in the momentum representation. The axes are 
initial and final momenta, k and k', from — k max to k max . On the 
upper-left is the original momentum basis two-body potential. The 
other three plots are the same potential after being converted to an 
oscillator basis and back. Notice the truncation of matrix elements 



list of 4NFs is in Fig. IA.41 Figure courtesy of Ref. [80] Il53 



lines are implicit. Figure courtesy of Ref. [87J 1154 



line) arising at v = 2. Figure courtesy of Ref. [80] Il56 



over the deltaless theory. Figure courtesy of Ref. [80J 1 1 56 
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CHAPTER 1 



INTRODUCTION 

1.1 History 

A complete understanding of the interaction between nucleons has been sought 
by physicists for more than seventy years, beginning with the pion exchange theory 
of Hideki Yukawa [1] in the 1930's, which set the paradigm of heavy boson exchange. 
Shortly after the discovery of the pion and confirmation of Yukawa's conjecture, 
nucleon-nucleon (NN) scattering experiments were interpreted as meaning that the 
NN interaction had a strong repulsive core, an intermediate ranged attraction, and 
strong tensor and spin-orbit forces [2]. Efforts to build an accurate NN potential that 
incorporated these features and was based on Yukawa's meson exchange idea were 
eventually quite successful by the late 1990's [3]. However, attempts to apply these 
meson-exchange interactions to many-body problems proved difficult. 

Figure 11.11 shows an example of such an NN potential in both coordinate and 
momentum representations. The hard repulsive core in the left panel and the corre- 
sponding strong high-momentum matrix elements in the right panel represent short- 
distance physics, the details of which are irrelevant to low-energy physics. The specific 
coupling of the low- and high-energy states in the high-energy matrix elements shown 
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Figure 1.1: The Argonne V18 nuclear potential in coordinate (left) and momen- 
tum (right) representations. Notice the hard repulsive core at short-distance/high- 
momentum and small mid-range attractive well. The typical spacing within nuclei is 
indicated by the arrow. 

here inflates the required basis size and hinders attempts to compute larger nuclei 
(A > 4) using ab initio methods. A softened potential, which has the low- and high- 
energy degrees of freedom decoupled, would enable a truncation of the two-body basis 
and ease the computational requirements of ab initio many-body problems. Such soft 
potentials were tried long ago but deemed insufficient, because they could not re- 
produce empirical nuclear matter saturation properties [I]. The problem was that 
many-body forces were not correctly considered. Subsequently, the need for many- 
body forces to accurately describe nuclear systems became apparent. For instance, 
the highly successful NN interactions could not account for the total binding energy 
of the triton [5] . However, it proved to be difficult to systematically account for even 
three-body forces. 
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Then in the early 1990 's Steve Weinberg proposed generalizing chiral perturbation 
theory [6] to develop an expansion consistent with the symmetries of the known un- 
derlying theory of the strong interaction, Quantum Chromodynamics (QCD)P Such 
a program would provide a systematic expansion of contributions, with controlled 
errors at every order, known as Chiral Effective Field Theories (xEFT's). The idea 
sparked an industry in developing and properly parametrizing such a theory of the 
nuclear interaction. Chiral EFTs provide a systematic and controlled expansion of 
not only the NN interactions, but also many-body forces with controlled errors, and 
in a natural framework for including consistent external currents (electromagnetic) 
and relativistic effects. 

Chiral EFTs have provided softer interactions because there is a lower cutoff of 
high-energy physics than in conventional potentials, but ab initio methods are still 
computationally very difficult. Simply using a lower cutoff to further soften xEFT 
potentials has resulted in pathologies in the potentials and significantly diminished 
their accuracy [7]. Recent work that softens these interactions using a unitary block- 
diagonal transformation referred to as V\ ow k have had many successes, but many-body 
forces have only been included approximately. The Similarity Renormalization Group 
(SRG), which generates unitary transformations that can be used to soften potentials, 
can also consistently treat many-body forces and is particularly suited to the needs 
of the xEFT program. 

Now, in an era of rapidly expanding computational capacity, nuclear structure 

physics is experiencing a renaissance and helping to renew collaborations across the 

^ue to asymptotic freedom, QCD is weakly coupled at high energies and thus easier to analyze. 
Unfortunately, it is strongly coupled at low energies and not suitable for perturbative calculations 
in the regime relevant to nuclear structure. 
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larger physics community. Of paramount importance to this program is the devel- 
opment of new techniques to make nuclear many-body problems more tractable in a 
given computational environment. The Similarity Renormalization Group promises 
to be a cornerstone of these techniques for many years to come. This thesis represents 
some of the earliest development of applications of the SRG to problems in nuclear 
physics. 

1.2 Modern Many-Body Nuclear Physics 

The UNEDF (Universal Nuclear Energy Density Functional) collaboration is a 
broad effort within the nuclear structure community to develop a wide-reaching and 
self-consistent model of the structure and behavior of large nuclei based on Density 
Functional Theory (DFT), which describes the bulk properties of many-body systems 
in terms of functionals of the nucleon densities. Because of its scaling properties, 
it is suited for large many-body systems. The DFT strategy has been employed 
successfully for atoms and molecules in quantum chemistry and also in nuclear theory 
using empirical EDF's that are fit to nuclear properties. The UNEDF collaboration 
is attempting to improve the nuclear DFT approach with microscopic input such as 
from Chiral EFT. 

Figure [L2l shows a chart indicating the relationships between the various scientific 
efforts involved in the UNEDF project. The work in this thesis is directly connected 
to the top two bubbles in the chart, and enables work on ab initio DFT's represented 
in the center bubble [8]. The upper left corner represents work on the basic nuclear 
forces in few-body systems and efforts to find optimal and practical parametrizations 
of the nuclear interaction. Some examples listed are Argonne fig (also denoted as 
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Figure 1.2: Scientific interconnections within the UNEDF collaboration. 



AV18) [9J, a phenomenological potential developed in the mid-nineties and widely 
used today, EFT, which refers to the Chiral Effective Field Theories described further 
in section [Dll and appendix \A\ and V\ ow k, which is a renormalization technique related 
to the SRG. The SRG is applied in this thesis to treat the "initial" interactions (such 
as AV18 and xEFT's) to ease computational requirements for calculations elsewhere 
in the UNEDF program. Like V\ ow k, it works to soften potentials and is generally 
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3 (NCSM), which utilizes a finite 



applicable to any interaction. However, unlike V[ ow fc, the SRG handles many-body 
forces in a straightforward way. 

The upper right corner depicts the first step in practical calculations with a given 
initial interaction. In this area of research, ab initio calculations of the lightest nu- 
clei are performed to validate the properties of the initial interactions. This is a 
computationally intensive program requiring large clusters of processors numbering 
into the tens of thousands to accurately compute the nuclei above 6 Li (A = 6). The 
examples listed here are the No-Core Shell Modei 
basis of harmonic oscillator wavef unctions, Green's Function Monte Carlo (GFMC), 
which stochastically computes coordinate-representation integrals numerically, and 
the Coupled-Cluster (CC) method, which is a many-body technique based on a wave- 
function ansatz organized in subclusters of the total system. The GFMC is currently 
limited to nuclei with A < 12 due to its scaling properties and is limited to using local 
potentials like the phenomenological potential, AV18. However, work is proceeding 
to address these challenges [HI [15]. Coupled-Cluster calculations [161 02] have the 
farthest reach up the chart of nuclides of all ab initio techniques and are an important 
source of validation of DFT efforts, but they are limited to a small subset of nuclei 
near closed shells. 

The few-body calculations done in this thesis make use of the NCSM and its 

antisymmetrized basis of harmonic oscillator (HO) wavef unctions. The NCSM basis 

provides for a variational calculation in both the size of the basis, N max , and the 

size of the oscillator parameter, hw. As N max is increased the results converge, and 

at a given N max the oscillator parameter hio must be optimized, balancing infrared 

2 The moniker "No-Core" simply means that all nucleons are active instead of the more traditional 
shell model which considers valence nucleons around an inert closed shell [TUJ [TTJ [T^l H3] • 
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and ultraviolet cutoffs. These features are explored in more detail in appendix [Bj In 
the literature, the name NCSM usually implies use of the Lee-Suzuki (LS) procedure 
to obtain an effective interaction that is not variational in the HO basis parameters. 
Here, we will forgo this method in favor of the SRG, which will provide a softened 
effective interaction and retain the variational features of the basis. 

1.3 Resolution and Renormalization 



Degrees of Freedom Enoro (Motf) 




Col leave Coo rdl n ales 




Figure 1.3: Two pictures illustrating common degrees of freedom (dof), interactions, 
and calculational techniques used in theoretical nuclear physics. 



Figure 11.31 displays an overview of important aspects of theoretical nuclear struc- 
ture physics. The picture on the left is the chart of nuclides on a log-log plot |18j . 
Above the diagonal of the chart are indicated various calculational techniques posi- 
tioned near their relevant regimes, and below are listed interactions that can be used 
in various parts of the chart. The right side of Fig. 11.31 is a picture from the UN- 
EDF collaboration showing some explicit degrees of freedom used in studying nuclear 
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phenomena [19] . Any calculation can be performed using any degree of freedom, but 
many choices may be intractable. We should pick the degrees of freedom most useful 
to the problem at hand. For instance, while Quantum Chromodynamics (QCD) is 
the underlying theory to all nuclear interactions, trying to calculate even the small- 
est nuclei with quark and gluon degrees of freedom is computationally impossible 
with current techniques (though calculating up to A = 4 is a long-term goal for the 
next generation of computers [20] ). Independent checks on calculations from different 
methods is very important to the UNEDF program, and a technique like the SRG is 
crucial to extending the overlap of techniques like DFT, Coupled-Cluster, and NCSM. 

In addition to the choice of gross degrees of freedom, the capacity to adjust the 
parameters of an interaction within a given formulation is desirable. The ab initio 
calculations in this thesis deal with explicit nucleon degrees of freedom with interac- 
tions formulated in terms of relative nucleon momenta. The computational difficulty 
of the few- and many-body problem can be eased with an appropriate treatment of 
this initial interaction. As the number of nucleons increases, the size of the basis 
required for convergence increases very rapidly. Thus, reducing the size of the basis 
required to encode NN physics will reduce the required size of the A-body basis in 
which the NN physics must be embedded. The SRG will provide a convenient way 
to accomplish this adjustment in a manner that preserves initially chosen physics. In 
addition it will also consistently treat many-body forces and other observables. 

Figure 11.41 presents a simple analogy of what the SRG can do for us. The picture 
on the left is a photomosaic image where the image of President Lincoln has been 
made up from hundreds of smaller photos that have been fit to their particular lo- 
cation in the whole image. The picture on the right is just the original photograph 
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Figure 1.4: Two pictures of our 16th president. 



that has been pixelated by standard photo software. Both of these images get the 
basic long-distance information correct - both show Lincoln's face about equally well. 
However, the mosaic image of Lincoln contains a lot of extra information that the 
viewer's brain must work to filter and process, using lots of memory and power. Not 
only is the image on the left harder to use but the extra information it contains is 
even incorrect. All the details in the individual photographs are not actually part 
of Lincoln's face and the viewer has to filter them out. The image on the right has 
much less information but still contains what is needed to recognize the face. Ide- 
ally, a calculation requiring only the long distance information of the right picture 
would benefit from a technique which transforms the mosaic into a version of the 
pixelated right-hand side. This is analogous to the renormalization group technique 
implemented here, in which irrelevant details of the short distance information are 
removed in favor of the dof's consistent with the relevant physical information, sim- 
plifying the interaction for subsequent calculations. The AV18 potential is like the 
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photomosaic on the left because it contains a lot of short distance details irrelevant 
to low-energy nuclear structure physics. The SRG is going to help filter out those 
details and retain the long-distance, low-momentum, information like the picture on 
the right Jf 

1.4 Chiral Effective Field Theories 

As mentioned above, Chiral Effective Field Theories (^EFT's) have evolved to 
play an important role in the development of a quantitative understanding of the 
nuclear interaction. Here we give a brief overview of their features, with various 
additional details involved in constructing xEFT's discussed in appendix lAl 

All these theories are called "effective" because they encode the details of short- 
distance (high-energy) physics in their coupling constants, which are called Low En- 
ergy Constants (LECs). These are the parameters that we must fit to make the 
xEFT's consistent with the data and predictive of low-energy physics. Ultimately, 
short-distance behavior of the strong interaction is governed by Quantum Chromody- 
namics (QCD). The degrees of freedom in a nuclear xEFT are not the same as those in 
QCD because the quarks and gluons are confined at the low energies we are consider- 
ing. However, xEFT's do not need to incorporate explicitly the details of high-energy 
physics in their Lagrangians because the high-energy behavior is parametrized in a 
general way. Thus, xEFT's exploit the fact that any of an infinite number of possible 
parameterizations is sufficient to predict the correct low-energy observables of nuclei. 

3 This analogy breaks down in the sense that a simple low-pass filter works well on the photomosaic 
but not on the potential. The first is a classical problem, and the second is quantum mechanical; 
the photo-mosaic is the observable, while the potential is not. This will be described in more detail 
in chapter [3l 
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Few-body forces are inevitable in any effective theory as a result of elimination 
of degrees of freedom. In any process involving three incoming nucleons there are a 
large number of possible intermediate states, some of which will have high energy and 
occur at very small distance and time scales. Because the details of this high-energy 
physics are irrelevant to low-energy observables, we can use a cutoff to eliminate 
explicit high-energy degrees of freedom. However to account for their effects, we must 
introduce local interactions, including many-body forces. The new LECs that come 
with such terms must be fit to some experimental observables. A major advantage 
of xEFT's is that they provide a consistent and systematic framework with which to 
formulate such low-energy few-body interaction terms in addition to the basic two- 
body (NN) interaction between nucleons. We will consider in appendix |A] the explicit 
form of three- and four-body forces (3NF and 4NF). These terms become important 
for reproducing the properties of larger nuclei (i.e. A > 3). 

Building any good EFT depends on following three basic rules [22] • First, we must 
identify the degrees of freedom we want the theory to encompass and write down the 
most general Lagrangian that obeys the symmetries of the underlying theory (QCD). 
In pionless EFT (EFT (if)), we use only nucleons (m^ ~ 1 GeV) as the degrees of 
freedom; any pion contributions will be parametrized in the LECs. With the pion 
mass (m n ~ 140 MeV) as the cutoff scale this theory will only be valid for nucleon 
momenta well below the pion mass. In %EFT, we will take the cutoff scale to be much 
larger, such as the nucleon or p meson mass. Then we must include pions explicitly 
and the spontaneous breaking of chiral symmetry found in low energy QCD. Hence 
the name "Chiral" in %EFT. The A-isobar (tua = 1232 MeV) might also be included 
in the theory due to its mass being so close to the nucleon. Second, we must declare a 
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regularization and renormalization procedure. Several different types have been used 
with xEFT's including dimensional, cutoff, and spectral regularization procedures 
combined with minimal and power divergence (only used with dimensional regular- 
ization) subtraction schemes. Due to its symmetry preserving properties, dimensional 
regularization is often preferred but has not been successfully applied within the non- 
perturbative resummations involved in calculating nuclear bound states in %EFT. 
Instead, cutoff regularization is used. Finally, we must choose a method of organizing 
the contributions to the potential. The various methods of organization are collec- 
tively referred to as power counting schemes. In any effective theory we integrate 
out higher momentum contributions, and introduce an expansion of correction terms 
to account for the removed degrees of freedom. Thus, we need a formal scheme to 
keep track of the importance of these various contributions and decide which terms 
are important at a desired level of accuracy. The power counting one uses will lead 
to a particular expansion, or hierarchy, of contributions to the total inter-nucleon 
interaction. 

A picture showing a particular power counting scheme for xEFT terms is given 
in Fig. 11.51 The vertical axis in the table is the relative importance of terms entering 
into the amplitude of a process. The order of the graph is expressed in powers of 
momenta (Q/A) where Q is a momentum characteristic of the transfer between inter- 
acting nucleons and A is some large momentum scale associated with the underlying 
physics and hence the scale at which the EFT's predictive power fails|f| Each row 
is seen as a correction to the previous row. Those higher-order terms require higher 

4 Examples of A are: for EFT(7f) A « m T = 140 MeV, for xEFT any one QCD scale such as 
m p ,mN ~ 1 GeV, and for the Fermi theory of weak interactions A = Mw,Mz ~ 80 — 90 GeV 
respectively. 
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Figure 1.5: A table showing the organization of EFT terms in Weinberg power count- 
ing. The figure is from [2Tj . 
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momenta in order to become significant effects and therefore are suppressed at the 
low scale of Q. The columns are separating basic types of terms for multi-body pro- 
cesses. For instance, three body forces (3NF) first contribute at the third-to-leading 
order although in a two-body problem they would never be relevant Note that this 
represents a particular power counting choice originally proposed by Weinberg [6]. 
The graphs displayed in the figure are considered as a perturbative expansion of the 
relevant contributions to the nuclear interaction, but at a desired order in Q/A all 
contributions must be summed non-perturbatively to calculate bound states (in other 
schemes only the leading order terms are summed). 

The various LECs in a xEFT are currently obtained by fitting to experimental 
data. The fit theory can then be used to make predictions about other observables. 
When more LECs enter the expansion, they must be fit to additional experimental 
observables. At the outset, this situation seems to limit the predictive power of 
xEFT's in general, but it is currently the most practical way to proceed with the 
xEFT program of many-body physics within large nuclei. The relevant question 
is whether we can reach a desired accuracy with a number of terms that we can 
handle. For example, hopefully nuclear matter will only need a small set of A-body 
forces (hopefully no more than 4NFs) for an accurate description and therefore a 
relatively small number of LECs must be fit from experiment. Ultimately one would 
like to have an ab initio calculation of nuclear LECs from the underlying theory, 
QCD. Work is currently underway to develop techniques to calculate xEFTLECs in 
this way, and recent progress makes this a plausible future [231 [201 EH E5] . In light of 

5 3N "tadpole" graphs don't contribute to NN amplitudes because N's are not in the EFT. The 
pair production that produces them requires a high energy and has been integrated out. Thus, these 
effects are encoded in the NN contact LECs. 
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these two approaches the xEFT program requires the development of techniques both 
for controlling errors induced by fitting procedures, and for self-consistent A-body 
renormalization capable of dealing with each term in the EFT expansion, preserving 
the input physics. The SRG is a prime example of a technique that satisfies both of 
these needs. 

1.5 Similarity Renormalization Group 

While xEFT's are already significantly softer than many phenomenological po- 
tentials, further softening of these potentials provides increased convergence in cal- 
culations of larger nuclei. A potential like Argonne v\% usually uses a momentum 
mesh out to « 30 fm -1 while xEFT's need around ~ 7fm _1 . Even with the softer 
xEFT potentials the task of calculating properties of 4 He is computationally difficult, 
requiring large spaces in the NCSM basis to achieve convergence (e.g., for A = 4 
at N XDax = 18, the basis dimension is 4750 states and the 4 He binding energy is 
just beginning to converge). However, typical momenta inside a nucleus are of order 
1 fm -1 so we should be able to soften these potentials much further and increase their 
convergence properties. 

A common established technique, called Vj ow fc [2HI EZ] , softens potentials and im- 
proves convergence by transforming potentials to a block-diagonal formjj In contrast 
the SRG, as typically implemented here, evolves the potential into a band-diagonal 



form [28]. A schematic is shown in Fig. 11.61 of the difference in form between a 

6 V\ow k is a momentum space implementation of the block-diagonal transformation idea of Lee and 
Suzuki |29j , though in the literature the term Lee-Suzuki usually refers to a different implementation 
in the NCSM basis. 
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Figure 1.6: A simple illustration of V^ ow fc's block diagonalizing and SRG's band diag- 
onalizing forms. 



potential renormalized by Viowfc (left) and by the SRG (right). Both of these renor- 
malization strategies use a unitary transformation to soften the initial potential and 
thus decouple high- and low-energy degrees of freedom. On the left, Vi ow fc trans- 
forms the potential so that all interaction strength is contained in the two diagonal 
blocks, PVP and QVQ (or simply P and Q blocks), defined by the cutoff scale A[\ 
In contrast the SRG smoothly evolves the matrix element strength towards a chosen 
form, here usually band diagonal in momentum with a width ~ A which indicates 
the amount of coupling left between neighboring states. However, the general form of 
SRG evolution is a matter of choice and we will see that a particular choice reproduces 
Viowk evolution. Note that the V\ ow k renormalization can be achieved through both 
differential and integral procedures, while the SRG currently only exists as a differ- 
ential procedure (as fully explained in chapter [2]). In addition, a V\ ow k treatment of 
7 Usually the V\ ow k Q block is set to zero in practice, but this is not required. 
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many-body forces is not straightforward and limits its practical implementation with 
regard to many-body forces. A quantitative investigation of the decoupling properties 
of the SRG (and the behavior of residual coupling) is presented in chapter [3j 

In the NCSM's oscillator basis, a block-diagonal transformation similar to that 
suggested by Lee and Suzuki [29] has successfully improved convergence properties for 
calculations of many-body observables, it has significant drawbacks which the SRG 
is able to address. This transformation is model-space dependent and not variational 
and extrapolating its behavior in the basis size, iV max , is difficult. The SRG's im- 
proved convergence properties are smoothly variational and therefore lend themselves 
to extrapolation in basis size. This thesis shows SRG calculations to have improved 
convergence over traditional Lee-Suzuki types of treatments. As mentioned above, 
many-body forces are a natural consequence of any renormalization program. Their 
hierarchy is articulated naturally in the xEFT framework as opposed to the more 
ad hoc parametrizations using phenomenological potentials such as IL-IX [30], which 
was designed for use with AV18. The SRG provides for the consistent treatment 
of many-body forces and is seen to work with this natural hierarchy as input and 
complements xEFT's as a general renormalization technique. 

In addition, the formal basis of the SRG allows a term-by-term application to 
analyze its transformation of the initial xEFT hierarchy. This program was begun in 
this thesis by exploring how the SRG preserves the hierarchy of input interaction terms 
in a one-dimensional model and can be easily extended to the three-dimensional case 
within the NCSM framework. A topic of future investigations involves an analytic 
implementation of the SRG to study its preservation of the initial hierarchy and how 
it adjusts LEC's in a realistic xEFT setting. Ultimately, this should provide a deeper 



17 



understanding of renormalization in the broader EFT hierarchy and may lead to 
predictions about the quantitative importance and relation of particular mechanisms 
and processes within nuclear systems. 

While this work focuses on the few-body ab initio program of nuclear structure 
physics, the SRG is a general A-body technique that is also being applied to infinite 
matter and density functional calculations. These methods are currently under heavy 
development and will also depend on the ability to soften the initial Hamiltonian 
input. The SRG is a prime mover in bringing the few-body and many-body regimes 
into quantitative contact with each other. Specifically, the SRG will expand the useful 
range of wavefunction based methods like the NCSM used in this thesis. This will be 
an important check on the validity of the DFT calculations that develop. 

1.6 Thesis Organization 

The following topics in this thesis are original contributions to the field: 

1. A quantitative study of decoupling afforded by the SRG in momentum repre- 
sentation NN forces and the calculational benefits it provides. 

2. The first studies using alternative forms of the SRG to test different types of 
renormalization for nuclear interactions. 

3. Development of a one-dimensional Jacobi coordinate harmonic oscillator shell 
model as a laboratory to explore the behavior of the SRG as it applies to the 
evolution of many-body forces. Within this model we studied: 

• Vacuum expectation value decomposition and analysis of the many-body 
force evolution. 
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• Basis dependence issues and other SRG generators in the oscillator repre- 
sentation. 

• Additional procedures such as fitting three-body forces to evolved two- 
body interactions and evolving individual operators in few-body systems. 

4. First evolution of 3NFs for 3 H and 4 He using an established No-Core Shell 
Model code, based on insight and experience gained from the one-dimensional 
model. 

Decoupling between high- and low-energy states is the feature of SRG evolution 
that leads to increased convergence with basis size in calculations of nuclear observ- 
ables. While first demonstrated for the NN interaction in Ref. [21], in this work we 
achieved a quantitative understanding of its roots and its behavior in the NN partial- 
wave momentum representation. Residual coupling in the interaction above momenta 
of interest is perturbatively small and extends universally to many-body calculations 
like NCSM binding energies. The form of the SRG's transformations can be tailored 
by prudent choice of the operators entering the equations. We show that such choices 
have little effect on the qualitative behavior of decoupling in NN interactions. 

The formal equations of the SRG indicate that decoupling should occur for many- 
body forces in a similar manner as in the NN case. While the SRG can be applied 
in any basis, momentum-basis calculations of three-body potentials are less straight- 
forward because of spectator nucleons (those not participating in an interaction). 
Instead, we choose to work in a discrete basis of harmonic oscillator wavefunctions 
inspired by the no-core shell model. We first built a model in one dimension (ID) 
for the purpose of rapid development and intuition in SRG calculations involving 
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many-body forces. Many different types of calculations have been performed in one 
dimension to gain intuition for their behavior in three-dimensional calculations and 
provide a sanity check when working with the more complex three-dimensional codes. 
It has also provided a simpler environment in which to explore issues of basis depen- 
dence in calculations in general and SRG calculations in particular. 

Finally, we apply the insights gained in ID to the full realistic three-dimensional 
calculation of the NCSM. Here the ID model has proved to be highly predictive of the 
behavior of decoupling in realistic calculations. The SRG improves convergence in 
calculations of 3 H and 4 He, and the input hierarchy of many-body forces is maintained 
as expected from the one-dimensional studies. These evolved many-body interactions, 
which have much improved convergence properties and well controlled errors, will be 
highly sought inputs to various ab initio many-body calculations from across the 
UNEDF collaboration. 

The general layout of this thesis is as follows: Chapter [2] gives a brief introduction 
to the formalism employed in making SRG calculations. In chapter [31 a quantitative 
study is made of the decoupling between high- and low-energy states and its universal 
nature. With an eye towards realistic many-body calculations, chapter H] develops a 
one-dimensional analog to the No-Core Shell Model (NCSM) to explore the behavior 
of SRG evolution in an A-body space. Many different types of calculations are tried 
here to gain intuition for the real three-dimensional case. Chapter [5] applies the insight 
gained from the one- dimensional model to realistic few-body calculations using the 
Jacobi coordinate NCSM. Finally, we conclude in chapter[6]and discuss possible future 
investigations. Several appendices explore technical details of the basis and potentials 
involved in this work. 
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CHAPTER 2 



SIMILARITY RENORMALIZATION GROUP 



The Similarity Renormalization Group (SRG) [32| 133] IM] provides a compelling 
method for evolving internucleon forces to softer forms (281 El] • While observables 
are unchanged by the SRG's unitary transformations, the contributions from high- 
momentum intermediate states to low-energy observables is modified by the run- 
ning transformation. These transformations soften initial interactions and can dra- 
matically reduce the computational requirements of low-energy many-body calcula- 
tions [281 EH ES] • At the same time, the SRG induces many-body forces in response 
to its transformation of the high-energy states in the Hamiltonian. In this chapter I 
present the formalism used for analyzing SRG evolution and the subsequently induced 
many-body forces. 

2.1 Derivation 

We apply the similarity renormalization group (SRG) transformations to inter- 
nucleon interactions based on the flow equation formalism of Wegner [33] • The evo- 
lution or flow of the Hamiltonian with a parameter s is a series of unitary transfor- 
mations 

H s = U S HU\ = T rcl + V s , (2.1) 
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where T re j is the relative kinetic energy, H = T vc \ + V is the initial Hamiltonian in the 
center-of-mass system, and U s is the unitary transformation. Equation (12 .ip defines 
the evolved potential V s , with T rc i here defined to be independent of s in all spaced . 
Then = ^ evolves according to 

as ds ds 

= ^UlU s HU} + U s HUlU s ^ 
ds ds 

= ^UlH s + H s U s ^, (2.2) 
ds ds 

which motivates us to focus on rj s = — ^-U\ =>- = U s ^^-. Therefore, 

ds ds 

^ = Vs H s + H.4 . (2.3) 
ds 

Using UgU] = 1, we can establish the following 

±(U,Ul)=^(iu,)ut + U,(±u!)=0 (2.4) 

=>» Vs + vt = . (2.5) 
Therefore Eq. (12.31) can be written as 

^ = r] s H s - H S T] S = [r] s , H s ] . (2.6) 

The choice of rj s , as long as it is antihermitian, is mathematically free and specifies the 
transformation. A natural choice would involve parts of the Hamiltonian itself which, 
along with a commutator structure, would satisfy the antihermiticity requirement. 
Choosing one commutator argument to be the evolving Hamiltonian, H s , and leaving 
the other choice open, we can write 

Vs = [G s ,H s } . (2.7) 

8 The center of mass kinetic energy, T cm , will not contribute because we are working in a trans- 
lationally invariant system and it will commute with the other operators in Eq. (|2.ip 
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Table 2.1: Various choices for G s considered to date. 



Name 


Description 


T re \ 


Relative kinetic energy between nucleons [28J 


H D 


Running diagonal part of the Hamiltonian advocated 




by Wegner: T re i + V s ai& s S3] 


Hb, fixed 


Fixed, initial diagonal of the Hamiltonian: T rcl + K=o,dia g 


HbD, sharp 


The running block diagonal part of the Hamiltonian. 




P and Q spaces defined by a sharp cutoff at A J3H] 


-f^BD, smooth 


Block diagonal as above but with a smooth regulator of sharpness n 




A strange choice with three different block diagonal regions 




defined by two different cutoff As [36] 


Hho 


The harmonic oscillator Hamiltonian in an A-body space, which is 




diagonal in the HO basis 


T re i + ar 2 


A hybrid H^ type choice for use in the oscillator basis [12] 


T re i + W 2 ) 


The two-body Hamiltonian, suggested as a way to renormalize 




only three-body forces 



The flow Eq. (12.61) now becomes 

^L = [[G S ,H S ],H S ] , (2.8) 
as 

where G s will be referred to as the generator, following conventional usage, though 
technically the generator of the group is rj s . Various possibilities for G s have been 
proposed [3U [281 ESj, though the most common in this thesis is G s = T rc \, for which 
the flow equation simplifies to 

^ = [[T rch V s ),H s ] . (2.9) 

Other choices that have been investigated to varying degrees are listed in Table 12.11 
Equations (12.81) or (12. 9p are operator equations, independent of any basis. They 
can be applied numerically by projecting onto any convenient basis. For two-body 
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systems, a partial-wave momentum basis with states labeled by the initial and fi- 
nal relative momenta, k and k', is convenient, accurate, and straightforward. In 
higher-body systems we will resort to a discrete basisj in which all the states can 
be enumerated up to some basis limiting parameter. The bulk of this thesis will 
employ a basis of harmonic oscillator functions known as the No-Core Shell Model 
(NCSM) [39|, HO, HI] . The NCSM usually uses a renormalization procedure inspired by 
a Lee-Suzuki type transformation [HI [T2J, [13] . Calculations in this thesis will not in- 
clude any such transformations; the SRG is applied to the "bare" or initial potentials 
as given by the various models/xEFT's available. The SRG provides renormalization 
in addition to that intrinsic to the formulation of input potentials (especially those 
from xEFT). 

Note that the dimension of s is 1/ energy 2 . We can define a momentum variable 
A = 1/s 1 / 4 , with units of fm _1 , that measures the spread of the off-diagonal strength 
in units relevant to the basis. Using A as the flow parameter Eq. (12.91) becomes 

^ = ^ [[T rel , V x ] , H x ] , (2.10) 

so that the flow goes from A = oo toward zero instead of starting at s = and going 
toward infinity. Most of our formalism is presented using s for simplicity (since the 
flow is linear in s) but most calculations are done using A. Chapter [3] will demon- 
strate how A is a physically relevant measure of the extent of evolution specifically as 
it applies to the decoupling of high- and low-energy degrees of freedom. Matrix ele- 
ments connecting states with (kinetic) energies differing by more than A 2 are largely 
suppressed and therefore decoupled. 

9 Work is underway to develop a three-body momentum representation code which incorporates 
SRG calculations [37]. 
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The earliest applications to nuclear physics have been in a partial- wave momentum 
basis using G s = T rol [28]. It is a straightforward matter to project Eq. (12. 9p into this 
basis. After some algebra we can show that, for each partial- wave, Eq. (12. 9p can be 
specified in the momentum basis as 

dV s {k',k) 



ds 



9 r°° 

■{k 2 -k' 2 ) 2 V s {k\ k) + - q 2 dq (k 2 + k' 2 -2q 2 )V s (k' , q)V s (q, k) . (2.11) 
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In fact, it is simpler to implement the commutator form of the flow equations on 
platforms like MATLAB or Fortran which can easily handle matrix operations and 
coupled differential equations. Once the Hamiltonian has been projected onto a cho- 
sen mesh (usually Gaussian quadrature), the matrix multiplications follow directly. 

Figures I2.1H2.4I document the evolution of various partial waves from several dif- 
ferent NN interactions in the momentum representation. Note the blue areas, which 
indicate attraction and contribute to binding phenomena seen at low energies, and 
the red areas which indicate repulsion and the hard core probed by large momenta. 
However, the off-diagonal red matrix elements which are coupling low- and high- 
momentum states are also contributing to low-energy observables. The major feature 
to notice looking across the film strips from, left to right, is that the potential is 
driven to a band diagonal form with a width proportional to the square of the mo- 
mentum parameter A. This behavior is a major subject of chapter [3J In all of the 
cases, the SRG is successful at transforming the potential to this band-diagonal form, 
regardless of the specific input physics. Rather the form of the SRG-evolved poten- 
tials is determined by the choice of G s , as shown formally in section 1X2"! The evolved 
potential approaches a universal soft form in a given partial wave, regardless of the 
initial Hamiltonian. This suggests a natural underlying hierarchy of renormalization 
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Figure 2.1: Film strips of evolution for representative partial waves in the N 3 LO (500 
MeV) interaction of Ref. in the momentum representation. Partial waves shown 
are %, 3 S U l P x , and 3 F 3 . 
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Figure 2.2: Film strips of evolution for representative partial waves in the N 3 LO (600 
MeV) interaction of Ref. in the momentum representation. Partial waves shown 
are %, 3 S U l P x , and 3 F 3 . 
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Figure 2.3: Film strips of evolution for representative partial waves in the N 3 LO 
(550/600 MeV) interaction of Ref. [20] in the momentum representation. Partial 
waves shown are 1 ) So, 3 Si, l Pi, and 3 F 3 . 
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Figure 2.4: Film strips of evolution for representative partial waves in the AV18 
interaction of Ref. [H] in the momentum representation. Partial waves shown are ^So, 
3 Sx, ^i, and 
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particularly suitable to the xEFT program, though it is equally successful at softening 



relatively hard potentials such as AV18 (see Fig. 12.41 
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Once we have obtained wavefunctions of the evolved Hamiltonian through the 
flow equations of Eq. (12. 8p . we can directly construct the unitary transformation, U s , 
at each s explicitly [31]. Given a complete set of eigenstates of H s , {\ip a (s))}, U s can 
be written as, 

U s = J2\Ms))(M0)\ ■ (2-12) 

Q 

Because of the discretization, the summation over a is finite and can be achieved 
through simple matrix multiplications using platforms such as Fortran and MATLAB. 
With U s , the wavefunctions can be evolved independently, 

u.iMV) = ^l^(s))(^(o)|^(o)) 

a 
a 

= \M*))- (2-13) 

This form is particularly convenient for evolving individual operators in parallel with 
the chosen initial interaction. In particular, any operator can be evolved directly as 

O s = U s OU\ , (2.14) 

as opposed to solving another flow equation, 

dO* 



[[Gs,H B ],0 8 ] , (2.15) 



ds 

just as for the Hamiltonian. These are identical strategies, but for some operators 

evolving the potential and using U s as in Eq. (12.141) is more numerically robust |43j . 

10 Note however that SRG evolution does in general introduce nonlocality to the initially local 
AV18, spoiling its use in the GFMC program. 
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Ultimately, having the explicit unitary transformation matrix facilitates rapid calcu- 
lation of many different operators with minimal coding. 

To carry out the SRG evolution, we use a built-in M ATLAB differential equation 
solver, such as the MATLAB function ode23 |44], which is an implementation of a 
Runge-Kutta differential equation algorithm [45]. We studied the running time to 
evolve the various SRG schemes (i.e., the choice of G s ) by plotting the time to run 
versus the evolution parameter, s. We find a linear relationship to s, indicating no 
stiffness, in every combination of potential and SRG scheme used to date. Details 
regarding these timing studies and other computational scaling issues are reviewed in 
Appendix [Ql 

2.2 Second Quantization 

To see how the two-, three-, and higher-body potentials are identified, it is useful 
to decompose H s in second-quantized form [36] which was first used for this purpose 
in Ref. [17]. We can write an A-body Hamiltonian in general as: 

H = Ti i aj i a i + ^2 V ijkA a ) a i a k + ^ J] Vifkimn4 a Wk a na m ai + • • • , (2.16) 

ij ijkl ijklmn 

where a\ and represent creation and destruction operators with respect to the 
vacuum in some single-particle basis. The quantities Tij, V^, and V^ lmn represent 
matrix elements of their respective operators. 

As a straightforward way to see how the SRG induces many-body forces in an A 
body space, one can evaluate the flow equations using the Hamiltonian in the form 
of Eq. fl2~T6D : 
dH. 



ds 



- ij ijkl 



Tija\ aj + v m,s a H a i a k + ■■■ 

ij ijkl 



(2.17) 
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For example consider the evolution of the potential, V s , with the Hamiltonian 
truncated to just the two-body interactions: 



ds 



T H a i a 3> Yl V k\Ln.A a War 
ij ijkl 



ijkl 



\f t t t ft t \ t t 



^ij^klmm^opqr ^ f Q*&H&u(%] &n(%rn.(^(^&r 



^ KJbj KAi^KAiJ^ KAi'yi^* l, YYL q p T Q 

— a],apa r aqalaja\aia n a m + aj,aj,a r a 9 aj„a|a n a m afa.,J , (2.18) 

resulting in the large strings of a's and a^s. 

These terms can be brought into normal order with respect to the vacuum to 
simplify them. For this we will employ Wick's theorem [IS], which states (using Ai 



to represent or a 



AiAjAk ■ ■ ■ A m 



N(AiAjA k Ai ■■■A m ) 

+N ^{AiAjA^Ai ■ ■ ■ A m ) + all single contractions^ 
/ i — i i — i 

+N I (AiAjA k Ai ■ ■ ■ A m ) + all double contractions 



+N (all full contractions) 



(2.19) 



and its corollary 



i 1 



NiAAjWAkAti = N(A i A j A k Ai) + N^A^Ai) + N{A i A j A k Ai) 



+N(A i A j A k A l ) + N(AiAjA k Ai) + N(AiAjA k Ai) 
+N(A i A j A k A l ) , (2.20) 
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which helps to simplify products of multiple, individually normal-ordered operator 
strings. There is one non-vanishing type of contraction in the vacuum 

<3] = , (2.21) 

and the other three, a\aj, a^L/, and ajaj, give zero. 

Returning to the result of Eq. (12.18p . we can now see what many-body contribu- 
tions it contains. The normal order of all operators in each term is the same for all 
terms and thus cancel, so no-five body operator survives at this order. All terms re- 
sulting from a single contraction also cancel each other because of the double commu- 
tator structure. For example N(a\oja\aja n a m ala p a r aq) and N(—a\a p a r a q a\a^\a\a n a m ) 
are equal and opposite. At the level of double contractions the commutator cannot 
cancel all the terms and three-body terms survive. For example the first term gives: 

ttt ft ( t t t ttf 

ttt XX ttf XX 

-\-Q,^Qi^Q, Qi n Cl r Qiqdjj t d rn p -\- CL^CL^ Q, CL rn Q, r Clqdji t d n p 

ttf XX ttt XX 

+a\a\a\a n a r aq8 j i8 mp + a\a\a\a m a r aq8ji5 np ) , (2.22) 

which can be simplified due to the symmetry of the potential matrix elements. In 
fact all of these terms are permutations of a single contribution in a diagrammatic 
representation shown in Fig. 12. 51 Equation (12.221) lists the explicit legs of the second 
diagram on the bottom, which is the lone two-body only contribution to the induced 
three-body force. Carrying out the same calculation for all triple contractions will 
result in enumerating the legs associated with the two diagrams in the top right of 
Fig. 12.51 which contribute to two-body evolution. In an A-body system this means all 
disconnected combinations of these two-body contributions. 

33 




Figure 2.5: A diagrammatic decomposition of the SRG induced forces (17j. A circle 
at a vertex denotes a commutator with T rc \. 

As follows from these normal-ordering calculations, A-body forces cannot be af- 
fected by the value of any initial higher-body terms; e.g., two-body evolution cannot 
be influenced by three-body forces. This can be seen in second quantized form by 
writing out commutators involving an A-body force. One can quickly see that it 
is impossible to contract such terms to have less than A pairs of a''s and a's. So 
for instance, an initial four-body force cannot contribute to the evolution of two- or 
three-body forces. This feature is true for calculations done in the vacuum. However, 
when the SRG is implemented in the background of a nuclear medium, the rules for 
contractions change. This has profound consequences for the use of the SRG in nu- 
clear matter calculations. Those calculations are outside the scope of this thesis, so 
the "in-medium SRG" will not be considered further here. 

The correspondence between equations like Eq. (12.221) and Fig. 12.51 can be written 
more concisely by expanding the commutators of the SRG flow equation, Eq. (12. 8|) . 
For instance the top row of Fig. 12.51 corresponds to the expansion of Eq. (12. 8p in 
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2-body space: 



[[T Teh V s ],H s 



ds 

= [(T Tel V s -V s T iel ),(T Iel + V s )} 



V { *\(T Tcl + V s ] 



V {2) T 

V s > 1 rel 



= Vf + vfWW-VWV®, (2.23) 

—(2) =( 2 ) 

where V s represents a single commutator with T TC \ and V s represents a double 
commutator with T re i, [[T re i, Vs ], T re J. These correspond to the vertices with one 
and two circles in the top row of Fig. 12. 5[ so the last line of Eq. (12.231) is represented 
by those diagrams. Similarly the bottom row represents the expansion in the 3-body 
space, 

dV s {3) =(3) r _( 2 ) 



ds 



v [ s + [Vf\ y?>] + [Vf\ k (3) ] + [vf\ K (2) ] + [vf\ v™] . (2.24) 



We stress that these diagrams are not Feynman diagrams, but merely a mnemonic 
to summarize the contributions to the SRG transformation. However, each vertex 
does represent a particular type of interaction and the diagrammatic treatment can 
be used to explore the interplay within the flow equations. An example of this kind 
of investigation will be presented in chapter HI 

Disconnected diagrams for each A-sector are canceled out by the evolutions de- 
fined in lower sectors. Consider the flow equation for the general three-body Hamil- 
tonian [28] (in a modified notation), 

dV s dV 12 dV 13 dV 23 dV 123 u , , , , 
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where V^ is the two-body potential between the ith and jth particles and V123 is 
the three-body potential. The relative kinetic energy T rc \ is the combination of the 
relative kinetic energy between the ith and jth particles, Ty = {ki — kj) 2 /2m, and 
the relative kinetic between that pair's center of mass and the third particle, = 
+ kj)/2 — kk) 2 /m. Because Tj and Vjk commute the commutator of T rc \ and 



is 



Prei, Vij] — [Tij, Vy] , (2.26) 
so the flow equation for the two-body system is 

^L = [[Tij,Vij],(Tij + Vij)]. (2.27) 

We can expand the three-body flow equation, shown in Eq. (12.251) . and using Eqs. 12.261 
and 12.271 see that the two-body evolution cancels out precisely, leaving 

^ = [[T 12 ,V 12 ],(T 3 + V 13 + V 23 + V 123 )] 

+ [[T 13 , Vi 3 ], (T 2 + V l2 + V 23 + V l23 )] 
+ [[T 23 , V 23 ], (Ti + V12 + V 13 + V123)] 
+ [[T rch V 123 ],H s ] . (2.28) 

So the three-body evolution contains no spectator graphs, but only connected commu- 



tators of the form [Vij, Vjk], [V^, Vyk], or [V^, Vijk]- The two-body evolution equations 
have completely cancelled out of the three-body sector and are completely determined 
by the two-body evolution. In other words, free-space A-body forces are defined by 
evolution in the A-body space and are not altered or evolved differently in higher 
sectors. 

We can also see these cancellations in second-quantized form by considering the 



one-loop two-body graphs on the top of Fig. [275] coup led with a single line representing 
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a spectator particle in a three-particle system. Such graphs cannot be generated from 
Eq. (12.181) . there not being enough contraction opportunities to generate the necessary 
creation/annihilation operators and momentum delta functions. 

So, the SRG induces A-body forces using connected terms involving lower-body 
forces; three-body forces can be initially zero but will grow due to two-body-only con- 
tributions. Eventually A-body forces are induced as one can see by taking successive 
infinitesimal steps in the evolution. In the first step, three-body forces appear which 
allow induction of four-body forces in the next step. These in turn allow higher- 
body forces ad infinitum. Including induced A-body forces is a unitary calculation 
in an A-body space. In a two body-space, three-body forces cannot be represented 
so evolution in the two-body space is only approximately unitary when the evolved 
Hamiltonian is used in the three-particle system. Because the induced three-body 
forces are not included, we speak of "missing induced A-body forces" . These induced 
forces appear in a hierarchy of decreasing strength, which will be explored quantita- 
tively in chapter HI 
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CHAPTER 3 



DECOUPLING IN THE MOMENTUM 
REPRESENTATION 

3.1 Benefits of Decoupling 

While observables are unchanged by the SRG's unitary transformations, the con- 
tributions from high-momentum intermediate eigenstate components to low-energy 
observables is modified by the running transformation. In particular, the SRG as 
implemented in Refs. [28l I3T] has the effect of partially diagonalizing the momentum- 
space potential to a width of order the evolution parameter. Because of this partial 
diagonalization, one anticipates a direct decoupling of low-energy observables from 
high-energy degrees of freedom. 

In Ref . [21] , evidence for decoupling at low momentum was shown for the Argonne 
^18 [H] potential in calculations of phase shifts and the deuteron. In this chapter, we 
extend the demonstration of decoupling to nucleon-nucleon (NN) potentials from 
chiral effective field theory (EFT) [IHJ ED] arid to few-body nuclei up to A = 6 to 
verify its universal nature and to show quantitatively that the residual coupling is 
perturbative above the energy corresponding to the SRG evolution parameter. 
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The practical test for decoupling is whether changing high-momentum matrix 
elements of the potential changes low-energy observables. The strategy here is to 
first evolve the initial potential Vnn with the SRG equations to obtain the SRG 
potential V s , where s denotes the flow parameter of the transformation. Then we 
apply a parametrized regulator to cut off the high-momentum part of the evolved 
potential in a controlled way. This cutoff potential is used to calculate few-body 
observables and their relative errors. By varying the parameters of the regulator and 
correlating them with errors in the calculated observables, we have a diagnostic tool 
to quantitatively analyze the decoupling. 
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Figure 3.1: Phase shifts and relative errors in the 1 So channel for SRG potentials 
evolved from the Argonne t^g potential of Ref. [H]. The upper row shows cut and 
uncut potentials and phase shifts for the unevolved initial potential. The lower row 
shows the same for the initial potential evolved to A = 2.0 fm~ . 
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Shown in Figure 13.11 is an example of the strategy employed in this chapter to 
study the decoupling benefits of the SRG. On the top row, the pictures from left to 
right are the uncut Argonne v±$ potential in the partial wave, the same potential 



row shows the same three plots except now the Argonne potential is evolved with the 
SRG to A = 2.0 fm -1 . In the top phase shift picture, the calculation from the cut 
potential fails everywhere; the low-energy observable is lost despite the low-energy 
matrix elements of the potential being preserved. In the bottom phase shift picture the 
low-energy phase shift is preserved because the SRG has renormalized the potential 
and preserved all the information relevant to the calculation. 

The phase shift, Se{Ek ), (where £ denotes a particular partial wave) is a nuclear 
scattering observable calculated non-perturbatively from the potential, via the rela- 
tion Re(k, k'\ Efo) = — tan ^ - [5T] , where Ri is the reaction matrix, 



Where Vi represents the potential in a particular partial wave, and k is the mo- 
mentum corresponding to the desired on-shell energy of the calculation. The P here 
indicates that we are using the Cauchy principal value prescription which keeps matrix 
elements real, a convenient property for numerical calculations. Cutting the potential 
as shown in Fig. 13.11 does not lose information for the first order term (linear in V) 
since it represents only one matrix element on the diagonal of the potential. However, 



cut at A = 2.1 fm 2 and the phase-shift calculated from both potentials. The bottom 




(3.1) 
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the second (and subsequent) terms involve sums over intermediate states so that ma- 
trix elements along the whole row or column involving the state k$ are necessary to 
the calculation. When these are arbitrarily cut off as in the top row, the phase shift 
will fail at any energy. 

However, in the bottom row of Fig. l3.1l the SRG has transformed, or renormalized, 
the potential so that all the relevant information has been placed inside a low-energy 
region of the potential. Now we can cut the high-energy matrix elements, which 
are decoupled from the low-energy ones, leaving the phase-shift intact for low ener- 
gies. Thus we need less of the total basis to represent the low-energy physics we are 
concerned with reproducing. 

The tool we will use to study decoupling is a smooth exponential regulator applied 
to the potential to cut off momenta above A: 

Vx,\(k, k') = e-^ A ^ nexP V x (k, £;')e- (fc ' 2/A2)ncxp , (3.2) 

where n exp takes on integer values. From the cut potentials we calculate observables 
such as phase shifts and ground-state energies and compare to values calculated with 
the corresponding uncut potential. If there is decoupling between matrix elements in 
a given potential (evolved or otherwise), we should be able to set those elements to 
zero in this systematic way and use the relative error in the observable as a metric of 
the degree of decoupling. By varying n exp we can identify quantitatively the residual 
coupling strength. 

Here we are working with NN interactions only and therefore the SRG transfor- 
mations are truncated at the two-body level, which means that they are only approx- 
imately unitary for A > 3 such that those observables will vary with s. In these cases 
decoupling is tested by comparing cut to uncut potentials at a fixed s. All two-body 

41 



observables calculated with the uncut H s are independent of s to within numerical 
precision. The actual numerical error depends on the details of the discretization 
(e.g., the number and distribution of mesh points, usually gaussian) and on the ac- 
curacy and tolerances of the differential equation solver. While in practice we can 
make such errors very small, to avoid mixing up small errors we will also compare 
cut to uncut potentials rather than to the unevolved (s = 0) potential for two-body 
observables. 

3.2 Mechanics of Decoupling 

The source of decoupling is the partial diagonalization of the Hamiltonian by 
the SRG evolution |52j. For the NN interaction, the flow Eq. (12.81) can be simply 
evaluated in the space of discretized relative momentum NN states [28]. For a given 
partial wave, with units where h 2 /M = 1, we define diagonal matrix elements of 
momentum k aj^ 

(k\H a \k) = (k\H D \k) = e k , (3.3) 

and 

(A;|r rel |fc} = e k = k 2 . (3.4) 
If we take G s = T re i, the flow equation for each matrix element is 

^- s (k\H s \k'} = ^2(e h + e k ,-2e q )(k\H a \q)(q\H a \lJ) 

= ~{e k ~ e k >f{k\V s \k') + ^(e fc + e w - 2e q ){k\V s \q){q\V s \k') .(3.5) 

i 

A useful image of the mechanism by which the SRG diagonalizes the potential is 



shown in Fig. 13.21 The pictures represent a snap shot of the different parts of the flow 

11 Note that the derivation here is only valid for a discretized momentum basis such as the gaussian 
mesh used here. 
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Figure 3.2: A snapshot of the contributions to the flow of the potential when using 
the G s = T re i choice of SRG. From left to right, the pictures are V\=2.5> 1st and 2nd 
rhs terms from Eq. (13. 5p . and Vx=i.5. The color scale for the middle pictures is scaled 
up for visibility. 



equations contributing to the evolution between two specific values of A. The plot on 
the left is the potential evolved down to A = 2.5 fm" 1 and the plot on the far right is 
the potential further evolved to A = 1.5 fm -1 . The two panels in the middle are the 
first and second terms of the right-hand side from the flow equations in Eq. (13. 5p . It is 
clear that the first term is the dominant contributor to the suppression of off-diagonal 
matrix elements outside the band of width A 2 . At the same time the second term 
is providing for the flow of physics information to low-energy states by increasing 
the magnitude of matrix elements associated with those states. While the first term 
achieves decoupling by suppressing off-diagonal matrix elements, the second term is 
preserving unitarity by transforming the information of high-energy matrix elements 
to low-momentum states. 

We can find a semi-quantitative approximation for the flow of off-diagonal matrix 
elements by keeping only the first term on the right side of Eq. (13.51) . Then the flow 
equation applied to individual off-diagonal matrix elements simplifies to 

±(k\H s \k') = ^(k\V s \k>) « -( Cfc - e k ,) 2 (k\V s \k') , (3.6) 
as as 
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which has the simple exponential solution 

(k\V s \k') « (k\V s=0 \k')e- s ^-^ 2 



(3.7) 



In Fig. I3.3[ we plot (k\V s \k'} and the approximation from Eq. (13.71) versus s for 
some representative off-diagonal points in two partial waves. In almost all cases the 
approximation gives a reasonable estimate of the monotonic decrease to zero; in the 
one exception there is a significantly more rapid decrease than predicted. 
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Figure 3.3: Absolute value of the matrix element (k\V s \k') for a representative sam- 
pling of off-diagonal (k, k') pairs as a function of s, compared with the simple solutions 
from Eq. (13.71) . which are straight lines (they agree at s — 0). The partial waves 3 Si 
and 3 Pi are shown. 



Equation (13. 7p shows that it is convenient to switch to the flow variable A = 1/s 1 / 4 , 
which has units of fm -1 , because it is a measure of the resulting diagonal width of 
V s in momentum space. More precisely, the matrix (k\V s \k') plotted as a function of 
kinetic energies k 2 and k' 2 will rapidly go to zero outside of a diagonal band roughly of 
width A 2 , which is verified by numerical calculations [521 [53]. For momenta within A 

44 



of the diagonal, the omitted quadratic part of the flow equation is, of course, essential, 
and drives the flow of physics information necessary to preserve unitarity. 

We can examine the flow of off-diagonal matrix elements as a whole by considering 
the trace of H 2 , Tr[# 2 ] = Tr[U s HUlU s HUl] =Tt[H 2 } using UW = 1 and the cyclic 
property of traces [33J: 

hi 

= ^EKW>I 2 + ^EKW>I 2 > (3.8) 

and therefore we get, 

^Ek«>i 2 = -^Ekw>i 2 = -2j2m s \i)-[ s (i\H s \i) . (3.9) 

In the case of Wegner's choice, G s = Hp, the flow equation for each matrix element 
is 

j- s (k\H s \k'} = ^(e fc + e k , - 2e q )(k\H s \q)(q\H s \k'} . (3.10) 
then (I3.9P simplifies to 

^Ekw>i 2 = ^E^-^iw^i 2 

= -2^e fc (e fe - e q )\(k\H s \q)\ 2 - 2^e,(e 9 - e fc )| (g|# s |£:>| 2 

= -2^(e k - e q ) 2 \(k\H s \q)\ 2 , (3.11) 

which shows a sum of positive definite terms whose derivative is negative, so that the 
off-diagonal matrix elements decrease with s. On the other hand using the choice 
G s = T rol and Eq. (13. 5p to simplify Eq. (13.91) we get 

^EK^'>i 2 = -4 y £ l <*fa-e 9 )\(k\H.\q)\* 

= -2j2(e k ~ e q )(e k - e g )\(k\H s \q)\ 2 , (3.12) 
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so that off-diagonal elements are not guaranteed to decrease monotonically if — e q 
and e k — e q have opposite signs (see Ref. [51] for details). However, this does not 
happen in the range of s that has been considered in the nuclear case because of 
the dominance of the kinetic energy. T is much larger than V on the diagonal and 
therefore T is effectively Hp in the nuclear case. We should see the difference when 
we evolve down to a scale at which T is no longer dominant. At such a scale new 
bound states appear and the choice G s = T re i cannot order them properly [54| . We 
don't see any such pathologies in the nuclear case because we have never needed to 
evolve so far in A as the scale at which the deuteron appears, ~ 0.5 fm -1 ^. 

3.3 Phase Shift Errors 

In the upper-left panel of Fig. 13.41 we show results for the 1 So phase shifts vs. 
energy calculated using the unevolved 500 MeV N 3 LO potential of Ref. [49] and the 
corresponding SRG potential evolved to A = 2.0 fm -1 and then cut using the regulator 
of Eq. (13.21) with n = 8. We do not explicitly show results from uncut SRG potentials, 
because they are indistinguishable from the unevolved results. 

The qualitative pattern is that when the regulator parameter A is greater than 

A, there is good agreement of phase shifts from uncut and cut potentials at small 

energies and reasonable agreement up to the energy corresponding to the momentum 

of the cut, E'lab ~ 2A 2 /m (with h = 1). When V^ rg is cut below A, there is poor 

agreement everywhere and the phase shift is zero above this energy (e.g., above Ei^ = 

100 MeV for A = 1.1 fm _1 ). Thus the decoupling of high and low momentum means 

that we can explicitly cut out the high-momentum part of the evolved potential 

12 In fact, we will not evolve this low also because we expect the hierarchy of many-body forces to 
break down at such small A's, as is discussed in detail in chapter's [4] and [5] 
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without significantly distorting low-energy phase shifts as long as we don't cut below 
A. Cutting out the high-momentum part of conventional nuclear potentials does cause 
distortions, which has led to the misconception that reproducing high-energy phase 
shifts is important for low-energy nuclear structure observables [3T] . 




Momentum of Cut, A [fm ] (n=8) Momentum of Cut, A [fm ] 



Figure 3.4: Phase shifts and relative errors in the 1 So channel for SRG potentials 
evolved from the N 3 LO (500 Me V) potential of Ref. [19]. The upper- left graph shows 
the phase shifts vs. energy for the uncut A = 2 fm -1 potential and several cut versions 
with n — 8. The other panels show the relative error as a function of the momentum 
cut parameter A at various energies E, A's, and n's, respectively. 
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The quantitative systematics of SRG decoupling are documented in the other 
panels of Fig. 13.41 where we look at the relative error as a function of the cutting 
momentum using log- log plots. In these error plots, three main regions are evident. 
In the region below the A corresponding to the fixed energy, the predicted phase shift 
goes to zero since the potential has vanishing matrix elements, so that the relative 
error goes to one. Starting at A slightly above the value of A, there is a clear power-law 
decrease in the error. In between is a transition region without a definite pattern. 

We focus here on the power-law region. In the lower-left pane, we find that this 
decoupling starts with a shoulder at momenta slightly above A. This effect saturates 
when A becomes comparable to the underlying cutoff of the original potential (see 
Fig. 13.61) . In the upper- right pane we see that the shoulder signaling the start of the 
power-law decrease is not affected by the energy, E. This holds for other values of 
A and n. In the lower-right pane we vary the exponent of the regulator, n, which 
changes the smoothness of the regulator. The smoothness affects the slope of the 
power law and the fine details in the intermediate region, but does not change the 
position of the shoulder near A. As discussed below, the power-law behavior in the 
relative error signifies perturbative decoupling with a strength given by the sharpness 
of the regulator used to cut off the potential. 

Indeed, the behavior of the errors in the decoupling region, where A > A, can 
be directly understood as a consequence of the partial diagonalization of the evolved 
potential. The calculation of the phase shift at a low-energy k 2 <C A 2 will involve 
an integral over p of V\^(k,p). But the potential cuts off the integral at roughly 
p 2 pa k 2 + A 2 < A 2 , which means that we can expand the difference in the uncut and 
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cut potentials: 

Xto + jpi)v(k,p). (3.13) 

Simple perturbation theory in 5V then predicts the dependence of the phase shift 
error to be 1/A 2n , which is the power-law dependence seen in Figs. 13.41 and 13.51 The 
accuracy of first-order perturbation theory is evidenced by the constant slope of the 
error curves, which translates into perturbatively small residual coupling. 

The detailed dependence on the energy and A is not so trivially extracted. How- 
ever, the weak dependence on < 100 MeV and strong dependence on A < 3 fm _1 
at fixed A seen in Fig. 13.41 implies that the integration picks up the scale A, so that 
the dominant error scales as (A/A) 2n . This is, in fact, observed numerically for inter- 
mediate values of A (e.g., for 1.8fm~ < A < 2.8 fm -1 when A = 3fm _1 ). 

We checked this decoupling behavior in different partial waves and for other N 3 LO 
potentials and found the same perturbative region in all cases. The plots of Fig. 13.41 
are reproduced in Fig. 13.51 for representative partial waves. The potential in the S 
waves typically passes through zero for momenta in the region of A = 2fm _1 (as in 
Fig. 12. ip . which might lead one to associate decoupling with this structure. The error 
plots for other partial waves that lack this structure show that it is a more general 
consequence of the SRG evolution. 

Note that in higher partial waves, such as 3 F 3 , the phase shift is already well 
decoupled and therefore doesn't benefit as much from the SRG evolution. Errors 
shown for this partial wave in the bottom strip of Fig. 13.51 are below 10~ 4 at the start 
of the decoupling region (note the shoulders at A in the center panel). Such partial 
waves are the large angular momentum components of the partial wave decomposition. 
Due to the centrifugal barrier term they probe long-distance forces and have strength 
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Figure 3.5: The phase shift errors computed in select partial waves. Other channels 
exhibit the same power-law dependence of the error for A > A. 
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only at low momenta. Therefore, low-momentum and high-momentum states are 
automatically decoupled in these channels and little change is provided by the SRG 
here. 




Momentum of Cut, A [fm ] 



Figure 3.6: The phase shift errors computed in the channel of the Argonne Vis 
potential. Here the onset of the power-law slope is much less saturated at higher A's, 
because this potential has higher-momentum components to be renormalized. 

As noted above, the shoulder in the relative error log-log plots corresponding to 
the onset of the decoupling region saturates at large A due to the high-momentum 
cutoff introduced into the initial potential during its initial formulation. In the plots 
of Fig. 13.41 the shoulders for the A = 4 and 5 curves occurred at about A = 3 since the 
xEFT potentials did not have off-diagonal matrix elements at momenta larger than 
about k 2 « 10 — 11 fm~ 2 . As counterpoint to this effect, Figure IBTBI shows the relative 
error plot, with a range of A's at E = 100 MeV and n exp = 8, using the Argonne Vi$ 
potential, which has matrix elements extending out to momenta of 30 fm - . Here 
one can see less saturation out to A = 10 as there are plenty of off-diagonal matrix 
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elements for the SRG to suppress between k 2 > 100 fm 2 and k 2 = 36 fm 2 , which 
corresponds to A = 6fm _1 . 
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Figure 3.7: Sample phase shift error plots using other choices of G s . On the left is 
G s = T 2 el and on the right is G s = H D for partial waves 1 S (top) and 3 Si (bottom). 



We repeated the phase shift calculations with other choices of the SRG generator, 
G s , that are diagonal in momentum space, including T 2 cl and Hd = T TC \ + Vb, where 
Vd is the (running) diagonal part of the bare potential. Samples of these results are 
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shown in Fig. 13.71 on the left (right) are plots using G s = T 2 (Hp) and on the top 
(bottom) are 1 Sq ( 3 5i). We found that these other choices for G s do not alter the 
qualitative features of the power-law behavior region of the previous error plots. This 
provides further evidence that the high- and low-energy decoupling results primarily 
from the partially diagonalized nature of the evolved potential. 

3.4 Decoupling and Deuteron Observables 

To test the generality of the observations made for phase shifts, the same de- 
coupling test techniques were applied to other low-energy observables such as the 
deuteron binding energy, radius, and quadrupole moment. The binding energy and 
momentum-space wavefunction were computed using standard eigenvalue methods. 
The computation of Qd and from the wavefunction uses [55] . 



oo 



Q, = ~ J a dk 



(3.14) 



and 



r d =- 



1/2 



(3.15) 



where u(k) and w(k) correspond to the S and D components of the deuteron wave- 
function respectively. We again computed relative errors in these observables and, 
as shown in Fig. 13. S\ the errors show the same behavior as observed for the phase 
shifts. That is, a power-law drop-off in the error begins at A just above A, with a 
slope determined by the sharpness of the regulator as given by n. 
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Momentum of Cut, A [fm ] Momentum of Cut, A [fm ] ] Momentum of Cut, A [fm '] 

Figure 3.8: The relative error vs. cut parameter A of the deuteron energy (left), rms 
radius (center), and quadrupole moment (left) of the deuteron with several values of 
the regulator parameter n indicated in the legends. In each case, the relative error is 
confirmed by the first-order perturbation theory shown in Eq. (13.131) . 

As with the phase shift, the analytic dependence of the error from cutting the 
potential can be estimated directly in perturbation theory. In this case, partial di- 
agonalization of the potential means that the deuteron wave function has negligible 
momentum components starting slightly above A. This in turn validates the expansion 
in Eq. (13.131) and the dependence of the errors on 1/A 2n . The numerical calculation 
of the error in perturbation theory is plotted in Fig. 13.81 and shows close agreement 
in the decoupling region A > A. 

3.5 Decoupling and Few-Body Energies with the NCSM 

The calculations described above have been only for two-particle systems. We 
can test whether the high-energy decoupling behavior extends to few-body systems 
by using No-Core Shell Model (NCSM) calculations of ground-state energies with the 
Many-Fermion Dynamics (MFD) code |56j . 

As noted before, the NCSM is a harmonic oscillator basis calculation of many- 
body systems where states can be explicitly enumerated. The size of the basis is 
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tracked by the parameter N max , the number of oscillators functions used in the basis. 
Calculations are variational in iV max , improving accuracy with larger bases, and in 
hu with a minimum occurring at the optimal value of hui for a given space and iV max . 
Here calculations are made for selected iV max and Kuj values. Also, note here that the 
NCSM is a black box calculational tool, which takes as input the harmonic oscillator 
matrix elements of potentials evolved and cut in momentum-space and provides NN- 
only calculations of few-body observables like the binding energies of the light nuclei. 
This is opposed to subsequent chapters where the NCSM machinery will be developed 
in order to access three-body interaction matrix elements. 

Conversion to the oscillator basis introduces a unique truncation of the interaction 
based on the values of iV max and hu (as detailed in Appendix [B]) . However, here the 
two-body Hamiltonian is being evolved in the momentum basis so that the evolved, 
and therefore decoupled, potential is unitarily equivalent to the initial A = 2 potential. 
We will take advantage of that decoupling to avoid the oscillator basis truncations 
at a given basis size, iV max , and we will use the same systematic cutting as in the 
previous sections to study the behavior of decoupling in few-body systems. 

In this section, only NN interactions have been considered, with the testing of 
decoupling for many-body forces deferred to later chapters. Because induced three- 
body (and higher) forces are not included, the converged values will differ significantly 
with A. In spite of this, we will still be able to see the evidence for NN decoupling in 
these systems. Also, note that the general features of the SRG exhibited in Section [3721 
implies that off-diagonal matrix elements of the three-body force will be suppressed 
as well, with decoupling as an expected consequence. The issue of three-body force 
evolution will be covered explicitly in chapters H] and 
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Figure 3.9: Calculations of the 4 He ground-state energy using the NCSM. On the left 
is the energy obtained from the NCSM for potentials evolved to several different A 
values as a function of the cut (regulator) momentum A with n = 8. On the right 
is the relative error of the energy for the A = 2 fm~ function of the cut 

momentum (with n = 8) for several different harmonic oscillator basis sizes. Also 
shown is the slope of the error in the decoupling region predicted from perturbation 
theory (dotted line). 



We first verified that the decoupling behavior already observed using a direct 
calculation of the deuteron wavefunction is reproduced using the MFD. We then 
calculated a series of larger nuclei, including 3 H, 4 He, and 6 Li, comparing results 
from uncut and a range of cut potentials evolved to different values of A. On the 
left panel of Fig. 13.9} the 4 He ground-state energy is plotted versus the regulator 
parameter A for several different values of the SRG flow parameter A. Each of the 
plotted points is at a basis size N max = 12 which is within several hundred keV of 
the energy from extrapolating to N max = oo. A similar plot for 6 Li is given in the left 
panel of Fig. 13.101 using a basis size iV max = 8, also within several hundred keV of the 
extrapolated energy for A = 1.5 fm" 1 but still several MeV off for A = 3.0 fm" 1 . 
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Figure 3.10: Calculations of the 6 Li ground-state energy using the NCSM. On the 
left is the energy obtained from the NCSM for potentials evolved to several different 
A values as a function of the cut (regulator) momentum A with n = 8. On the right 
is the relative error of the energy for the same A's as a function of the cut momentum 
for the same A values but with two values of n. 



Both examples show that when the potential is cut with A comparable to A or 
lower, the converged energy is significantly different from the asymptotic uncut value, 
while it approaches that value rapidly as A moves above A. This means that, for 
smaller A, more high momentum matrix elements can be discarded without a loss of 
accuracy. This decoupling explains the greatly improved convergence with basis size 
seen in the NCSM for corresponding A values [55] . As noted above, the uncut (A — > 
oo) energies vary for each A because the SRG evolution includes the NN interaction 
only; the closeness of the results for A = 2fm _1 and 3fm _1 for 4 He is coincidental 
(see Ref. [35] for further discussion about the running of the energies). 

The quantitative behavior of the relative error parallels that observed for two- 
body observables, as seen on the right panels of Figs. 13.91 and 13.101 In all cases, for 
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a fixed value of A the power decrease in the error starting with A slightly above A is 
clearly seen, even though there are fewer digits of precision in the NCSM results (so 
the relative error is in the range 10~ 6 -10~ 5 at best). The same perturbative residual 
coupling is seen for different basis sizes, with the slope given by the dependence 1/A 2n , 
although the onset of the decoupling region shifts to higher A until the calculation 
is near convergence (see Fig. 13.91) . Similar results are found for other nuclei and for 
other values of A. 

3.6 Block Diagonalization 

As demonstrated in sections 13.31 - 13.51 decoupling between low-energy and high- 
energy matrix elements is naturally achieved in a momentum basis by choosing a 
momentum-diagonal flow operator such as the kinetic energy T rc i or the diagonal 
of H s ; either drives the Hamiltonian toward band- diagonal form. Renormalization 
Group (RG) methods that evolve NN interactions with a sharp or smooth cutoff in 
relative momentum, known generically as V\ ow k, usually rely on equations based on 
the half-on-shell invariance of the two-nucleon T matrix [5*71 15"5] . These approaches 
achieve a block- diagonal form characterized by a cutoff Abd (see left plots in Figs. l3.lTl 
and 13. 12[) using a Lee-Suzuki type transformation which unitarily transforms matrix 
elements into a low-momentum P space and a high-momentum Q space. As usually 
implemented they set the high-momentum matrix elements to zero but this is not 
required. 

The SRG can also produce block-diagonal decoupling similar to the sharp Vjowfc 
form by choosing a block- diagonal flow operator [581 159] , 

G - = ( P Y QH.Q ) s ^ • < 316 > 
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with projection operators P and Q = 1 — P. In a partial-wave momentum representa- 
tion, P and Q are step functions defined by a sharp cutoff Abd on relative momenta. 
This choice for G s , which means that r] s is non-zero only where G s is zero, suppresses 
off-diagonal matrix elements such that the Hamiltonian approaches a block-diagonal 
form as s increases (or A decreases). 




Figure 3.11: Comparison of momentum-space V\ ow k (a) and SRG (b) block-diagonal 
potentials with Abd = 2fm _1 evolved from an N 3 LO 3 Si potential [19"] . 

The plots in Figs. 13.111 and 13.121 show a comparison between the potentials renor- 
malized via the Vjowfc procedure (left) and those evolved with the block-diagonal SRG 
of Eq. (13.161) (right). The initial potential used here is again the N 3 LO potential from 
Ref. [49J. Figure [3.111 uses two-dimensional contour plots as usually presented and 
Fig. 13.121 uses a three-dimensional surface plot for comparison and reference. Both 
plots use a value Abd = 2 fm" 1 with the SRG potential being evolved to A = 0.5 fm" 1 . 
The agreement between V\ ow k and SRG potentials for momenta below Abd is striking. 
While the color scale is a bit saturated in the former (because we tried to keep the 
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Figure 3.12: Comparison of momentum-space V\ ow k (a) and SRG (b) block-diagonal 
potentials with Abd = 2fm _1 evolved from an N 3 LO 3 Si potential [^9] . 



Q-space visible), the latter helps give perspective. A similar degree of universality 
is found in the other partial waves. Deriving an explicit connection between these 
approaches is the topic of an ongoing investigation. 

One difference between these procedures is that VJowfe can be formulated as an 
integral transformatio rfl while the SRG is currently only differential; here V\ ow k is 
implemented as a single transformation on the initial matrix, while the SRG flows 
to the desired amount of transformation. An integrated form of the SRG would be 
useful in understanding the exact correspondence between the two procedures, and 
to what extent V\ ow k is a subset of these block-diagonal SRG forms. 

If one considers a measure of the off-diagonal coupling of the Hamiltonian, 

Ti[{QH s P)\QH s P)} = Tt[PH s QH s P} > , (3.17) 
13 V\ ow k has been expressed in both integral and differential forms [57] , 
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then its derivative is easily evaluated by applying the SRG equation, Eq. (12.81) : 

^Tr [PH S QH S P] 

= Tt[P Vs Q(QH s QH s P - QH S PH S P)} 

+ Tr[(PH s PH s Q - PH s QH s Q)Q Vs P] 
= -2Tr[(Q Vs P)\Qr] s P)) ^ . (3.18) 

Thus, the off-diagonal QH S P block will decrease in general as s increases [581 159] . 

The evolution of the "off-diagonal" matrix elements (meaning those outside the 
PH S P and QH S Q blocks) can be roughly understood from the dominance of the 
kinetic energy on the diagonal. Let the indices p and q run over indices of the mo- 
mentum states in the P and Q spaces, respectively. Analogous to the approximation 
in sec. 13.21 we can replace PH S P and QH S Q by their eigenvalues E p and E q in the 
SRG equations, yielding EH] 

-^j-h pq w r] pq E q - E p r] pq = -(E p - E q ) r] pq (3.19) 

and 

Vpg ~ Eph pq — h pq E q = (E p — E q ) h pq . (3.20) 
Combining these two results, we have the evolution of any off-diagonal matrix element: 

i-h pq ^-{E p -E q fh pq . (3.21) 

In the NN case we can approximate the difference of eigenvalues by that for the 
relative kinetic energies, giving an explicit solution 

Us)« V(0)e" s(£p - e ' )2 (3.22) 
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with e p = p 2 /M. Thus the off-diagonal elements go to zero with the energy differences 
just like with the SRG with T rel ; one can see the width of order 1/y/s = A 2 in the k 2 
plots of the evolving potential in Figs. 13. 131 and l3~T4l While in principle the evolution 
to a sharp block-diagonal form means going to s = oo (A = 0), in practice we need 
only take s as large as needed to quantitatively achieve the decoupling implied by 
Eq. fl3~22|) . 



k' 2 (1m 2 ) k' 2 (fm 2 ) k' 2 (fm 2 ) k' 2 (frrf 2 ) 




Figure 3.13: Evolution of the 3 Si partial wave with a sharp block-diagonal flow equa- 
tion with Abd = 2fm _1 at A = 4, 3, 2, and lfm -1 . The initial N 3 LO potential is 
from Ref. The axes are in units of k 2 from 0-11 fm -2 . The color scale ranges 
from —0.5 to +0.5 fm as in Fig. 13.111 
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Figure 3.14: Same as Fig. 13. 131 but for the x Pi partial wave. 
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This is illustrated in Figs. 13.131 and 13.141 which show the evolution with A of two 
representative partial waves ( 3 Si and 1 Pi) starting with the N 3 LO potential from 
Ref. [15] . One can clearly see the transformation to a block-diagonal form. Notice 
the "neck" between the P and Q spaces. Its width, especially visible in the X P\ 
channel, is proportional to A 2 . Presumably, though needing further investigation, 
this finite and smooth interface between the two spaces helps control errors, due to 
the mixing of P and Q space wavef unctions, that have plagued traditional Lee-Suzuki 
transformations [60] . 




Figure 3.15: Phase shifts for the 3 Si partial wave from initial potentials N 3 LO and 
AV18 and the evolved sharp SRG block-diagonal potential with Abd = 2 fm~ at 
various A, in each case with the potential set identically to zero above Abd- 

Definitive tests of decoupling for NN observables are now possible for Viowfe po- 
tentials since the unitary transformation of the SRG guarantees that no physics is 
lost. For example, in Fig. 13.151 we show 3 Si phase shifts from an SRG sharp block 
diagonalization with Abd = 2 fm _1 for two different potentials. The phase shifts are 
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calculated with the potentials cut sharply at Abd- That is, the matrix elements of the 
potential are set to zero above that point. The improved decoupling as A decreases 
is evident in each case. By A = 1 fm" 1 in Fig. I3.15[ the unevolved and evolved curves 
are indistinguishable to the width of the line up to about 300 MeV. 
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Figure 3.16: Errors in the phase shift (two partial waves) at = 100 MeV and 
three deuteron observables for the evolved sharp SRG block- diagonal potential with 
Abd = 2 fm -1 for a range of A's and a regulator with n = 8. Two partial waves are 
shown. 
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In Fig. l3.16l we show a quantitative analysis of the decoupling as in section [3731 The 
figure shows the relative error of the phase shift at 100 MeV calculated with a potential 
that is cut off as before at a series of values A cut . We observe the same universal 
decoupling behavior: a shoulder indicating the start of a perturbative decoupling 
region, where the slope matches the power 2n fixed by the smooth regulator. 

Figure 13.171 shows a limiting factor with respect to the decoupling generated by 
the block-diagonal SRG. The onset of the shoulder in A cut decreases with A until 
it saturates for A near Abd, leaving the shoulder at A cut ~ Abd- Thus, as A — > 
the decoupling scale is set by the cutoff Abd- Of course, the P space inside Abd is 
not transformed (diagonally or otherwise) by the SRG and the potential cannot be 
decoupled below this scale. 

Taking the block diagonal idea further we can apply more general definitions of 
P and Q. To smooth out the cutoff, we can introduce a smooth regulator j\, which 
we take here to be an exponential form: 

/a(*) = e"( fc2 /A| D )" , (3.23) 

with n an integer. For V[ ow fc potentials, typical values used are n = 4 and n = 8 (the 
latter is considerably sharper but still numerically robust). By replacing H^ d with 

G s = f A HJ A + (1 - f A )H s (l - / A ) , (3.24) 

we get a smooth block-diagonal potential. 

Fig. 13.181 shows an example of the smooth block- diagonal SRG with Abd = 2 fm _1 
and n = 4. As we evolve the potential down toward A = 1, the overlap of the P 
and Q spaces becomes significant and the potential develops new structure instead of 
becoming simpler. This sort of behavior indicates that there is no further benefit to 
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Figure 3.17: Decoupling error plots using the sharp block-diagonal generator on two 
partial waves ( 3 Si and 1 P±) in two different initial potentials (AV18 and N 3 LO (500 
MeV)). 
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Figure 3.18: Evolution of the 3 Si partial wave with a smooth (n = 4) block-diagonal 
flow equation with Abd = 2.0 fm" 1 , starting with the N 3 LO potential from Ref. |49j . 
The flow parameter A is 3, 2, 1.5, and lfm -1 . The axes are in units of k 2 from 0-11 
fm -2 . The color scale ranges from —0.5 to +0.5 fm as in Fig. 13.111 



evolving in A very far below Abd; i n fact the decoupling worsens for A < Abd with a 
smooth regulator. 
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Figure 3.19: Phase shift errors from potentials evolved with the smooth block-diagonal 
SRG. The block parameter is A B d = 2. Here the regulator of Eq. (I3.23P uses n = 4. 
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Figure 13. 1 91 shows the phase shift relative error plots for the smooth block-diagonal 
SRG applied in two partial waves of the N 3 LO potential. Note that as A is run down 
below the block parameter Abd the error worsens, signaling a collision of high- and 
low-energy degrees of freedom due to the overlap of P and Q space wavefunctions as 
defined by the smooth block transformation. 
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Figure 3.20: Evolution of the 3 Si partial wave with a second-order exact block- 
diagonal flow of equation 13.251 with Abd = 2fm~ 1 at A = 4, 3, 2, and lfm -1 . The 
initial N 3 LO potential is from Ref. [H]. The axes are in units of k 2 from 0-11 fm~ 2 . 
The color scale ranges from —0.5 to +0.5 fm as in Fig. 13.111 

Another type of SRG that is second-order exact and yields similar block diago- 
nalization is defined by 

rj s = [T, PV S Q + QV S P] , (3.25) 

which can be implemented with P — > /a and Q — > (1 — f\), with f\ either sharp 
or smooth. Figure 13.201 shows an example of evolution of the N 3 LO potential from 
Ref. [19] using this choice. 

We can also consider bizarre choices for f\ in Eq. (13. 24p . such as defining it to 
be zero out to Ai owcr , then unity out to Abd, and then zero above that. This means 
that 1 — f\ defines both low and high-momentum blocks and the region that is driven 
to zero consists of several rectangles. Results for two partial waves starting from 
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Figure 3.21: Evolved SRG potentials starting from Argonne Vi 8 in the 1 So (left) and 
x Pi (right) partial waves to A = 1 fm -1 using a bizarre choice for G s (see text). 

the Argonne v±s potential are shown in Fig. 13.211 Despite the strange appearance, 
these remain unitary transformations of the original potential, with phase shifts and 
other NN observables the same as with the original potential. These choices provide 
a proof-of-principle that the decoupled regions can be tailored to the physics problem 
at hand. 

In the more conventional SRG, where we use G s = T re i, it is easy to see that 
the evolution of the two-body potential in the two-particle system can be carried 
over directly to the three-particle system. In particular, it follows that the three- 
body potential does not depend on disconnected two-body parts [2S1 HZ] • If we could 
implement rj s as proposed here with analogous properties, we would have a tractable 
method for generating V\ ow k three-body forces. While it seems possible to define 
Fock-space operators with projectors P and Q that will not have problems with 
disconnected parts, it is not yet clear whether full decoupling in the few-body space 
can be realized. In fact, if P and Q are A-dependent than the two-body potential 
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will not be completely fixed by evolution in the two-body space. This is an essential 
ingredient in the many-body calculations done in this thesis. 
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CHAPTER 4 



ONE-DIMENSIONAL MODEL 
IN THE HARMONIC OSCILLATOR BASIS 



As shown formally in Section [221 the SRG induces many-body forces as it evolves 
the Hamiltonian in an A-body space. For the SRG to be a useful tool, we must 
develop methods for calculating these induced many-body interactions and establish 
the conditions under which an initial hierarchy of many-body forces is maintained. 
In this chapter, we study a one-dimensional system of bosons as a proof-of-principle 
of a practical method to evolve and evaluate such forces, and establish a road map 
for full three-dimensional calculations carried out in chapter Beyond this, a one- 
dimensional system provides a straightforward laboratory that is indispensable in 
gaining deeper understanding of the SRG's effect on many-body potentials. 

Most previous applications of the SRG to nuclear structure have been in a momen- 
tum basis, where decoupling between low-energy and high-energy matrix elements is 
naturally achieved by choosing a momentum-diagonal flow operator such as the ki- 
netic energy T re \. However, evolution in a three-body space requires proper handling 
of the "dangerous delta functions" associated with the spectator particle to the two- 
body interaction. In a momentum representation, these delta functions are difficult 
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to resolve in a practical mesh and require solving separate flow equations for each 
combination of three-body matrix elements. 

Since we can evaluate the SRG flow equations Eq. (12. 8p with G s = T ve \ in any con- 
venient basis, an alternative is a discrete basis that allows direct application of the 
SRG flow equations in each A-body sector. We adopt this approach in the present 
chapter, using harmonic oscillator wave functions as our basis and mimicking the 
formalism used in the no-core shell model (NCSM) [TT1 [T2| [T3] to create properly 
symmetrized (for bosons) matrix elements in relative (Jacobi) coordinates. The re- 
striction to one dimension makes the construction particularly straightforward and 
requires only moderate matrix sizes. We use bosons in this model for easy comparison 
with existing model analysis. However, the boson ground states also coincide with 
fermion ground states when the flavor degeneracy is greater than the number of par- 
ticles, because the overall anti-symmetrization is realized by the flavor wavefunction. 

The choice of the harmonic oscillator basis is important because it is the only or- 
thogonal basis in which the center-of-mass can be factorized. With a translationally- 
invariant basis we have many fewer basis states and can achieve better convergence 
in a given A-bodj system. Other bases can be translationally invariant, like a sys- 
tem of correlated gaussians, but are non-orthogonal. Others, like Coupled-Cluster 
calculations are orthogonal but do not have a straightforward truncation in iV max . 

We use simple flavor-independent potentials that imitate the short-range repulsion 
and mid-range attraction characteristic of realistic local nuclear potentials. Previous 
studies of the SRG imply that properties of the transformations are primarily due to 
the matrix structure (G s , H s , choice of basis, etc), so we will be able to directly carry 
over some of our observations to three dimensions. Because the NCSM formalism is 
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already developed in three dimensions, we will see in chapter [5] that the generalization 
to three-dimensional fermionic calculations with spin-isospin degrees of freedom and 
using realistic nuclear interactions is algebraically straightforward, though far more 
computationally intensive. 

Another form of the NCSM uses oscillator wavefunctions associated with lab- 
frame coordinates and momenta instead of the relative Jacobi coordinates. This basis 
uses Slater determinants, making it much easier to antisymmetrize. This method 
is often known as "m-scheme" [10] as it is organized according to the m quantum 
number of the oscillator states. Similarly, a lab-frame NCSM scheme referred to 
as jj-coupled [M], takes advantage of the rotational symmetry of specific nuclei to 
reduce the size of the basis used in CC calculations. The major benefit of lab-frame 
coordinate based oscillator bases is the ability to generalize to arbitrary A due to its 
reliance on a Slater determinant to obtain an antisymmetrized basis as opposed to 
the recursive procedure used in Jacobi systems that builds on previous antisymmetric 
subsystems. On the downside, this results in poorer scaling of basis sizes for the A- 
body Hamiltonian matrices. As a check of the one- dimensional Jacobi coordinate 
oscillator basis we built a one-dimensional analog of the m-scheme NCSM, though it 
was unwieldy for our current purposes and not developed further. The details of this 
model are discussed in appendix [F] 

4.1 Symmetrized Jacobi Harmonic Oscillator Basis 

The No-Core Shell Model is a wavefunction-based method used to make ab initio 
calculations of nuclear bound-state properties. The basis it uses is built from the 
eigenfunctions of the harmonic oscillator Hamiltonian (described in Appendix |B]) up 
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to a maximum oscillator number, iV max . The functions are products of a gaussian and 
polynomials of order less than the value iV max . The other free parameter in defining 
this basis is cu, which sets the scale of the oscillator wavefunctions. This is usually 
multiplied by h, which is set to one in this unitless model. 

Each boson (fermion) wavefunction in an A-body system must be fully symmetric 
(antisymmetric). Therefore, we must identify and include all the symmetric states 
which the system may occupy. A general and systematic way to build this basis is to 
list all the possible states and build a symmetrization operator in this full basis. The 
eigenstates of this operator with eigenvalue one will be the exclusive set of symmetric 
states relevant for the system up to the stated iV max ; those with eigenvalue zero are 
irrelevant and can be thrown out. Also note that bases of smaller iV max are a subset 
of larger bases. The details involving the organization of states for this method is 
discussed in Appendix IC1 

The initial two-body potential, usually given in the relative momentum basis, can 
be converted to the relative oscillator basis simply by inserting a (semi-l 14 l) complete 
set of states to expand the momentum states in terms of oscillator wavefunctions. For 
three-body calculations, this interaction must be embedded in the three-body sym- 
metric basis obtained from the symmetrizer diagonalization. Calculations in higher- 
body bases require a general iterative procedure. The A-body basis, being symmetric 
already, is taken as an explicit set of states to which is added the additional Jacobi os- 
cillator coordinate. A new symmetrizer is built and the A + 1-body basis is obtained. 
The potential in the A-body basis is embedded in the new space. Diagonalization 
14 see Appendix IBl for a discussion on the incompleteness of the oscillator basis. 
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or evolution can proceed as usual. The details of this iterative process, including 
potential embedding and symmetrizer construction, are given in Appendix O 




Figure 4.1: A snapshot of the embedding process involving the initial (unevolved) two- 
body potential. The far left is a two-body, one-dimensional, potential in momentum 
space [62] with axes labelled by k 2 , the middle shows that potential converted to 
the two-body symmetric oscillator space, and the far right shows it embedded in the 
symmetric three-body basis. Each box is a matrix element in the oscillator basis, and 
the matrix elements are organized in order of increasing energy. The oscillator basis 
axes are unmarked due to the complicated organization of those bases (see text). For 
visibility, a small basis with iV max = 12 is shown. 



Figure I4TT1 shows an example of the process to embed a two-body momentum basis 
potential into the oscillator bases for A = 2, and 3. The left shows an initial two- 
body potential in the momentum basis. The center shows the same potential in the 
two-body oscillator basis and is physically identical to the momentum basis version 
except for truncation errors induced by oscillator basis conversion (see appendix IB"]) . 
The far right shows the potential embedded in the three-body basis. The axes here 
are unmarked because the order of the states is partially arbitrary. Here, as is usual, 
they are organized in order of increasing energy, labeled by Na (for an A-body state), 
up to the maximum iV max . Any degeneracies in Na are labeled by an arbitrary index 
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i A . Note the many matrix elements on the far right plot which are zero (green in 
this color scheme) indicating transitions which cannot be described by only two-body 
forces and a spectator particle. These matrix elements represent unique three-body 
transitions and, while some of these may remain zero due to symmetry properties, 
they are now available to the SRG during evolution. Any change in the value of these 
matrix elements represents three-body forces being induced during evolution in this 
basis. 

Once we have constructed a complete symmetrized basis for A particles (speci- 
fied by the value of iV max ) and evaluated the Hamiltonian matrix elements in this 
basis, applying the SRG transformations is immediate. That is, we have coupled, 
first-order differential equations for each matrix element of the Hamiltonian, with 
the right side of each equation given by a series of matrix multiplications. This is 
efficiently implemented in any computer language with matrix operations and differ- 
ential equation solvers. Just as in the momentum basis, individual matrix elements^ 
of the Hamiltonian obey the SRG's differential equations: 

^(N' A i' A \(V A ) s \N A i A ) = (N' A i' A \[[T,H s ],H s ]\N A i A ) 

= (N' A i' A \TH s H s \N A i A ) + (N' A i' A \H s H s T\N A i A ) 

-2(N' A i' A \H s TH s \N A i A ) . (4.1) 

We have defined dT/ds = so that all of the flow occurs in the matrix representation 
of the potential, (V A ) S . Using the matrix representations of T and H s in the \N A i A ) 
basis, the right side of Eq. (14 .11) is simply a series of matrix multiplications. The initial 

condition at s = is the initial Hamiltonian, (N A i' A \T + V A \N A i A ), which can have 

15 The ^4-body states are labeled by the total energy, N A and a quantum number, i A , introduced 
to keep track of states degenerate in N. In three dimensions this list will include angular momenta 
and isospin. See appendix [C] for the complete notation. 
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few-body components in Va- We consider here both two-body-only and two-body 
plus a three-body component (for A > 3). 




Figure 4.2: Same as in 14.11 but each matrix is now evolved to A = 2 in each basis 
shown. Again, iV max =12. 

Figure 14.21 shows an example of evolution in the oscillator bases for A = 2, and 
3. The left shows the evolution, to A = 2, of a two-body potential in the momentum 
basis. The center shows the evolution of the potential in the two-body oscillator 
basis and is physically identical to the momentum basis version again except for 
errors induced by oscillator basis conversion. The far right shows the evolution in the 
three-body basis and includes the evolution of two body forces but now also includes 
induced three-body forces not available in either of the two-body spaces. 

The SRG induces few-body forces as it evolves an initial interaction in a few- 
particle space. To study the contributions of different few-body forces we must isolate 
these components of the full interaction. The two-body force evolution keeps the A = 
2 binding energy invariant under evolution and defines the evolved two-body forces. 
We can isolate the three-body force from the two-body matrix elements by embedding 
the evolved two-body-only force in the three-particle space and subtracting it from 
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the full two-plus-three-body evolved interaction. Because we are in a symmetrized 
basis, this operation is relatively straightforward. However, we must be careful to 
keep track of the correct combinatoric factors associated with embedding A-body 
forces in higher spaces. In our MATLAB implementation these procedures take only 
a few lines of code in addition to the fully unitary calculation. 

The NCSM, both one- and three-dimensional, is a variational calculation in the 
two parameters, A max and frw. Largely this can be understood in terms of the 
truncation effects discussed in Appendix [Bl Essentially, conversion to an oscillator 
basis imposes both ultraviolet and infrared cutoffs Auv and Air determined by the 
parameters of the basis A max and hu. As derived in the appendix, increasing A max 
improves both cutoffs by raising Auv an d lowering Air, but changing hw shifts both 
cutoffs in the same direction. So we will be looking for a sufficiently large value of 
A max and the optimal value of fkv for a given A max . 

4.2 One Dimension Potential 

The bulk of the calculations in this chapter adopt a model from Ref. [52] that 
uses a sum of two gaussians to simulate repulsive short-range and attractive mid- 
range nucleon-nucleon two-body potentials: 

V W (x) = J^e-* 9 /^ + J^_ e -^| ( 4. 2) 

or 

V W(p,p') = J^ e -(^')M/8 + J^ e -( P - P ')V|/8 (4 3) 
2tt\/2 27T\/2 

The parameters used in Ref. [52] were chosen so that the one-dimensional saturation 
properties correspond to empirical three-dimensional properties, but we also want 
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Table 4.1: Parameter sets for the two-body potential of Eq. (14.31) . 



name 


Vy 


0"1 


v 2 


^2 


v a 


12. 


0.2 


-12. 


0.8 


v, 


0. 


0.0 


-2.0 


0.8 



to explore a range of parameters to test what behavior is general and what relies on 
specific features. Otherwise, the units for these model interactions are straightforward 
with the mass (same for all particles) and h taken to one, though they are displayed 
explicitly in some equations. We start with the parameters listed in Table 14.11 The 
potential V a is from Ref. [62] and is plotted in Fig. 14.31 We will fix the range of the 
attractive part and vary the relative strength and range of the repulsive parts and vice 
versa. We also vary the purely attractive potential Vp, which was used in Ref. [63] and 
is also plotted in Fig. 14.31 The eigenvalue problem for the relatively small matrices 
considered here can be solved by any conventional matrix diagonalization program 
(MATLAB was used here). 

To test that the symmetrized harmonic oscillator basis was correctly constructed 
for A = 2, 3, and 4 (see below), we first diagonalized the Hamiltonian using the purely 
attractive gaussian two-body potential Vp. The normalization is such that Vp(x) 
becomes a delta function with strength V 2 as o 2 — ► [63J (note the numerical factors 
from the Fourier transform because of our normalization of the Jacobi momenta). 
This limiting case has a known analytic solution for the (only) bound state of A 
bosons. For finite a 2 , we were able to confirm the accuracy of the diagonalizations 
as a function of the basis size iV max by comparison to coordinate-space stochastic 
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Figure 4.3: Potentials (dotted line, with axis on right) and probability distributions 
(other lines, with axis on left) for the lowest two-body bound state as a function of 
x — \x\ — %2\ at different stages in the SRG evolution (A = l/s 1//4 ). The left plot is 
V a and the right plot is Vp (see Table 14.11) . 



variational method (SVM) calculations using a published code [64] adapted to one 
dimension [65J [66] . 

We also explore the evolution of Hamiltonians with an initial three-body force. 
We choose a regulated contact interaction in the three-particle momentum space, 

V®(p, q,p', q') = c E f A (p, q)f A (p', q') , (4.4) 

where ce is the strength of the interaction and the form factor f\ depends on the 
Jacobi momenta as 

fA(p, q ) = e-^y^ n . (4.5) 

The regulator cutoff A sets the scale of the fall-off in momentum and n determines the 
sharpness of this fall-off. This form is analogous to the regulated three-body contact 
interactions used in chiral effective field theory [67] discussed in Appendix |X] We 
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Table 4.2: Ground-state energies for two-body potentials from Table fl~Tl with various 
strengths of the initial three-body potential Eqs. fl4.4|) - fl4.5p with A = 2 and n = 4 
for A = 2, 3, and 4. 





Ce 


E2 


E s 


Ei 


v a 


-0.10 


-0.920 


_3 


223 


-7.125 


v Q 


-0.05 


-0.920 


_2 


884 


-5.832 


v Q 


-0.01 


-0.920 


-2 


628 


-4.906 




0.00 


-0.920 


-2 


567 


-4.695 


v a 


0.01 


-0.920 


-2 


507 


-4.494 


v a 


0.05 


-0.920 


-2 


278 


-3.798 




0.10 


-0.920 


-2 


027 


-3.179 


v, 


-0.10 


-0.474 


-3 


379 


-8.412 


Vp 


-0.05 


-0.474 


-2 


283 


-5.727 


Vp 


-0.01 


-0.474 


-1 


792 


-4.183 


Vp 


0.00 


-0.474 


-1 


708 


-3.846 


Vp 


0.01 


-0.474 


-1 


626 


-3.517 


Vp 


0.05 


-0.474 


-1 


370 


-2.451 


Vp 


0.10 


-0.474 


-1 


240 


-1.874 



have not explored in detail the impact of adjusting A and n but have focused on how 
the SRG handles a varying strength eg. All results here are for A = 2 and n = 4. 

A sampling of ground-state energies are given in Table 14.21 A few simple features 
to note are the fact that the varied three-body force has no effect on the two-body 
binding energy and that a negative ce is added attraction and thus the 3- and 4-body 
bound states are deeper. For completeness we also include in Table 14.31 the optimal 
hu for each of these interaction choices. 

Most of the figures in this chapter show calculations with N max = 28. With this 
basis size, the ground-state energies are generally converged to one part in 10 4 , which 
is more than sufficient for our purposes. As usual, increasing N m&x leads to rapidly 
increasing matrix sizes and computation times; times for A = 3 with Af max = 32 are 
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Table 4.3: Optimal hws for the potentials in Table l4~2l with A = 2, 3, and 4. 





CE 


hLu(E 2 ) 


hw(E 3 ) 


tvulEA 


v a 


-0.10 


4.4 


5.8 


5.5 


v a 


-0.05 


4.4 


6.1 


5.8 


v a 


-0.01 


4.4 


5.1 


5.6 


v a 


0.00 


4.4 


5.0 


5.3 


v a 


0.01 


4.4 


4.9 


5.0 




0.05 


4.4 


4.7 


4.4 


v a 


0.10 


4.4 


4.5 


3.8 




-0.10 


1.1 


0.1 


0.1 


Vp 


-0.05 


1.1 


0.1 


0.1 




-0.01 


1.1 


0.5 


0.5 


v P 


0.00 


1.1 


1.6 


2.1 




0.01 


1.1 


0.5 


0.5 




0.05 


1.1 


0.2 


0.3 




0.10 


1.1 


0.1 


0.1 



a factor of 3 longer than with iV max = 28 and with iV max = 40 the time increases 
by another factor of 10. These timing and scaling properties of the one dimensional 
NCSM are explored in depth in Appendix [Dl 

4.3 Evolution of Many-Body Forces in Bound States 
4.3.1 Two-body Results 

We first consider the bound state of two identical bosons using the potential V a . 
Because the SRG is a series of unitary transformations, we expect that the binding 
energy will not be changed by evolving the two-body interaction in the two-particle 
space. Indeed, we find it to be constant to high accuracy. The ground-state wave 
function, however, changes dramatically, as seen from the probability densities plotted 
in Fig. 14.31 The initial probability density exhibits a sizable "wound" near the origin 
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that is filled in as A decreases. By A = 2 there is no signature of a repulsive core (and 
the wave function is modified out to larger x). This is the same pattern seen for the S- 
wave component of deuteron wave functions starting from three-dimensional nucleon- 
nucleon S-wave potentials with strong repulsive cores such as Argonne vi$, with the 
"uncorrelated" final wave function at A = 2 roughly comparable to A = 1.5 fm" 1 for 
the deuteron [3Tj . 




Figure 4.4: Even part of the SRG potential [V} (p,p') + V} (p, — p')] in dimensionless 
units as a function of p and p' for A = oo, 5, 3, and 2 (where A = 1/s 1 ' 4 ). The initial 
potential is V a from Table 14. 11 

The evolution of the potential in the momentum basis, shown as a color contour 
plot in Fig. 14.41 also demonstrates this behavior. (The even part of the potential is 
shown, which is the analog of the S-wave part.) The initial potential is dominated by 
strongly repulsive matrix elements coupling low and high momenta. The evolution 
in A band diagonalizes the potential to a width in p 2 of roughly A 2 while a soft 
attractive part emerges in the low- momentum region. The pattern in Fig. 14.41 reflects 
increasing non-locality as A is lowered, which in turn reduces the wound in the wave 
function. From the probability density and the momentum space plots we estimate 
that evolving to halfway between A = 2 and 3 for V a corresponds roughly to the A 
scale typically used in nuclear structure calculations (around 2fm _1 ). 
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4.3.2 Three-body Results 

To calculate properties of the three-particle system we construct the Hamiltonian 
in the basis of symmetric three-particle eigenstates as described in Sec. 10.41 The SRG 
evolution of the potential in the three-particle space leaves the ground state energy 
invariant if the full Hamiltonian is kept, because the transformations are unitary. 
However, the Hamiltonian matrix elements in this space do not follow simply from 
the pairwise sum of the two-body potential matrix elements; as the SRG evolves a 
three-body force is induced even if its initial strength is zero. 
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Figure 4.5: The lowest bound-state energy E 3 for a three-particle system as a function 
of A with the initial two-body-only potential V a . The (red) curves with squares include 
the full evolution of the Hamiltonian while the (black) curves with circles use the 
two-body potential evolved in the two-particle system. The right frame shows two 
additional results from varying o\ and Vi from the values in Table 14.11 



The effect of this full three-particle space SRG evolution is shown in Fig. 14.51 
for initial two-body potential V a and with initial = (non-zero values of eg 
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are considered in the next section). We plot the ground-state energy for the three- 
particle system both with the initial two-body interaction embedded in the three- 
particle symmetric space and then evolved (the red curve with squares) and also 
with the two-body interaction evolved in the two-particle space before embedding in 
the three-particle space at each A (the black curve with circles). We can see that 
the energy evaluated with the two-body interaction alone deviates noticeably as A 
drops below 10. This variation is the signature that the two-body transformation 
is only approximately unitary in the three-particle sector. The error reaches a peak 
in A between 2 and 3 and then decreases. The same pattern has been observed 
for NN potentials in three dimensions [35] and remains qualitatively the same when 
parameters in the potential are varied (e.g., see the right plot in Fig. 14.51) . We made 
the same calculation using the purely attractive initial two-body potential Vp, which is 
shown in Fig. 14.61 Here the induced three-body force has the opposite sign and there 
is no maximum, which implies that the qualitative pattern of evolution is dictated 
by the interplay between attractive long-range and repulsive short-range parts of the 
potential. These features are explored further in Sect. 14.41 

4.3.3 Results for A = 4 and A = 5 

Next we turn to A = 4 and A = 5, where we expect to see the effects of induced 
three-, four- and five-body forces. The key issue is the relative sizes of these contri- 
butions; we are looking to test whether an initial hierarchy of few-body interactions 
is preserved and therefore can be truncated with a controlled error. In Chapter we 
will address the same question in the realistic NCSM. 
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Figure 4.6: Same as Fig. 14.51 but with initial potential Vp. The right frame shows two 
additional results from varying and V2 from the values in Table 14.11 



In applying the SRG in the four-particle space we have three different calculations 
of the ground-state energy to compare. The first is the two-body potential embedded 
successively in the three- and four-particle spaces and then evolved in the four-particle 
space. The resulting unitary transformations will induce three- and four-body inter- 
actions that leave the eigenvalues invariant. We can also evolve in the two-particle 
space before embedding in the three- and four-particle spaces and diagonalizing. As 
we saw before in Fig. 14.51 and see now in Fig. 14. 7\ the two-body-only evolution de- 
viates because the Hamiltonian is not evolved by an exactly unitary transformation. 
Finally, to find the relative size of the three and four-body interactions we can evolve 
in the three-particle space, thereby inducing only three-body forces. Note that the 
two and three-body forces must be embedded differently in the four particle space 
because they have different combinatoric factors associated with them, i.e., there are 
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Figure 4.7: The lowest bound-state energy E4 for a four-particle system as a function 
of A with the initial two-body potential V a and different initial three-body force 
strengths (ce = ±0.05). 



(2) = ^ P a hs and ( 3 ) = 4 triplets. So, the proper mixture of two- and three-body 
force contributions to the four-particle system interaction is V = 6V^ + AV^. 

All three of these calculations for A = 4 are shown in Fig. 14.71 for the two-body 
potential V a and several choices of the initial three-body force. The magnitude of eg 
was chosen so that the fractions of the A = 3 and A = 4 ground-state energies from 
the three-body interaction are roughly comparable to the corresponding fractions for 
nuclei using typical realistic NN potentials. The qualitative behavior is similar for 
other choices of ce and V^. In all plots the curves for the two-body-only (black 
line with circles), the two-plus-three (red line with squares), and the full two-plus- 
three-plus-four interaction (blue line with diamonds) show the hierarchy of different 
few-body interaction components. The difference between the square and diamond 
lines represents the contribution of the four-body force, and the difference between 
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the circle and square lines is the contribution of the three-body force alone. In the left 
and right panels, where the initial three-body force has been switched on, we can see 
that it is perturbatively small from the balanced attractive/repulsive contributions. 

The four-body contribution is at most ten percent that of the three-body, which 
is itself small compared to the two-body contribution except when the latter gets 
small for small A (note the expanded scales on the figures). Considering calculations 
with different ce values, we see that the A dependence of the induced four-body part 
depends on the interplay of initial and induced forces. In some cases noticeable (but 
small) evolution starts at A = 10 while in other cases it is deferred until much smaller 
A. Regardless of the details, we stress that there is no sign that induced many-body 
forces have rapid growth with A or exhibit unusual scaling. 

We repeated for A = 4 our test of decoupling that was shown in Fig. l4.211 for A = 3. 
A similar pattern of decoupling is found, namely an increased degree of decoupling 
until a A corresponding to the minimum of the two-body-only ground-state energy of 
the A = 4 system, after which it deteriorates. 

In Fig. 14.81 we show results for the SRG evolution, with initial potential V a and 
no initial three-body interactions, in a five-particle system for several values of iV max . 
The right panel, with iV max = 28, shows the converged result. We see a decreasing 
hierarchy of induced many-body contributions for all A; the five-body contribution 
is essentially negligible (or not distinguishable from numerical noise). Differences at 
the lower iV max sized spaces arise both because the space needs to be large enough 
for convergence to the exact energy eigenvalues but also because the initial evolution 
of the potential needs a sufficiently large space. Decoupling may improve this feature 
but is dependent on the type of SRG used and the basis in which it is implemented. 
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Figure 4.8: The lowest bound-state energy E 5 for a 5-particle system as a function of 
A with an initial two-body-only V a potential for several values of iV max . 

4.4 Diagrammatic Analysis of Many-Body Force Evolution 

In Ref . |47] , a diagrammatic approach to the SRG equation was introduced, which 
organized the independent evolution of two- and three-body (and higher-body) poten- 
tials. This formalism is necessary in a momentum basis to avoid "dangerous" delta 
functions from spectator particles. In this section we examine how induced three- 
body interactions evolve in our one-dimensional laboratory and make connections to 
the diagrammatic expansion. 

As raw material for this analysis, we plot in Fig. 14.91 the error of the evolving two- 
body-only binding energy while varying several parameters of the initial two-body 
interaction V a for A = 3 ground-state energies. In the left plots we vary the range 
(top) and the strength (bottom) of the attractive part of V a . In the right plots we 
vary the range (top) and the strength (bottom) of the repulsive part. In Fig. 14.101 we 
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Figure 4.9: Differences of two-body-only and two-plus-three-body A = 3 ground-state 
energies as a function of A. Each of the parameters of the potential V a are varied 
in each plot as the other parameters are held constant. The upper panels vary the 
ranges while the lower vary the strengths; the left vary the attractive part and the 
right vary the repulsive part. 
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present similar plots for a simpler system that starts with the initial attraction-only 
two-body potential Vp. Note again that the plots here show the induced error between 
the two-body-only and two-plus-three-body binding energies. 

Certain qualitative features are found as expected in these figures. Shorter ranges 
imply enhanced coupling from low to high momentum and therefore we anticipate 
that the evolution will start sooner (i.e., at higher A). This is seen clearly on the 
left in Fig. 14.101 and for the variation of the (shorter-ranged) repulsive potential in 
Fig. 14.91 (top right plot). There is also an unsurprising increase in the magnitude of the 
induced three-body interaction at each A with increased magnitude of the potential, 
as seen on the right in Fig. 14.101 and on the bottom left in Fig. 14.91 For a more 
definitive analysis we need to recall the discussion from Ref. |47j . 

The SRG evolution equation for the three-particle sector in the notation of Ref. [UJ 

is 

^ + ^ = ^^ + \vt\v^) + \vf\v^) + ^\v^) + \vf\m) 
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ds 






dVs 
ds 







+ 



Figure 4.11: A diagrammatic decomposition of the SRG Eq. (14.61) . A circle at a vertex 
denotes a commutator with T re i. 



where each bar denotes a commutator with T rc \. We remind the reader that dT re \/ds = 
by construction. A diagrammatic decomposition of this equation is shown in 
Fig. 14. 1 11 In the two-body sector, the equation reduces to the first term on the 
left and the first and third terms on the right (the first row in Fig. 14.11]) . These terms 
keep two-particle energy eigenvalues invariant under evolution. In the three-particle 
sector, Eq. (14. 6 p results in not only these two-body graphs with a disconnected spec- 
tator but additional graphs involving connected combinations of two and three-body 
interactions. The diagrams with two-body interactions and a disconnected spectator 
line satisfy the two-body evolution equations, and so will cancel out of the full three- 
particle-sector evolution equation. Thus the evolution of the three-body interaction 
is dictated by the connected diagrams (the second row in Fig. 14.11"]) . In summary, the 
evolution of the A-body potential in the A-particle system is given by 

dV^ 



ds 



[Va, H : 



s A 



(4.7) 



where the "A" subscript on the right side means the fully connected A-particle terms. 
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To make a connection between the individual terms in the three-body interaction 
evolution and the running of the ground-state energy, we need to derive the evolution 

(3) 

equations for the expectation value of Vy in the ground state. Denoting the ground- 
state wave function for the A-particle system by \ipf), it evolves according to 

l#> = *Wio>, = r,.\tf) , (4.8) 

where U s is the SRG unitary transformation at s and 

Vs = ^fU} = - V l . (4.9) 
Then the matrix element of an operator O s evolves according to 

^{tf\o.\tf) = (^\^r - [v.,o.]\tf) ■ (4-io) 

If the operator O s evolves according to O s = U s O s= oU}, then the matrix element 
vanishes, as when O s = H s . 

However, if we wish to see how one part of H s evolves, such as the expectation 
value of we obtain 

£<1>}\v}*>\tf) = (1>t\^- ~ h,vPM?) , (4-11) 

which does not give zero in general because V s ^ ^ U s v}= Ul. In the two-particle 
case, the analog of Eq. (14.111) gives d(V^)/ds = ([r) s , T re i]). In the three-particle 
case, we can expand Eq. (14.111) as 

^IK (3) I^> = (MVs,H s ] 3 -[rj s ,V®M s ) 

= (^|[Fi 3) ,T rel ] + [vfM%+ Fl 2) ,K (3) ] + WfM 2) ] + Ff } ,K {3) ] 

-[vf\vp}-[vf\v^m 

= (Mvf\n s ] + [Vf\v}\-[vf\v}^s) 

= (Mvf\v}%-[vf\v®m, (4.12) 
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where V s and V s are the commutators V s = [T rel , V s ] and V s = [T rol , V} ]. In 

(3) 

the third line, the expectation value of the commutator, [V s ,H S \, vanishes identi- 
cally. 

The term with the subscript "c" has had the two-body disconnected diagrams 
removed. In our MATLAB implementation, this subtraction is achieved by first 
embedding the two-particle-space evolved version of this commutator in the three- 
particle space. Computing [V s , V^} in the two-particle space alone involves only 
the one-loop two-body interactions, so embedding in the three-particle sector results 
in only the disconnected parts. This disconnected part can then be subtracted from 
the total three-particle sector version of the same commutator, leaving only the three- 
particle fully connected part. 

It is most useful for our analysis to convert from derivatives with respect to s to 
derivatives with respect to A using 4- = — ^zk- In Fig. 14.121 we show the ground- 
state expectation values of the right side of Eq. (14. lip , which are broken down into 
the two terms from the right side of Eq. (14.121) for A = 3 and various potentials. It 
is apparent that the drivers of three-body matrix element evolution depend on the 
interplay between long- and short-range, attractive and repulsive parts. The lower 
right panel of Fig. 14.121 shows an increasing attractive strength of the three-body 
force when starting from an attractive-only two-body potential. In this regime the 
dominant contribution to the evolution of the three-body potential matrix element is 
the tree-level two-body connected part [V s , V s ] c - This observation accounts for the 
behavior in the right graph of Fig. I4.10[ where the size of the error scales (roughly) 
like (V^ 1 ) 2 . Varying the long-range attraction strength in Fig. 14.91 shows a similar 
effect. 
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Figure 4.12: Contributions from individual terms to the A = 3 ground-state expecta- 
tion value d{V^)/d\ for several different two- and three-body potentials, as indicated 
in the plots. We emphasize that A < 2 is very small, comparable to A < 1.5 fm _1 for 
NN forces in analogous calculations with the NCSM [55] . 



More generally, the impact on ground-state energies of the induced three-body 
interaction depends on the details of the correlations in the wave function. The other 
plots in Fig. 14.121 for the more realistic initial two-body potential, V a , show the in- 
terplay between two- and three-body contributions to the three-body matrix element 
evolution. The three-body contribution to the three-body evolution stays small until 
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the longer-range attractive part of the potential begins to affect the evolution. Most 
of the change is from A = 8 to A = 3, which is dominated by (V , V^). Thus, the 
feedback of the three-body potential depends on the initial conditions, but in general 
insures that the binding energy contribution stays small. 
We can repeat the above analysis for A = 4 

±mvpm = (4%vf\v(%+[vf\vp] c 

+ F; 3) ,K (3) ]e-Fi 4) ,K (4) ]l^ 4) ) (4.13) 

where we find no fully connected terms with only two-body forces. Again, discon- 
nected terms involving two and three body potentials cancel out in the lower sectors. 
The leading terms are commutators with one V s y ' and one Vs , followed by con- 
nected terms quadratic in and one term quadratic in Vs A \ All terms are small 
and additional cancellations among them further suppress the four-body contribution. 
Thus, the initial hierarchy of many-body forces implies that induced four-body (and 
higher-body) forces will be small. 

In Fig. 14.131 we plot these contributions to the evolution of the four-body expec- 
tation value. Again it is more useful to convert the derivatives in s to derivatives in 
A. The interplay of contributions is much more complicated here. We can see the 
cancellations between the two terms involving V s and V s , and the term quadratic 
in V s ^ . The total derivative of Vj^ is kept small, along with the in term quadratic in 
Vs'' 4 '*, until below A = 2. No terms quadratic in V s ^ appear as no connected diagrams 
can be constructed for the four-particle evolution. This keeps the evolution of V s ^ 
low because they are one more order removed from the initial interaction. 

Also, note the differences in combinatoric factors when comparing this analysis in 
different sectors. Between the A = 3 and 4 sectors we can only consider the [V^\ Vs 3 ^} 
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Figure 4.13: Contributions from individual terms to the A = 4 ground-state expecta- 
tion value d(yl )/d\ for several different initial three-body potentials, as indicated 
in the plots. We emphasize that A < 2 is very small, comparable to A < 1.5 fm _ for 
NN forces in analogous calculations with the NCSM [35] . 



contribution. Diagrammatically speaking, in the A = 3 sector, this can only include 
vf^ and Vs 3 ^ fully contracted with each other. However, in the A = 4 sector the 
connected diagrams are only those contractions__with one loop. There are 3! = 6 
ways to choose the pair legs composing the loop 
6 times larger in four-body sector. This does not negate the effects of the absence 
of two-body contributions since they would also have large associated combinatoric 
factors. 

16 There are 6 ways to choose the legs on the receiving term of the diagram but this is true in both 
sectors and this 6 cancels out of the comparison. 



So the |Vi 3 \ V^l contribution is 
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4.5 Fitting Three-Body Forces 

Figures 14.141 and 14.151 show a simple fitting procedure for the three-body force. 
In this calculation the initial interaction was evolved only in the A = 2 basis and 
then embedded up to the A = 4 basis for calculations of the three- and four-particle 
binding energies. This corresponds to the black curves (circles) in the figures, which 
are all identical. Before, the interaction was also evolved in the A = 3 basis to obtain 
a renormalization that is unitary for A = 3 and used to compute the induced four- 



body forces, 
from xEFT, 



3ere an independent operator, inspired by a three-nucleon contact term 







is used to fit the missing three-body interaction to the three-body 
binding energy: 

yg) = D oe -l(pl+p 2 2)/^} n e ~{(p?+p?)/^} n (4.i4) 

where pi and P2 are the usual Jacobi momenta, the two regulator parameters, A 
and n e xp, set the cutoff and sharpness of the regulator and are kept constant over the 
range of A in each panel shown, and the strength, Do, was fit at each value of A. Then 
this fitted term plus the evolved two-body only interaction is embedded to A = 4 to 
obtain the four-body binding energy. This corresponds to the red lines (squares). 

As can be seen in the figures, this fitting procedure reduced the error due to 
missing three-body forces by 80 to 90 percent. In Fig. 14. 141 the regulators sharpness 
exponent, n cxp , is held constant and the range, A, is varied. A significant dependence 
on A is observed indicating significant momentum dependent structure in the evolved 
two-body force. In Fig. 14.151 now the range is held fixed but sharpness is varied. 
Here very little dependence is found indicating less sensitivity to the exact shape of 
the regulator. 

17 In the xEFTparadigm, any regularized delta function will work for the fitting performed here. 

98 




Figure 4.14: Ground-state energy of the A = 4 system with NN interaction evolved 
and then a three-body interaction term (with an exponential regulator) fit to the 
A = 3 ground-state energy The calculation is repeated with a variation of the range 
parameter, A 



-4.5 




> 2 - body only 
i 2+3flt - only 
unitary 



Figure 4.15: Same as Fig. 14. 141 but at one value of A and a range over the regulator's 
sharpness parameter, n exp 



The fits shown in Figs. 14.141 and 14.151 show that this one term can accommodate 
much of the missing three-body force. A natural progression would suggest adding a 
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second term with gradients as in a Taylor expansion: 

yg = D 2 e-^ + ^ A ^[pl +pl +p[ 2 +^ ]e -[( P ?M 2 )/AT (415) 

where D 2 is the new low-energy constant that also must be fit to data. In this one- 
dimensional model, what we choose as data is arbitrary. A particular fitting strategy 
used for xEFT's uses the average of 3 H and 3 He binding energies and the triton (3 
decay. In an attempt to be consistent with that work, we have chosen to fit our two 
parameters to the three-body binding energy (as we did with the one term fitting 
above) and the first excited state in the three-body spectrum as another piece of 
three-body data. 
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Figure 4.16: Same as Fig. 14. 141 but now the A = 3 ground-state energy and a second 
piece of data, the first excited state for A = 3, are fit to two three-body interaction 
terms. The first term is the same as in Fig. 14.141 and the second includes a gradient 
correction. 

Figure 14". 1 61 shows the predictions for the A = 4 binding energy after fitting the two 
terms to two different three-body quantities. As in Fig. 14. 14} n cxp , is held constant 
and the range, A, is varied. Now we have an expansion in the strength, D, which 
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we fit to three-body quantities. The resulting fit three-body interaction can then 
be used to make predictions of four- and higher-body observables. Here we display 
the prediction of the A = 4 binding energy after fitting both Do and D2. The 
improvement is inconclusive. This is an ongoing investigation and we must next 
check several features of the system and procedure to determine how much correction 
we should expect. We may need to choose a different second value for fitting if the 
one chosen now is already well fit by the first term. A logical next step is to reproduce 
the fitting approach previously employed in realistic calculations by varying A as a 
function of A. 

Fitting three-body interactions in this manner is another tool to study the form 
of induced many-body forces. It is complimentary to the vacuum expectation value 
analysis done in section 14.44 providing another quantitative tool to understand the 
SRG's behavior. In addition, well controlled fitting procedures will be indispensable in 
estimating the effects of some induced A-bodj forces that cannot be induced explicitly 
due to the size of a given A-body basis. 

4.6 Evolving Individual Operators: A Study of Cutoff Be- 
havior 

An aspect of the SRG which contributes to its versatility and power is the ability 
to evolve individual operators as well as the Hamiltonian. The majority of work in this 
area is contained in Ref . O H2] , but some preliminary work was done for few-body 
systems using a simple operator for the momentum distribution in the wavefunction 
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of few-body systems |68j : 

.4 



n (k) = 5^ / ^a(^i^2 • • • kA)S(k — kj)S(y^ j kj)dk\dki . . . cUca 
i=i J j=i 

A 

= A 2p A (kik 2 ■ ■ .kA)8(k — kA)8C^^kj)dkidk 2 . ■ -dkA 
J 3=1 

= A J ^\[p\V2 ■ . .pA-i)8(y^-^-(p ~ PA-i))dpidp 2 ■ ■ ■ dp A -i 

(4.16) 

where A is the number of particles, fc's are lab-frame momenta and p's are Jacobi 
momenta. The second line is a simplification due to the symmetry properties of the 
present basis. The wavefunctions are symmetric with respect to all the particles so 
the total is simply A times the momentum distribution of one particle. We choose 
the A'th particle (the last one) due to a special property of the last Jacobi coordinate 
that it can be related directly to the last single-particle coordinate when the center of 
mass momentum is zero. Using that property, the third line is converted from single 
particle coordinates to the Jacobi coordinates of the basis. 

Figure 14.171 shows the evolution of the momentum distribution within the three- 
body wavefunction for a series of iV max 's. In these plots the operator itself is not 
actually evolved. The expectation value of the evolved operator on the evolved wave- 
functions is simply a verification of the unitarity of U s , with the curves for all A's on 
top of the A = oo curve. Here the interesting part is using the bare operator with the 
evolved wavefunctions (or vice versa) to compute the momentum distribution. The 
horizontal axis is in single particle momenta while the cutoffs are in terms of Jacobi 
momenta. A conversion factor of ps \f~A is needed when analyzing the figures. 
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Figure 4.17: A sequence of plots showing the evolving components of the 3-body 
wavefunction for increasing iV max . Each of these plots uses hw = 5. 



The three panels in Fig. 14.171 show convergence as A max is increased. The series of 
curves in each plot shows that the evolving wavefunction has a considerable decrease 
in high-momentum components as A is lowered. At large iV max , the oscillator basis 
ultraviolet cutoff (Auv) is large and sufficient high-momentum information is included. 
Hence the convergence in the unevolved wavef unctions. At larger N max , the evolved 
wavefunctions are more converged than the unevolved because the high-momentum 
components have been suppressed and therefore the wavefunction is unaltered by Auv- 
One can see A uv explicitly in the unevolved wavefunctions using Auv ~ V ^ A max 
with Huj — 5. At A max = 12 we get Auv ~ 8 fm _1 which is about the location of the 
right-most shoulder of the unevolved wavefunction. Again at A max = 28 the cutoff 
is Auv ~ 12fm _1 corresponding to that shoulder. The shoulder in the momentum 
distribution is the clearest high-momentum structure to compare between the panels, 
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hence in this instance and at constant A, we have ignored the Jacobi conversion 
factor. 

As iV max is lowered we can see the infrared cutoff (Air ~ •» / J™ ) become signifi- 

v y ''max , 

cant. The oscillatory behavior in the JV max = 12 panel is the ringing from the Air cut- 
off inherent in the small oscillator basis. Given a wavefunction in position space (recall 
coordinate space plots of the wavefunction in Fig. 14. 3[) which is cutoff at some point 
x = 27t/Air will show ringing artifacts from the Fourier transform of the sharp cut. 
For a square function the minima should be Air apart. Due to the single-to- Jacobi 
conversion factor, the minima of the ringing should be a/AAir rs a/ (A)hw/N max 
apart from each other. 




Figure 4.18: A sequence of plots showing the evolving components of the 3-body 
wavefunction for fixed A r max and varied f\w. 



Figures f4. 181 and 14. 191 further explore the behavior of the cutoffs, Air and Auv- In 
the former, all the panels are at the same relatively high N max = 28 with hw varied. 
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Figure 4.19: A sequence of plots showing the evolving components of the 3-body 
wavefunction for fixed hui and varied N max . 

The behavior of the ultraviolet cutoff is evident with the right-most shoulder given by 
y/huN max . The infrared is clearer here, with the spacing between oscillations given 
by f=s a/ AHlu I N max . This is also evident in Fig. 14. 191 where iV max is varied with a fixed 
fko. The spacing decreases by a factor \f2 from N m£ix = 6 — > 12 or iV max = 8 — > 16 
This is surprisingly accurate in the regions with relatively uniform oscillation spacings 
at midrange momenta. This is indicative of how much high-momentum details have 
been transformed by the SRG. We are left with a smooth gaussian-like wavefunction 
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which is chopped at large distance by the infrared cutoff, leaving straightforward 
ringing effects. 




k k k 

Figure 4.20: A sequence of plots at high iV max showing the evolving components of 
the many-body wavefunction for A = 2,3, and 4. 



The momentum distribution plots of Figure 14.201 shows the same quantity as the 
previous plots but at a single large 7V" max and for different numbers of particles, 
A = 2,3,4. The same renormalization pattern for high-momentum components is 
evident regardless of the considered one-dimensional system. 

4.7 Convergence in the Oscillator Basis 

In Fig. 14.211 we test SRG decoupling [52J within the harmonic oscillator basis. 
The left panel shows the results when only the two-body evolved potential is used to 
calculate A = 3 binding energies. For the right plot, the potential was evolved in the 
full three-body space, allowing three-body forces to be induced. The Hamiltonians for 
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Figure 4.21: Decoupling in the three-particle system using the choice G s = T ie \. The 
initial V a potential is evolved to each A shown in a basis with iV max = 40. On the 
left only the two-body potential is evolved and embedded while the right involves the 
full unitarily transformed potential. Matrix elements of the potential are set to zero 
if one or both states have N > N cut and the resulting Hamiltonian is diagonalized to 
obtain the ground-state energies plotted. 



selected A values are diagonalized in bases of decreasing size, as measured by "N cut " , 
which is the cut-off applied to the potential to study its decoupling properties. The 
potential is set to zero for matrix elements for which one or both states has N > N cut . 
The degree of decoupling is measured by the point of departure, as N cut is lowered, 
from the asymptotic energy for iV max = 40. As the potential is evolved from the 
initial potential (A = oo) down to A = 2.5, decoupling is achieved for smaller spaces, 
which means convergence is reached for smaller basis sizes. This is the same pattern 
as found for realistic NN potentials In the left panel of Fig. I4.21[ the ground- 
state energies asymptote to different values because the evolution is not completely 
unitary. One might imagine that this would affect the decoupling, but we show this 
is not the case here. In the right panel, the induced three-body interaction is kept, 
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so the curves asymptote to the same energy at large N max , while the same pattern of 
decoupling is observed. 

We repeated for A = 4 our test of decoupling that was shown in Fig. l4.211 for A = 3. 
A similar pattern of decoupling is found, namely an increased degree of decoupling 
until a A corresponding to the minimum of the two-body-only ground-state energy of 
the A = 4 system, after which it deteriorates. 




N cut N cut 

max max 

Figure 4.22: Decoupling in the four-particle system using G s = T re \. The initial V a 
potential is evolved to each A shown in a basis with iV max = 40. On the left only the 
two-body potential is evolved and embedded while the right involves the full unitarily 
transformed potential. Matrix elements of the potential are set to zero if one or both 
states have N > N cut and the resulting Hamiltonian is diagonalized to obtain the 
ground-state energies plotted. 



We note that the decoupling benefits afforded by evolution in the oscillator basis 
are less straightforward than in the two-body momentum basis studied for NN in 
[52] . In particular the cut-off errors increase for A smaller than the point at which 
the two-body-only binding energy is at a minimum (i.e., for A = 1.5 in Fig. I4.2ip . 
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In the oscillator basis the Hamiltonian is not being driven towards a diagonal form 
by the choice G = T Te \, since it is not diagonal in this basis. In addition, the cuts 
made in the oscillator basis do not correspond to analogous cutting done to study 
decoupling in the momentum representation. However, control over cutting in the 
oscillator basis is the ultimate goal of the improved convergence in calculations with 
SRG evolution. Therefore, when working in the oscillator basis, one should look for 
a choice of G which optimizes the decoupling in that basis. 
One possibility is the harmonic oscillator Hamiltonian, ifho, 

H ho = T rcl + -mu 2 r 2 (4.17) 

which in the oscillator basis is simply a diagonal matrix with each state's value equal 
to Nfkjj for that state. Such a matrix is easy to build from the lists of state quantum 
numbers readily available from the oscillator basis code. This choice of G should drive 
the Hamiltonian to a diagonal form in the oscillator basis. 

Figures 14.231 and 14.241 shows the decoupling analogous to Figs. 14.211 and 14. 22} but 
now using G s = H^. Here the SRG drives the Hamiltonian towards a form diagonal 
in iV max (and hence, energy), which is block diagonal in the oscillator basis. Again 
the calculation on the right has the full three-body (or four) Hamiltonian evolved 
before being cut to measure the convergence. Here, with G = H^ Q evolution, there is 
no saturation effect in the increased convergence. Instead the convergence continues 
to steadily improve with decreasing A. However, the plot on the left side of the figure 
has the two-body (two+three-body) Hamiltonian evolved before being embedded to 
the A = 3 (A = 4) space and cut for study. Here, as in the analogous T re \ evolution 
plot, we expect to see the asymptotic values of the different A curves vary according 
to the missing induced three-body (four-body) forces. But in the H^ evolution, the 
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Figure 4.23: Decoupling in the three-particle system with G s = if ho- The initial V a 
potential is evolved to each A shown in a basis with iV max = 40. On the left only the 
two-body potential is evolved and embedded while the right involves the full unitarily 
transformed potential. Matrix elements of the potential are set to zero if one or both 
states have N > N cut and the resulting Hamiltonian is diagonalized to obtain the 
ground-state energies plotted. 



variance is dramatically larger, due to spurious bound states produced during the 
evolution of two-body (two+three-body) only matrix elements. These spurious states 
are a problem for the use of ifho as the SRG generator and they are further explored 
in Appendix [El 

In addition to the effects of the generator itself on the Hamiltonian, the method 
of cutting used here was applied directly from that used to study decoupling in the 
momentum basis. While cutting in iV max is a variational technique it is not an 
exclusive test of high-low energy decoupling because a cut made in the oscillator 
basis imposes both IR and UV cutoffs while a cut in momentum space is simply a UV 
cutoff. Therefore the decoupling studies performed here in the oscillator basis cannot 
be directly connected to the decoupling studies of chapter [3J However, cutting in 
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Figure 4.24: Same as for Fig. 14.231 but for the A = 4 binding energy. On the left only 
the two+three-body potential is evolved and embedded while the right involves the 
full unitarily transformed potential. 

-^Vmax does establish the convergence with basis size and is a direct measure of success 
in the broader program of controlling resource requirements. The issue of oscillator 
space decoupling is further informed by the discussion on general oscillator basis 
features in appendix iBl 
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CHAPTER 5 



EVOLUTION OF MANY-BODY FORCES IN THE 
NO-CORE SHELL MODEL (NCSM) 

Our calculations are performed in the Jacobi coordinate harmonic oscillator (HO) 
basis of the No-Core Shell Model (NCSM) [10], which is a direct extension of the 
one-dimensional oscillator basis built in chapter HI This is a translationally invariant, 
anti-symmetric basis for each A, with a complete set of states up to a maximum exci- 
tation of iV max hw above the minimum energy configuration, where to is the harmonic 
oscillator parameter. Further details of the basis construction are explained in sec- 
tion Note that the realistic NCSM was built by Navratil [TU] and collaborators 
long before this thesis work. The work done here consisted of the code modifications 
necessary to implement the SRC The procedures used here build directly on those 
of chapter 0] where we studied the convergence benefits of SRG evolution along with 
a general analysis of the evolving many-body hierarchy. 

Just as in chapter H] with the one-dimensional model, we start by evolving H s 
in the A = 2 space, which completely fixes the two-body matrix elements. Next, 
by evolving H s in the A = 3 space we determine the combined two-plus-three-body 
matrix elements. We can isolate the three-body matrix elements by subtracting the 
evolved two-body elements embedded in the A = 3 basis. Having obtained the 
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Table 5.1: Definitions of the various calculations used to study SRG evolution. 



NN-only 



No initial NNN interaction and 

do not keep NNN-induced interaction. 



NN + NNN-induced 



No initial NNN interaction but 

keep the SRG-induced NNN interaction. 



NN + NNN 



Include an initial NNN interaction and 
keep the SRG-induced NNN interaction. 



separate NN and NNN matrix elements, we can apply them to any nucleus with their 
appropriate combinatoric factors exactly as described in the one-dimensional case. 
We are also free to include any initial three-nucleon force in the initial Hamiltonian 
without changing the procedure. We summarize in Table [5] the different calculations 
to be made for 3 H and 4 He to confront these questions. If applied to A > 4, four- 
body (and higher) forces will not be included and so the transformations will be only 
approximately unitary for those nuclei. The question to be addressed is whether the 
decreasing hierarchy of many-body forces is maintained as observed in chapter H] or 
the induced four-body contribution is unnaturally large. 

If needed in the future, we could evolve 4-body matrix elements in A = 4 and will 
do so when nuclear structure codes can accommodate them. The prospects for such 
work are straightforward. For A = 4 the Jacobi NCSM used here can be modified to 
produce the necessary matrix elements for evolution by the SRG. These matrices are 
very large (see appendix [D] on basis scaling issues) and may require more sophisti- 
cated techniques for computing the many simultaneous differential equations involved 
in the evolution. Implementing these algorithms should be straightforward, though 
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computationally intensive. For A > 4 we would have to implement the SRG in a dif- 
ferent type of NCSM which uses single particle coordinates and Slater determinants 
to antisymmetrize the states. This scheme is also well developed and, due to its basis 
organization, easy to scale in A but accordingly computationally intensive. 

The initial (A = oo) NN potential used here is the 500 MeV N 3 LO interaction from 
Rcf. [4Q]. The initial NNN potential is the N 2 LO interaction [67] in the local form of 
Ref. [69] with the additional LEC constants, cd and ce, fit to the average of triton 
and 3 He binding energies and to triton beta decay according to Ref. [70]. The values 
used here are cd = —.205 and ce = .029. We expect similar results from other initial 
interactions because the SRG drives them toward near universal form [2B]; a survey 
will be given in Ref. [7_T] . NCSM calculations with these initial interactions yield ener- 
gies of -8.473(4) MeV for 3 H and -28.50(2) MeV for 4 He compared with -8.482 MeV 
and —28.296 MeV from experiment, respectively. So there is a 20keV uncertainty in 
the calculation of 4 He from incomplete convergence (shown in the third decimal place 
- in parentheses) and a 200 keV discrepancy with experiment. The latter is consistent 
with the omission of three- and four-body chiral interactions at N 3 LO [72]. These 
provide a scale for assessing whether induced four-body contributions are important 
compared to other uncertainties. 

It is remarkable how well the one-dimensional model predicts behavior in these 
three-dimensional calculations. Often serious qualitative differences appear between 
one- and three-dimensional systems. In this case several factors conspire to damp 
those differences. First, the three-dimensional potential is usually stated in a partial 
wave basis, which breaks up the three- vector information based on angular momentum 
states. So, when we work in a given partial wave, we are working with an effective 
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one-dimensional potential, one which we were able to successfully mimic in chapter HI 
Secondly, and more central to this thesis, the SRG is a simple unitary transformation 
and its behavior is governed only by the qualitative structure of the matrices it is 
working to evolve, like the Hamiltonian and the generator. Qualitatively mimicking 
those, the outcome of SRG evolution will be predictive. This leaves us with a powerful 
analytic tool in the one-dimensional model of the NCSM to quickly explore other 
methods to improve many-body calculations beyond just SRG implementation. 

5.1 Induced Many-Body Forces 
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Figure 5.1: Ground-state energy of 3 H as a function of the SRG evolution parameter, 
A. See Table [5] for the nomenclature of the curves. 



In Fig. 15.11 the ground-state energy of the triton is plotted as a function of the flow 
parameter A. Evolution is from A = oo, which is the initial (or "bare") interaction, 
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toward A = 0. We use iV max = 36 and fko = 28MeV, for which all energies are con- 
verged to better than 10 keV. We first consider an NN interaction with no initial NNN 
("NN-only"). If the Hamiltonian is evolved only in an A = 2 system, higher-body 
induced pieces are lost. The resulting energy calculations will only be approximately 
unitary for A > 2 and the ground-state energy will vary with A (squares). Keeping 
the induced NNN yields a flat line (circles), which implies an exactly unitary trans- 
formation; the line is equally flat if an initial NNN is included (diamonds). Note that 
the net induced three-body is comparable to the initial NNN contribution and thus 
is of natural size. 
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Figure 5.2: Ground-state energy of 4 He as a function of A. Two hus (left and right) 
and two iV max 's (overlayed) are shown for comparison. 



In Fig. 15.21 we examine the SRG evolution in A for 4 He with Ml = 36 MeV. 
The two- and three-body matrix elements were evolved in A = 2 and A = 3 with 
-Wmax = 28 and then truncated to iV max = 18 (solid) and iV max = 16 (dashed) at each 
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A to diagonalize 4 He. The NN-only curve has a similar shape as for the triton and 
when the induced NNN is included, the evolution is close to unitary, the deviation 
being an indication of the induced four-body force. The pattern only depends slightly 
on strength and type of an initial NNN interaction. 

In both cases with and without an initial three-nucleon force, the dotted lines 
represent the converged values for their respective initial Hamiltonians. At large A, 
the discrepancy is due to a lack of convergence at iV max = 18, (the level of truncation at 
this calculation) but by A < 3 fm _1 SRG decoupling takes over and the discrepancy is 
due to induced four-body forces. This is illustrated nicely by the overlay which shows 
that, for A < 3 fm _1 , cutting from iV max = 18 to iV max = 16 introduces no more basis 
truncation errors, where above A = 3fm _1 there is a significant truncation effect just 
from 18 to 16. This indicates that the SRG has achieved convergence well inside an 
N max =16 basis for this calculation. 

Also, in Fig. 15.21 are shown the same overlay comparison for two hu) values, that 
which is optimal for the initial Hamiltonian to compute 4 He (fkv = 36) and one that is 
nearer the optimal for that potential evolved to A = 2 fm" 1 as expected from previous 
SRG studies. We can see the same converged results at low A in both fta>'s despite 
the different bare answer. This is consistent with being well converged and thus flat 
in hhj. 

The induced four-body forces then make up the difference to the converged value 
and contribute about 50keV net at A = 2fm~ . This is small compared to the rough 
estimate in Ref. [72] that the contribution from the long-rang ed part of the N 3 LO 
four-nucleon force to 4 He binding is of order a few hundred keV. If needed, we could 
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Figure 5.3: Binding energy of the alpha particle vs. the binding energy of the triton. 
The Tjon line from phenomenological NN potentials (dotted) is compared with the 
trajectory of SRG energies when only the NN interaction is kept (circles). When 
the initial and induced NNN interactions are included, the trajectory lies close to 
experiment for A > 1.7 fm -1 (see inset). 



evolve 4-body matrix elements in A = 4 and will do so when nuclear structure codes 
can accommodate them. 

The impact of evolving the full three-body force is neatly illustrated as a Tjon 
plot in Fig. 15.31 where the binding energy of 4 He is plotted against the binding energy 
of 3 H. The experimental values of these quantities are known to a small fraction of 
a keV and define only a point in this plane (at the center of the X, see inset). The 
well-known Tjon line (dotted) is the approximate locus of points for phenomenological 
potentials fit to NN data but not including NNN [73] . The SRG NN-only results trace 
out a trajectory in the plane that is analogous as it transforms NN forces unitarily 
but cannot account for the evolving 3NFs in the A = 2 space. In contrast, the 
short trajectory of the SRG when including the induced NNN interaction (shown 
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for A > 1.8 fm -1 ) highlights the small variations from the omitted four-nucleon force. 
Note that a trajectory plotted for NN+NNN-induced calculations would be a similarly 
small line at the N 3 LO NN-only point where the initial NNN is zero. 



5.2 Convergence in -/V max : G = T 2 
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Figure 5.4: Ground-state energy of 3 H as a function of the basis size iV max for an 
N 3 LO NN interaction [H] evolved with the SRG using G s = T re i to selected A's. The 
same calculation is shown at two values of fkv. 



Figure 15.41 shows the convergence properties of the three-body Hamiltonian as 
it evolves. The initial interaction used here is an N 3 LO NN potential with no 
initial NNN. The figure shows the same plot for two different hw's. On the left is 
hu = 28 MeV which is the optimal value for the initial interaction, and on the right 
is the same plot for t\w = 20 MeV which is more optimal for significantly evolved 
potentials (i.e., those evolved to A = 2.0 fm" 1 ). The calculations are performed by 
first evolving the initial Hamiltonian in the three-body space to the A indicated then 

119 



truncating the potential at each iV max -cut by setting matrix elements to zero above 
the extent of that iV max -cut. The curve for A = lOOfm -1 is essentially unevolved and 
is approximately equivalent to the initial interaction. 

Notice two major features as the initial Hamiltonian is evolved down from A = 
100 fm" 1 . First, the converged value at each A is the same - the evolution is unitary. 
As the SRG evolves the Hamiltonian in the three-body basis, it evolves the two-body 
matrix elements just as in the A=2 basis, but to keep the transformation unitary 
it must induce three-body forces to account for the high-momentum information it 
has transformed. These induced three body forces have been kept in the calculation 
of Fig. 15.41 and hence the ground-state energy is preserved by the transformation. 
Secondly, the convergence in iV max improves as A decreases down to some saturation 
point. The SRG is suppressing the coupling between states of large and small N max . 
Here we have made the choice of generator G s = T TC \, which was diagonal in the 
momentum representation, and should afford the same decoupling benefits it did 
there. However, in the oscillator basis, T re i has both a diagonal piece and a slightly 
off-diagonal piece that couples states at iV max and N max + 43 So in this basis the 
convergence saturates as the Hamiltonian evolves toward the shape of T re i and retains 
a certain amount of coupling. However, the evolution is itself basis independent and 
this apparent coupling is really due to the structure of the oscillator basis and the 
way we are required to cut it. Converting these Hamiltonians back to the momentum 
representation (or working there in the first place) should return to the situation 
where the potential asymptotically approaches the diagonal. Work along these lines 
is in progress [37] . 

18 see Eq. (|C.18[) for an example of this in the one-dimensional basis and a sample picture of T le \ 



can be found in Fig. IE. 11 

120 



The situation is already improved from the left panel to the right due to a more 
optimal value of fhv for the evolved Hamiltonians. The point of saturation is as far 
down as A = 1.0 fm _1 . One would like to optimize each point on this plot in fku. Since 
each combination of iV max and ftw represents a considerable amount of computation, 
a useful tool would be a simple orthogonal transformation of the Hamiltonian from 
one value of hu to another without building up a separate basis for each parameter 
set. Until such a tool is coded and tested, the alternative is to work at sufficiently 
large iV max that the energy is converged and relatively flat in (insensitive to) hu. This 
would result in a curve that may be less smooth and harder to extrapolate but still 
variational, converging at the large- iV max and optimal fkv for the bare potential. 
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Figure 5.5: Ground-state energy of 3 H as a function of the basis size iV max for an 
N 3 LO NN interaction [19] plus an initial NNN interaction [70j [71] evolved with the 
SRG using G s = T re \ to selected As. 



The initial Hamiltonians used in Fig. 15.41 are NN-only, so to increase the accuracy 
in calculating observables we must include initial three-body forces in our inputs. One 
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may be concerned that the exact form of any initial three-body forces will affect the 
form of the three-body forces induced to maintain unitarity. This question is tested 
in Fig. 15.51 which repeats the calculation of Fig. 15.41 but now including an initial NNN 
force at N 2 LO from Ref [67] as stated at the beginning of the chapter. We note that 
the shapes of the convergence curves in both figures are identical, indicating that 
the induced three-body forces are of natural size and are insensitive to the details of 
any initial three-body forces. Therefore the SRG is insensitive to the particulars of a 
given interaction, which is consistent with the observation that different Hamiltonians 
evolve to a near universal form. 
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Figure 5.6: Ground-state energy of 3 H as a function of the basis size N max for an N 3 LO 
NN interaction [49J with and without an initial NNN interaction [TOj [72] . Unevolved 
("bare") and Lee-Suzuki (L-S) results with Ml = 28MeV are compared with SRG at 
fku = 20 MeV evolved to A = 2.0 fin -1 . 

In Fig. 15. 6[ we again show the triton ground-state energy as a function of the 
oscillator basis size, iV max , to make a comparison between various calculations. The 
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lower, solid (upper, dashed) curves are with (without) an initial three-body force 
(see Table [5]). The convergence of the bare interaction is shown along with the SRG 
evolved to A = 2.0 fm _1 . The oscillator parameter hw in each case was chosen roughly 
to optimize the convergence of each Hamiltonian. Here, we compare to a Lee-Suzuki 
(L-S) effective interaction, which has been used in the NCSM to greatly improve 
convergence [751 176] . These effective interactions consist of a block-diagonal type 
unitary transformation to achieve an effective interaction in the oscillator basis at low 
-^Vmax • The resultant matrix is dependent on the model space of a given nucleus in 
contrast to the free-space nature of the SRG. Notice that the L-S effective interactions 
are not variational in iV max due to their model dependence and therefore it is difficult 
to extrapolate their behavior. 

The SRG calculations are variational and converge smoothly and rapidly from 
above with or without an initial three-body force. The dramatic improvement in 
convergence rate is seen even though the xEFT interaction is relatively soft. Thus, 
once evolved, a much smaller iV max basis is adequate for a desired accuracy and 
extrapolating in iV max is also feasible. 

Figure 15.71 illustrates for 4 He the same rapid convergence with iV max of an SRG- 
evolved interaction. Here we show results in two different fouj's (in parenthesis) for a 
Lee-Suzuki calculation, and Hamiltonians evolved to A = 2.0 fm -1 and A = 1.5 fm~ . 
We can see the rapid convergence both in the main plot and more clearly in the 
inset. It is again smooth and variational despite being an approximately unitary 
calculation for this sector. In this case the asymptotic value of the energy differs 
slightly because of the omitted induced four-body contribution. The SRG-evolved 
asymptotic values for different hu differ by less than lOkeV, so the gap between the 



123 




max 



Figure 5.7: Ground-state energy of 4 He as a function of the basis size iV max for an 
N 3 LO NN interaction [32] with an initial NNN interaction [23 [13]. Unevolved (bare) 
results are compared with Lee-Suzuki (L-S) and SRG evolved to A = 2.0 fm" 1 at 
hw = 28 and 36 MeV. 



converged bare/L-S results and the SRG results is dominated by the induced four- 

nucleon forces rather than incomplete convergence. Convergence is even faster for 

lower A values [71], ensuring a useful range for the analysis of few-body systems. The 

induced four body forces here are about 30-60 keV, which is less than the 200 keV 

discrepancy with experiment considered to be due to the omission of initial 4NFs from 

xEFT. However, because of the strong density dependence of four-nucleon forcea^f , 

it will be important to monitor the size of the induced four-body contributions for 

heavier nuclei and nuclear matter. It may be necessary to evolve unitarily in the 

four-body space to include induced four-body forces, and even initial xEFT terms in 

this sector. However, it is doubtful that higher body forces (A > 4) will be significant 

given the strong four particle clustering nature of larger nuclei. 

19 This follows from the power-counting of the %EFT lagrangian. The coupling constant must have 
one more power of A -3 , or density, due to the extra factor of N^N. 
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5.3 Convergence in iV max : G = H] 



ho 



As shown in Figs. 15.41 and 15.51 the choice of SRG generator, G s = T IC \, is not the 
most efficient choice in the harmonic oscillator basis. Due to the organization of the 
states and their complicated momentum dependence, T re i is not diagonal in this basis 
and therefore cannot be expected to completely diagonalize the Hamiltonian. We can 
make a different choice of generator, as already mentioned in chapter HI G s = H^o, 
the harmonic oscillator Hamiltonian. This choice is diagonal in the oscillator basis, 
having the eigenvalues E n = (N max +3/2)huj of each state along the diagonal. In this 
section we check, in the three-dimensional case, the behavior of this choice of G s . 
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Figure 5.8: Ground-state energy of 3 H as a function of the basis size iV max for an 
evolving N 3 LO NN interaction without an initial NNN interaction. Three values 
of oscillator parameter, hu = 28, 20, and 12 MeV are shown for comparison. 



Figure shows the analog of Fig. l5.4l but using H ho as the generator. We see, just 
as in the one-dimensional case, that the convergence improvement does not saturate 
at a specific point in the evolution but keeps improving as the SRG evolves it further 
towards the diagonal in the oscillator basis. On the left and center are fiio = 28 
MeV, optimal for the unevolved Hamiltonian, and hw = 20 MeV for comparison to 
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the T re i case. The optimal fko for i?h -evolved potentials has not been fully explored 
for specific A's, though given the amount of increased convergence it is likely to shift 
further. We can see in the center panel with hw = 20 the convergence is not improved 
by the same proportions as it was between the T re i plots, but on the right we show an 
even smaller value, hw = 12, than was used for T re i. Here we see more improvement in 
the convergence with a more optimal hw for evolved potentials. However, at hw = 12 
we start with a less optimal oscillator representation of the Hamiltonian and therefore 
converge to an under-bound value. Again, this is where a hw switching tool would 
come in useful as discussed in the last section. The hw dependence would of course 
have to be tested for each value of A as discussed in Appendix In Fig. I5.9l we show 
again that the specifics of the initial three-body force do not alter the qualitative 
pattern of evolution. 




4 8 12 16 20 24 28 32 36 4 8 12 16 20 24 28 32 36 4 8 12 16 20 24 28 32 36 

N cut N cut N cut 



Figure 5.9: Same as in Figure l5\8l but now including the initial NNN force [7DJ [71]. 
Note the qualitative form of convergence is not affected. 

Unfortunately, the same spurious states appear here in the approximately unitary 
calculations as did in the one-dimensional version of G s = H^o, leaving concern about 
the renormalization by unnatural operators. Specifically, H^ Q contains the long-range 
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operator r 2 which is not present in the Hamiltonian and therefore one might expect 
it to induce unnaturally large components in various parts of the Hamiltonian as it 
tries to cancel spurious r 2 physics. On the other hand the choice G s = T rc i does not 
improve convergence very well for long ranged operators like, in fact, r 2 while the H^o 
SRG does. The focus of present investigations [12] has been to consider a hybrid G s 
of the form 



that is essentially T re \ but incorporates a small admixture of the r 2 operator to renor- 
malize long distance operators without contaminating the Hamiltonian's flow. This 
work has progressed well in the one-dimensional model, and the prospects are promis- 
ing for such a hybrid G s to renormalize both long and short distance operators with- 
out inducing large spurious many-body forces. Calculations to verify this behavior 
in three-dimensions have been slowed by coding complications that will be resolved 
soon. 

5.4 Radius Calculations 

In addition to binding energies, we would like to understand how the SRG affects 
other observables, especially long-ranged ones like the root mean square radius, 



where |0o) is the ground-state wavefunction of the v4-body system, Xi are the single 
particle coordinates, and x cm is the center of mass coordinate. To obtain r rms in terms 
of Jacobi coordinates we can use the identity [66] 



G s = T rc i + ar 2 



(5.1) 






(5.3) 
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where the are single particle coordinates and r\ is the first Jacobi coordinate. The 
second equality results from our definitions of the Jacobi coordinates and the fact that 
we are working in a symmetrized basis. This replaces the sum by a multiplicative 
factor, the total number of pairs in the A-body system. 

This calculation is performed directly within the three-dimensional NCSM code. 
We can evolve a two-body Hamiltonian and then use the code to compute r for 3 H 
or 4 He, and we can also evolve the three-body and embed it to compute r for 4 He. 
As of this writing however, complications with the NCSM code have prevented us 
from producing the unitarily evolved version. We cannot yet evolve A-body forces 
and compute r for the A-body nucleus. So, the pictures in the figures for 3 H and 
4 He are missing the induced three- and four-body forces respectively. This is of high 
priority in the near future. 

In Fig. 15.101 we show the proton rms radius of the triton as a function of iV max 
for two values of fko, 20 and 28 MeV (top and bottom). We plot several values 
of A along with the bare and Lee-Suzuki effective interaction for comparison. Note 
the poor convergence in basis size. Unfortunately, evolution with the SRG does 
not significantly improve convergence for this long ranged quantity, though when 
comparing the two values of hu the convergence is improved from 28 to 20 MeV. 

In Fig. 15.111 we show the same radius calculation for 4 He but for three ftw's from 
left to right 20, 28, and 36 MeV, which are optimal values for the evolved and bare 
Hamiltonians respectively. Here both NN-only and NN+NNN-induced calculations 
are shown, but again the full unitary evolution calculation (including induced four- 
body forces) is not yet available. Again there is some improvement in convergence for 
lower values of hu. Those bases have a lower infrared cutoff, Air, and thus keep more 
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Figure 5.10: The 3 H radius as a function of basis size iV max for two choices of hui, 
20, and 28 MeV, corresponding to the optimal values for the bare and some evolved 
Hamiltonians. The straight dashed line indicates the best converged value at the 
optimal hu (28 MeV) for the LS effective potential at N max = 36. 

long-distance information. Decoupling due to the SRG is then more beneficial to the 
convergence of the radius calculation. Unfortunately, the ultraviolet cutoff, Auv; is 
also lower and limits the accuracy of the converged binding energy, in turn limiting 
the accuracy of the radius. The black and red solid curves at hu = 20 approach the 
converged value much slower than at 28 or 36. At the same time, the green and blue 
curves are better converged (to the wrong value) at hu = 20. Of course, a higher 
iV max would help the situation by pushing A IR down and Auv U P- 
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Figure 5.11: The 4 He radius as a function of basis size iV max for three choices of 
hu>, 20, 28, and 36 MeV (corresponding to the optimal values for the bare and some 
evolved Hamiltonians) . The solid lines are calculations that include induced NNN 
forces. The dashed lines are NN-only calculations. The straight dashed line indicates 
the best converged value at the optimal f\w (36 MeV) for the LS effective potential 



Also, notice in Fig. 15.111 the spread in converged value (at iV max = 18) as a 
function of A. Here we present the bare operator acting on the evolved four-body 
wavef unctions. Not only are the effects of induced three- and four-body forces on the 
wavefunction contributing to the spread, but the evolution of the operator, necessary 
to preserve unitarity in this calculation, is not included. At hv = 20, this seems 
to be a much bigger effect than that of the three- and four-body parts the evolved 
wavefunction, considering a lack of any clear hierarchy between the intervals from 
the NN-only to the NN+NNN-induced and the NN+NNN-induced to the converged 
value. At the higher hu's, the r 2 operator is evolved much less though the converged 
value is more accurate. Here the three- and four-body hierarchy in the wavefunction 
is more evident. 
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As discussed in the previous section, a possible strategy around this has focused 
on using a hybrid generator, G s = T rc \ + or 2 , where a is a constant determining the 
mixture of the two terms. With this choice, one hopes to balance the renormalization 
of short- and long-range parts of the potential and increase convergence for all op- 
erators [12]. Results of this work in one-dimension are promising with regard to the 
ability to tune the SRG to given operators and basis structure. As noted in several 
places, progress in the three-dimensional code has been slowed by coding issues. Most 
of these relate to the same technical issue of explicitly building the r 2 operator in the 
NCSM basis. Upon resolution of this problem, many more calculations will become 
available quickly. 

It must be stressed that difficulties in finding an optimal SRG generator stem not 
from the SRG itself, but from the particular basis we have been working in. The 
harmonic oscillator basis is convenient to work in, both because it is variational in 
two parameters with straightforward cutoff implications, and because the center of 
mass separates trivially from the calculations. However, it is not the best choice for 
describing certain physical properties, especially long-ranged operators. The oscillator 
functions simply do not have the same fall-off as an operator like r 2 . Other bases 
that may prove useful to this work include hyper spherical harmonics, and a basis of 
correlated gaussians as used in the Stochastic Variational Method (SVM) [M]. These 
bases should be more efficient at resolving long-distance observables, such as nuclear 
radii, and the choice of SRG generator may not be as constrained. However, work 
using these bases has not been as developed as that of the NCSM, especially work to 
implement the SRG, but progress is on the horizon. 
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CHAPTER 6 



CONCLUDING REMARKS 



6.1 Recapitulation 

For many decades progress in nuclear structure theory has been hampered by 
the complexity of the nuclear interaction. A strong repulsive core, significant tensor 
force, and poorly understood many-body forces have all contributed to difficulty in 
accurately calculating nuclear observables consistently These features are associated 
with high- momentum/short distance structure of the nuclear interaction, the details 
of which are irrelevant to the low-energy structures of bound nuclear systems. More 
precisely, infinitely many different choices of high-energy details can produce the 
same low-energy results. The difficulty in many-body calculations stems from the 
fact that high- and low-energy states are coupled. These coupling matrix elements 
of the potential require an increased basis size to represent the relevant low-energy 
physics. Such extra size increases dramatically with the number of particles in the 
system under consideration. This situation is ripe for a method, such as the SRG, 
that damps high-energy details without adverse effects to the physics relevant to 
low-momentum interactions. 
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In recent years many inroads have been made into the problem of expensive com- 
putation being spent on irrelevant interaction components. Foremost is the develop- 
ment of xEFT's, a model independent formulation of pion exchange interaction the- 
ory in the mold of Yukawa, but now informed by the underlying principles of QCD. 
Among its successes, these potentials are able to parametrize many-body forces con- 
sistently and robustly with controlled errors at each order. However, while they are 
much softer than their phenomenological predecessors, they are still computationally 
expensive for computing all but the lightest nuclei. Also of note is the V[ ow fc renor- 
malization technique, which has provided a proof of principle that one can unitarily 
transform a given potential to decouple the irrelevant high-momentum degrees of free- 
dom and retain the relevant physics inside a smaller space. However, this technique 
has not left a clear path to the consistent treatment of many-body forces inherent in 
the xEFT (or any other) potential on which it is used. 

The Similarity Renormalization Group has proved to be a powerful tool for ad- 
dressing the above general problems. Its flowing unitary transformations leave the 
input physics unaffected while transforming the potential to a softened form. Imple- 
mentation of the SRG is simple in any basis and on any given free-space Hamiltonian 
while allowing freedom of choice as to the general form of the flow operators. Fur- 
thermore, the SRG provides a self-consistent treatment of many-body forces easily 
understood in either first or second quantized forms. It has been demonstrated to be 
robust in its scaling of both evolution times and residual errors. With these features, 
the SRG is well suited to address the needs of the nuclear structure community. 
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The predominant improvement provided by SRG evolution of a given potential is 
the decoupling between low- and high-energy states in the two-nucleon (NN) inter- 
action as explored in chapter [31 The degree of decoupling achieved by evolution is 
studied by measuring errors induced by an arbitrary cutoff on evolved potentials. In 
this way we can see the clean perturbative scaling of the residual coupling above a 
given amount of evolution denoted by the momentum scale A or s = 1/A 4 . Decoupling 
in the two-body sector was found to propagate universally to many-body calculations 
and was checked explicitly for few-nucleon systems. In addition, an investigation into 
the behavior of other SRG forms, set by the choice of the generator G s , showed no 
significant impact on the general decoupling properties expected at a given A. 

The formal discussion in section I3T21 shows that, in principle, the same mechanism 
of decoupling also improves the convergence properties of the three-body (and higher) 
forces in an appropriate three-body (A-body) basis. In addition the SRG induces 
three-body forces to account for the changes in the two-body potential and keep the 
evolution fully unitary in the three-body system. In fact, it induces A-body forces in 
an A-particle space and we must verify that these induced forces can be controlled in 
order for the two-, three-, and four-body evolved potentials to be useful to many-body 
calculations. 

These issues were explored in chapter 2] through the use of a one-dimensional 
model in a symmetrized harmonic oscillator basis inspired by the highly successful 
No-Core Shell Model. We demonstrated that the induced three-body forces were as 
large as expected from previous calculations, and that induced four- and higher-body 
forces exhibited a hierarchy of decreasing importance. We covered several other topics 
relevant to realistic calculations such as the convergence properties of the three-body 
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evolution, the possibility of fitting omitted induced forces (thereby avoiding an explicit 
evolution in a higher-body space), and evolution of external operators involved in 
scattering processes. All of these calculations serve to build experience and intuition 
for the realistic three-dimensional case where angular momentum and other details 
complicate the calculations. 

In chapter [5] we achieved the first evolution of realistic three-body forces using 
the antisymmetrized Jacobi harmonic oscillator basis of the No-Core Shell Model. 
Here we duplicated the calculations tested in one-dimension, showing induced three- 
and four-body forces as a function of the evolution parameter, A, and studied the 
improved convergence in 3 H and 4 He energies due to the decoupling achieved in the 
three-body basis. The same hierarchy of induced forces was exhibited here. We also 
took a look at a long-range observable, the rms radius of these nuclei, to extend work 
being done by others in one dimension. 

Our experience from the ID calculations guided progress rapidly and that model 
continues to be a testing ground for new calculations. All the features of many-body 
evolution observed in the one-dimensional case were born out semi-quantitatively with 
the realistic interactions. For instance, all of the induced three- and four-body forces 
were in the same proportions to those in one-dimension. Often one is wary of trusting 
the results of a model calculation in one dimension to provide qualitative predictions 
of three-dimensional physics. Here, the basis structure is very similar to that of the 
three-dimensional problem which is already broken up in a partial-wave expansion. 
The SRG is only affected by the general form of the matrices, Hamiltonians and 
kinetic energies, involved in its flow. 
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6.2 Plans for Future Investigations 

A broad range of topics should be considered, or reconsidered, in light of the 
advantages presented by the Similarity Renormalization Group. Here we discuss a 
few of the most pressing issues. 

While this thesis has shown the induced many-body forces to follow a decreasing 
hierarchy in many-body calculations, they will continue to be a source of error in 
all sectors from light nuclei all the way up to nuclear matter. Interactions evolved 
in a three- or four-nucleon basis are unitarily equivalent for their respective nuclei. 
However, the use of those potentials, embedded to higher-body spaces, is complicated 
by the A-body forces that were not induced in the evolution at a lower space. For 
example, in chapter [5] results showing increased convergence in 4 He also show a small 
variability, with A, in the converged value. This cutoff dependence is small in this 
calculation, smaller than the error with respect to the experimental value, but that 
space is only one nucleon higher than that in which the SRG was applied. More 
calculations are needed to test how much the cutoff dependence will be reduced by 
including just the induced three- and four-body forces. 

A first priority to achieve calculations of larger nuclei is to make the A = 3 evolved 
potentials available in a usable form to those groups who can perform calculations 
of larger nuclei. The coupled cluster community uses the single-particle Slater deter- 
minant basis to perform their calculations. They have so far achieved much success 
pushing the boundaries of their calculations with NN-only evolved potentials, reach- 
ing as high as 56 Ni. Likewise, NCSM calculations which use the Slater determinant 
basis can also reach much higher nuclei and would be complementary calculations. 
Codes already exist to convert the matrix elements in the Jacobi basis used here 
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to the Slater determinant basis, or "m-scheme", so progress on this front should be 
rapid. 

While induced many-body forces in all sectors will be a source of error, this thesis 
has shown them to follow a decreasing hierarchy in many-body calculations. To avoid 
applying SRG evolution in increasingly large A-body bases, one might control the 
cutoff dependence by fitting the omitted induced forces to terms inspired by the xEFT 
expansion. An implementation of this idea was made in chapter |4] with encouraging 
results and further progress should come easily. Through such a fitting program 
the pattern of induced many-body forces may even lend itself to extrapolation and 
increased control on the errors present in approximately unitary calculations. 

Another avenue currently under development to test the viability of SRG evolved 
interactions in larger nuclear systems is to calculate nuclear matter properties using 
the Hartree-Fock approximation. Here the SRG is applied in the second quantized 
form demonstrated in section 12. 21 Instead of evolving the free-space interaction, the 
SRG is applied in the medium where normal ordering has different consequences with 
respect to the induced many-body forces. The induced many-body forces can then 
be controlled in a qualitatively different way than the free-space calculations. 

As mentioned previously the major goal of the UNEDF collaboration is the devel- 
opment of a density functional for nuclei which incorporates microscopic inputs such 
as those from xEFT's. In this ongoing project, the SRG is playing a vital role by 
softening initial interactions [77]. 

Fundamental to the idea of convergence within a particular basis size is the ques- 
tion of correspondences between various cutoffs in different bases. While the cutting 
done in the momentum representation in chapter [3] is a straightforward ultraviolet 
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cutoff, the same procedure in the oscillator basis requires a different interpretation as 
briefly discussed in appendix IB.31 Furthermore, many basis choices exist, each being 
more suitable to a different physics task. While the NCSM is a relatively straight- 
forward and familiar basis, other bases like hyperspherical harmonics and correlated 
gaussians are promising choices for describing other nuclear features. In fact, an im- 
plementation of the SRG in three-body momentum representation is currently under 
development, and should provide a check on the evolution done here in the NCSM [37] . 

One of the major advantages of xEFT's is their natural framework for the inclusion 
of external operators, for example electromagnetic interactions. The SRG can also 
be used for consistent, independent, evolution of these and other operators. Work 
on this topic is already proceeding [13] and involves, among other things, the one- 
dimensional tool developed for chapter H] for preliminary studies of operator evolution 
in a many-body space. 

While reproducing the static properties of nuclear bound states is an important 
task for testing our quantitative understanding of the nuclear interaction, we must 
address the dynamics associated with reactions. One current avenue is to couple 
the NCSM technique with the Resonating Group Method, which allows access to 
information on scattering amplitudes between states within the NCSM basis. These 
calculations will contribute greatly to the knowledge of how nuclei interact with one 
another and are a vital diagnostic tool for many experimental setups, especially the 
inertial confinement fusion experiments coming online at the National Ignition Facility 
(NIF) at Lawrence Livermore National Laboratory. This work is an exciting new 
direction for research using the SRG. 
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The SRG can be applied term by term, as demonstrated by the vacuum expec- 
tation value analysis in section 14. 4[ to explore the interplay between different parts 
of the flow equations. This analysis needs to be continued in the realistic NCSM 
where we have access to the different contributions to the three-body force in xEFT. 
Further down the road, one can envision an analytical application of the SRG to 
gain general insight to the SRG's treatment of xEFT potentials (or any other physics 
one wishes to address). The diagrammatic decomposition of SRG evolution would 
inform and provide checks on the results of such work, but understanding the effects 
of SRG evolution on individual coupling constants may have significant advantages 
to the formulation of xEFT's. Specifically, an analytical understanding of the SRG's 
behavior may provide a practical means to extrapolate its behavior with regard to 
induced many-body forces. 

The Similarity Renormalization Group has enjoyed a meteoric rise in its provision 
of benefits to nuclear few- and many-body problems. This thesis has provided a 
review of some of its basic features and initial applications to address calculations. It 
appears much more is to come. 
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APPENDIX A 



CHIRAL EFFECTIVE FIELD THEORIES 



The original efforts to build Chiral Effective Field Theories (xEFT's) with more 
than one nucleon were based on a suggestion by Weinberg [6] to use a non-perturbative 
treatment of the power counting rules from the perturbative Chiral Perturbation The- 
ory (xPT). A non-perturbative development is necessary in order to describe nuclear 
bound states. Weinberg used a power counting scheme to build the nuclear po- 
tential and used the Lippmann-Schwinger equation to iterate it non-perturbatively. 

The motivation for this non-perturbative scheme is the reproduction of the deuteron 
binding energy, scattering lengths, and other low energy observables. The above dis- 
cussion of renormalization scales and momentum contributions is important to this, 
but it does not solve the bound state problem completely. Naive dimensional anal- 
ysis (NDA) can only give an estimate of the effective coupling constants, referred to 
as Low Energy Constants (LECs). Only a careful calculation and renormalization 
procedure can fix their values and even signs. The power counting rule used here is 
based on an assumption of naturalness. That after renormalization, the LECs will 
be of order one. Also, the LECs are assumed to serve as counter terms to loops in 
the irreducible diagrams. In practice, in xEFT these absorptions are achieved by a 
fit to data. 
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Shortly after Weinberg's proposal there was an attempt to develop a perturbative 
treatment of pions in an EFT. This was formally worked out by Kaplan, Savage, and 
Wise [78], and is known as KSW power counting. While still non-perturbative in the 
contact interactions, the authors conjectured that the pion dynamics might be weak 
enough to treat perturbatively. Unfortunately, it turned out that the breakdown scale 
in this approach is too low to be of practical use in nuclear structure, though new 
variations have been proposed [79|. This approach will not be discussed here further. 

A.l Basics of power counting 

To explain the method of power counting due to Weinberg and used in Ref. |80j . 
we will start by writing out part of a chiraly symmetric Lagrangian for xEFT 



£(0) = Id-.^TT-lMV^tft 

2 2 



N 



^C s (N j N)(N j N) - ^C T (N j aN)(N j aN) + ... 



£(i) = Art 



4 Cl M 2 - ^M 2 tt 2 + |^ 2 + J|(cVr ■ Pn) 

C4 



2F2 



tijk tabc 0-iT a (Vj 7T 6 )(Vfc 7T C ) 



N 



° -(N j N)(N j arN) -Viz --E (N j N)(N j rN) ■ (N j rN) + ... (A.l) 



4iV " 2 

where the superscripts refer to the chiral dimension, Aj, of a vertex i. The operators 

u and r are Pauli spin and isospin matrices. The epsilons are Levi-Civita symbols 

necessary to achieve proper tensor multiplication. The constants and qa are 

associated with basic pion interactions and are determined elsewhere 20 !. The LEC's 

here are q, Cs,t, D, and E. Notice E is the LEC corresponding to the first 3NF contact 

20 F-x and <?a are parametrizations of weak nuclear interactions between pions and nucleons. They 
can be measured in processes such as pion decay and nuclear (3 decay. 
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term. This Lagrangian is complete up to N 2 LO including only isospin symmetric 
terms. 

In this power counting approach we calculate the T- matrix directly via the Lippmann- 
Schwinger equation (given here schematically), 

T = V + VG T , (A.2) 

where we will define the potential, V, to be the sum of all irreducible diagrams up to 
a certain order in v as defined below. An irreducible graph is defined as one which 
has no purely nucleonic intermediate states. The need for this definition will become 
apparent momentarily. The intermediate state Green's function, Go will play a crucial 
role in the power counting and will be defined also. 

All irreducible graphs can be categorized by their contribution ~ 0(Q/A) u gov- 
erned by the power counting rule 



v 



2(N + L-C- 1) + ^V l A i and A* = d; + m/2 - 2 . (A.3) 



The variables N, L, C, Vj, Aj, di, and rij are the numbers of nucleons, loops, connected 
pieces, vertices of type i, chiral dimension of vertex type i, derivatives or factors of 
at vertex i, and nucleon lines at vertex i respectively. At any given order only 
a finite number of time-ordered irreducible graphs are needed to compute the exact 
potential at that order. Chiral Symmetry sets a lower bound at Aj > 0. The simplest 
graph we can have is one with no nucleons and two pions (pion propagator ~ mfy or 
two nucleons and one derivative provided by a pion interaction or four nucleons and 
no pions (contact interaction). Those are the lowest-order options in the Lagrangian 
in Eq. f lA.lj) . Given the bounds on Aj, the lower bound on v is 0, thus we have a 



lowest order in the expansion and a valid power counting scheme. 

148 



Let us apply Eq. (1A.3j) to a few example graphs from Fig. II .51 in section [L~4l Take 
the simple 1PE graph at LO. The graph has two vertices, each with two nucleon lines 
and one derivative (from pion interaction) giving Aj = 0. There are 2 nucleons in 
this graph, no loops, and 1 connected piece. The grand total is v = (LO) which 
matches with the table. For a different example, consider the IPE-contact graph for 
a 3NF at v = 3 (N 2 LO). There are two vertices, one with two nucleon lines and one 
pion based derivative gives A« = 0, and the other with four nucleon lines and the 
pion interaction derivative (the term with the LEC "D") gives A,; = 1. So, with the 
sum over vertices giving one, three nucleons, no loops, and one connected piece we 
have a total of v = 3. So, all graphs are organized according to their contribution Q v 
given by Eq. (|A.3[) and Eq. flA. ID . Note that no contributions survive of order v — 1. 

The derivation of Eq. (|A.3[) appears in many forms throughout the literature [SB 
1821 155] . Start by considering a general graph and sum up momentum contributions 
via naive dimensional analysis. The momentum contribution, Q, for that graph is 4 
for each loop integration (L) , —2 for each internal pion (P), —1 for each internal 
nucleon (/), +1 for each vertex derivative(D). We obtain the expression 

v = 4L - I - 2P + v i d i ■ ( A - 4 ) 

i 

We can change the form of this using the well-established topological identities 

L = P + I-^Vi + l, (A.5) 

i 

and 

2I + E = Y,V tniJ (A.6) 

i 

where rij is the number of nucleon lines at vertex of type % and E is the number of 
external nucleon legs (E = 2N). The particular form of Eq. (1A.3I) presented here also 
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requires some non-trivial transformations associated with the formulation of graphs 
in an energy independent way [81] such as the method of unitary transformations that 
separates the purely nucleonic part of the Fock space from the rest [M]. Finally one 
arrives at the convenient expression of Eq. (1A.3|) for v. The variables N, L, and C are 



easy to count and Aj provides us with a way of organizing vertex types as discussed 
above. 

A. 2 Many-body forces 

One of the beautiful aspects of using field theoretic techniques to describe the 
nuclear interactions is the consistent and straightforward generalization to few and 
many body forces. In chapter [H we briefly mentioned why few-body forces are guar- 
anteed to exist in the xEFT expansion. We argued that the renormalization intrinsic 
to an EFT is responsible for integrating out the high-energy effects, so that in a 
few-body problem, certain components of the interaction will be hidden from the ex- 
periment (i.e. they will occur at distance scales much smaller than we can resolve). 
These effects can be divided roughly into three types. First, the A resonance can 
play a role in the interactions and is discussed in the next section. Second, simple 
high-momentum nucleonic intermediate states like a box diagram with a large mo- 
mentum running around the loop. Third, relativistic corrections which modify the 
kinetic energy of the intermediate nucleon states (i.e. through pair production) allow 
for irreducible diagrams with high-energy intermediate states that are (at low lab en- 
ergy) highly virtual and will occur over very small distances. As stated above, these 
are the major contributions that cannot be resolved at the relevant lab energies and 
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(a) (b) (c) (d) (e) 

Figure A.l: 3 and 4 nucleoli amplitudes that occur at v — 2. These terms will cancel 
each other out in an energy-independent formulation. Figure courtesy of Ref . [80] • 



therefore must be renormalized into a 3N contact force. Thus a 3NF appearing in the 
EFT Lagrangian is inevitable due to any of the above three effects. 

A simple NDA consideration leads us to expect that the basic 3NF terms should 
become relevant at v — 3. To see this, consider the 3NF contact term compared to 
the NN contact term at LO {y — 0). All other things being equal, the only difference 
should be that the 3NF has an extra pair of nucleon legs (NN^) each of order 3/2 for 
a total of 3 extra orders of the momentum scale Q. Thus we should expect 3NFs to 
become relevant at v = 3. 

If we use the power counting of Eq. f1A.3h to build the simplest three- and four- 
body diagrams we find the first contributions appearing at v — 2 which seems to 
contradict the above NDA argument. However, these first 3NF and 4NF terms (shown 
in Fig. lA.ip will not actually contribute. If we draw out all possible graphs from the 
possible topologies at this order we will find that they cancel each other out. This 
vanishing is evident when we formulate the nuclear force in an energy-independent 
way [851 [86]. 

Finally, at v = 3 we will find the leading 3NF terms. With E = 3, L = 0, C = 1, 
and one A = 1 vertex we can draw all such diagrams as shown in Fig. IA.21 The 
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(a) 



(b) 



(c) 



Figure A. 2: The leading order 3NFs at v = 3. Figure courtesy of Ref. [80 



contributions to the potential from these three graphs are given by [80 

2 



- e 

i+j+k 



1 ( 9 A 



2 \2Fj ;^ + m/^' 



with 



and 



and 



y lir, cont 



A Cl Ml 2c, 



p2 p2 



Qi ■ Qj 



X 



i+j+k 



8F2 q 3 ^ + M 2 



(A.7) 



(A.8) 



(A.9) 



KSt^E^;-^)- (A.10) 

where z, j, refer to nucleons, and $ = p^' —pf, Pi (pi ') is the initial (final) momentum 
of the nucleon i. These will involve vertices with Aj — 1. These expressions for 
the potential look complicated, but they must reproduce the structure of the pion 
interactions that they are parametrizing. Thus the intricate structures of spin and 
isospin matrices. 

At v = 4 we find the first corrections to the 3NF terms and the leading order 
4NF's. Some of these diagrams are shown in Fig. IA.3I while a more complete list 
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Figure A. 3: Sample 3NF corrections and leading 4NFs at v = 4. A more complete 
list of 4NFs is in Fig. IA.41 Figure courtesy of Ref. [80] 



of 4NF diagrams can be found in Fig. IA.4I from ref. [87], which also presents the 
analytic expressions for the potential contribution from these graphs. It has also 
been shown that all disconnected diagrams at this order vanish completely in an 
energy independent formulation [84] . 

In many-body systems we must also consider disconnected graphs that involve 
fewer-body interactions and spectators. This is accounted for in Eq. flA.31) by the 
value of C. For instance a contribution to a three-body system might include a 
diagram with a 2N interaction and one spectator nucleon. Another might be a four 
nucleon system with a 3NF interaction and one spectator or two disconnected NN 
interactions (such as the v — 2 interactions in Fig. IA.ll) . Such disconnected diagrams 
are expected to be two orders of momentum more important than their corresponding 
tree graphs, though many formally vanish and must be carefully worked out [SO] . 

It should be clear that the complexity and sheer number of graphs increases quickly 
with the order in momenta {v) under consideration. However, numerical calculations 
indicate [87] that graphs with v > 4 have a negligible contribution, though these 
results are very preliminary and not conclusive. 
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A. 3 The role of the A resonance 

The A-isobar is a baryon with mass = 1232 which is only 300 MeV above 
the mass of the nucleoli. The A is composed of 3 non-strange quarks and being in 
the quartet of isospin 3/2 and spin 3/2 can have an integral charge from —1 to +2. 
For historical reasons, one might also hear the A referred to as an excited state of 
the nucleon especially given the strong resonance in nN scattering. This view comes 
out of the simple quark model where the A is a very short lived (decays via strong 
interactions) particle whose main decay channel is to a pion and nucleon. 

So far we have ignored the question of the effects of an explicit A degree of 
freedom in the Lagrangian. Whatever effects there might be would be treated as short 
ranged and encoded in the LEC's of our existing theory. If they are not sufficiently 
short the LEC's will become unnaturally large. But some work has been done on a 
Weinberg-like power counting scheme with the assumption m& — ~ 277%^-. This is 
a simplification however since — is yet another small momentum parameter 
in which to make a possible expansion. 

An example Lagrangian for a xEFT involving an explicit A is [5U] 

£ = A f (id - Am) A+-^- (rfSTA + h. c.) • Vtt-LVAWA- (rfSTA + h. c.) , 

(All) 

where Ha and Dt are LECs and Si and T a are spin and isospin matrices. These 
matrices are 2x4 because they must describe the interactions between nucleons and 
deltas, the doublet and quartet channels respectively, in spin and isospin space. 

The first contributions involving explicit A's are NN graphs involving two pion 
exchange and shown in Fig. IA.5I One can see the A's dominance in the LECs at 
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Figure A. 5: Leading NN two pion exchange graphs with explicit A (double solid line) 
arising at v — 2. Figure courtesy of Ref. [5U] 



this order. The potentials corresponding to the one-A 2PE box graphs are identical 
in form to their analogous graphs without As. Matching up the LECs, we can see 
c 3 = 2c 4 = g 2 A /2{m& — m N ) and as above the mass scale appearing is (ttja — rn N ) 
rather than m w . 




Figure A. 6: Leading order 3NFs with explicit A arising at v — 3, an enhancement 
over the deltaless theory. Figure courtesy of Ref. [80] 



The 3NF terms with explicit A resonances show up at v = 2, instead of non-A 
3NF terms starting at v — 3. There is basically one major contribution at this order; 
the first graph in Fig. IA.6l is analogous to the leading 3NF graphs that had reducibility 
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cancellations. Now because of the existence of a A in the intermediate state, these 
cancellations don't occur. The other two 3NF graphs at this order formally cancel 
out due to Pauli exclusion [67] . 

In summary, the A can be included explicitly to achieve higher accuracy or it can 
be excluded for simplicity leaving those higher energy details to be encoded in the 
LECs of some A-less xEFT. 
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APPENDIX B 



GENERAL FEATURES OF THE HARMONIC 
OSCILLATOR BASIS 



B.l One-Dimensional Harmonic Oscillator Wavefunctions 



The one-dimensional harmonic oscillator functions in momentum space are given 
by [88] 



1/4 



iTmhuj2 2n n\ 2 
where H n are the Hermite polynomials 



exp 



V 



2mftw 



H„ 



V 



VmhuJ 



(b.i; 



H (y) = l 

Hi{y) = 2y 

H 2 (y) = Ay 2 -2 

H 3 (y) = 8y 3 -12y 



H 4 (y) = iV - A8y 2 + 12 



H 5 {y) = 32^-160^ + 1202/ 



(B.2) 



or more formally 



tf n ( y ) = (-l)«exp^exp 



(B.3) 
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We will often make the approximation 

H n {y) « (2y) B , (B.4) 

since this is the dominant term in the polynomial for large y. 

These oscillator functions are displayed in Fig. IB . 1 1 on the left, in a standard 
MATLAB color scheme, with reds positive, blues negative, and green zero. The 
vertical axis is increasing oscillator number, n, downward. The horizontal axis is 
momentum, p, from negative to positive k max . On the right, the truncated momentum 
delta function, 

max 

Yl Mp)Mp') « S{p - p')Q( Pmax - p) (B.5) 

n=l 

is a demonstration of the incompleteness of the finite oscillator basis. In the limit 
-Wmax — ► oo the sum should go to a delta function, S(p —p'). Here, because of the 
finite nature of the oscillator basis, the delta function is truncated at a point we will 
refer to as p max (referred to in other contexts as the ultraviolet cutoff, Auv)- The p max 
of the truncation and the extent of the last oscillator function in the basis correspond 
to one another. 

The coordinate space version of IB. II is 

mfiwr 2 



^ (r)= '^^) 6XP 



2h 2 



H„ 



mhw 
r 

h 2 



(B.6) 



Note that the only difference in form is that the factor of mhw is flipped in relation 
to the variable r as it was to p in IB. II This will be important later in an analytical 
discussion of infrared and ultraviolet cutoffs due to basis conversions. 
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Figure B.l: On the left are shown the harmonic oscillator functions. On the right is 
the truncated momentum delta function 5(p — p')0(p m ax ~ p)- 111 this MATLAB color 
scheme the red is positive and blue is negative. 

B.2 Truncation 

When converting a momentum space potential into the oscillator basis as in 

/km ax 

*l>n(p)ll>n>(p')VM(p,p')dpdp' (B.7) 
k max 

or 

/kmax 

(n\p) (p\V {2) \p) (p'\n')dpdp' (B.8) 

the potential becomes truncated due to the finite nature of the oscillator basis. When 
we try to convert back to the momentum basis we find that we have effectively 
multiplied by the truncated delta function of Eq. IB. 51 

This truncation is shown in Fig. IB.2l for different sized oscillator bases, each using 
the same value for the oscillator parameter, hw. Notice that as iV max increases we 
retain more of the original potential, and parts of it smooth out to their original form. 
There are also residual oscillatory effects left at the edges of the converted potential. 
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Figure B.2: A two-body potential in the momentum representation. The axes are 
initial and final momenta, k and k', from — k max to k max . On the upper-left is the 
original momentum basis two-body potential. The other three plots are the same 
potential after being converted to an oscillator basis and back. Notice the truncation 
of matrix elements above each p max of the oscillator basis. 



As stated above the truncation at p max corresponds to the extent in momentum 
space to the highest harmonic oscillator wavefunction in the oscillator basis, V'iVmax (p)- 
Two variables affect the value of p max : the oscillator parameter fad and the basis size 
A max • We will consider the dependence on each of these variables in order. 
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A quick look at the harmonic oscillator wavefunctions in Eq. IB. II shows a simple 
proportionality between p and ftw. The change inp max should scale like the square root 
of the change in ftw. Numerically this can be verified by looking at the furthest extent 
of oscillator functions for different ftu/s. Shown in Fig. IB.3l are the oscillator bases for 
-/Vmax =28 and /iu;=2,4,8 and 16. An increase in hu by a factor of four doubles p max . 
In Fig. IB. 31 looking down the vertical axis in the plots of shows a similar dependence 
on iV max . The curve formed by the edge of the oscillator wavefunctions looks like a 
square root function implying p max « n 1 ^ 2 . 




-40 -20 20 40 -40 -20 20 40 

Figure B.3: Harmonic oscillator functions up to iV max =28 for twj — 2, 4, 8, and 16. 
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We can measure these dependencies more quantitatively by picking out the mesh 
point at which each oscillator function passes a certain minimum threshold: 

|^n(Pmax)| = 4>min ■ (B.9) 

In Fig. IB.4I we show an example of such a measurement. The value VVn = -01 is 
used to plot the variation of p max with hu on the left and with n on the right. In this 
log- log plot the dependence on hw is a clean power law with a slope of .5 as expected 
from the oscillator wavefunctions dependence on p/y/hcJ. 




1 io 1 10 10 20 30 40 

hffl %x "W 

Figure B.4: Variation of p max with ftw on the left and with n on the right. 



B.3 Cutoff Derivation 

Now, taking a more analytical approach to understanding the cutoffs inherent in 
the finite oscillator basis, let us first consider the momentum space harmonic oscillator 
function of Eq. IB. 11 This function is essentially a gaussian shifted to the right by 
the largest iV max polynomial term. The local maximum will be correlated with the 
high-momentum cutoff, Auv> or p max . Making the substitution, U = p max / \/mfkj we 
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can write the large-n approximation by 

if}{U) = A(2U) n exp~ u2/2 . (B.10) 
Differentiating with respect to U, we get 

^ = (2U) n ~ 1 (n - (2U)U) exp~ u2/2 . (B.f f ) 
The maximum, at ^ = 0, yields 



U ~ y/n ► Auv = Pmax ~ \J W^max fa ■ (B.12) 

which is the large iV max behavior for the ultraviolet cutoff caused by converting to the 
harmonic oscillator basis. As expected, when either N max or hw is increased, then the 
scope of the basis also increases and we keep more of the momentum space potential. 

For the infrared cutoff, Air, we can make the same analysis, but now using the 
coordinate space oscillator wavefunctions where p 2 /mhw — ► mhwr 2 . We again ask 
about the A^ max dependence of the local maximum of ip{U). Here we get a maximum 
extent in coordinate space, 

^ (B.13) 



which implies an infrared cutoff 



Air-a/^t^ (B.14) 

V -'''max 

setting a bound on the size of objects that can be accurately described in a given 
A^max hw basis. When hw grows large, individual oscillations are large and lose reso- 
lution on the small details in the momentum basis potential corresponding to large 
r structures. However, high A^ max polynomials have many small oscillations at low 
momenta compensating for the large hw value. Thus, both Air and Auv are pushed 
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back (down and up respectively) by increasing N max which one expects as the basis 
is extended towards completeness. 

B.4 Three-Dimensional Harmonic Oscillators 

To build the NCSM in three-dimensions, we must use the isotropic Harmonic 
oscillator eigenf unctions [89] . 




where v is defined as 

with /i the reduced mass (m is now reserved for the magnetic moment quantum 
number) and uj is again the oscillator parameter. The normalization constant, Nja, is 




ipkim{r, 9, 4>) 



N kl r l e-» r V +1/2) (2z/r 2 )>UM) , 



(B.15) 




(B.17) 



The 



L k ^\2vr 2 ) 



(B.18) 



are generalized Laguerre polynomials are given by 




x a e x d n 
n\ dx 




(B.19) 



The first few generalized Laguerre polynomials are: 




1 




— x + a + 1 




(a + 3)x 2 (a + 2)(a + 3)x (a + l)(a + 2) (a + 3 
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The first several spherical harmonic : 

YnOA) = / ^ )( (^^| ^(cos^ , (B.21) 
are given by 

W0) = Y 2 ~\6A) = sin 2 ^ 

Y^\e,<j>) = iJi&nOe-W Y^\9,ct)) = §«/g an 9 cob Be'* 

^(M) = \\Jl™s6 Y 2 °(6,<P) = i^/f(3 cos 2 9- 1) (B.22) 

y-i (0, 0) = ^ ^ sin 0e** K, 1 (0, </>) = =£ sin cos Be* 

At large r, we can simplify the expression for ipkim as 

tuJrM) « y(^)' +3/ Ve-^ 2 (z/r 2 ) fc yr(^0) 

2k-\-£ 

- (sr ^iM) «-«-™- (B, 3 ) 



If we make the substitution 2/c + £ — > n this has exactly the form of Eq. (1B.6|) and we 



can make the same analysis of the cutoff as done in section IB~3l for the one- dimensional 
oscillator functions. Indeed, the quantity n = 2k + I is the quantum number for the 
energy of a state in the three-dimensional harmonic oscillator trap. The eigenvalues 
are 

E n = (2k + £ + 3/2)hw (B.24) 

with a degeneracy of (n + l)(n + 2)/2. 

So, the three-dimensional oscillator basis has the same variational properties as 
the one-dimensional discussed above. The exact values of hw and convergence in A^ max 
(maximum allowed n = 2k + £) will vary due to different potentials and dimensional 
factors such as the cube in the normalization factor. But the cutoff behavior is the 
same between the one- and three-dimensional cases. 
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B.5 Variational Properties 

Any harmonic oscillator basis such as the NCSM, both one- and three-dimensional, 
is a variational calculation in the two parameters, A max and frw. Larger A max allows 
inclusion of more of the momentum basis potential, increasing Auv- Larger fkv ex- 
tends the reach of a given polynomial of order A max but damps the resolution of the 
polynomials oscillations, raising both Auv and Mr- Smaller fkv increases the reso- 
lution but hampers the reach of the A max th polynomial, resulting in lowering both 
Auv and Air. So bigger A max is always better (except for increasing computational 
requirements) by raising Auv and lowering Air, but the optimal hw must result from 
a tuning of Air and Auv at a particular Af max . 




Figure B.5: Plots showing the variational nature of the harmonic oscillator basis and 
the effect of SRG transformations within it. 



These variational properties are shown nicely in Fig. IB. 51 The left plot shows the 
increasing convergence with the size of the basis, iV max . Here several select values 
of hit) show that one curve is optimal for most of the range in A max considered. As 
Amax grows larger, more polynomials are available for the expansion of momentum 
basis wavefunctions and the dependence on fko flattens out. The center panel also 
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shows this by plotting instead verses hui where the minimum is clear and the increase 
in iV max is apparent but less quantitative. Again the dependence on ftw flattens as 
JV^ grows. The overall goal in these calculations is to balance convergence in iV max 
(bringing a reduced sensitivity to hu) with the increasing computational requirements 
the larger basis demands. This is further explored in Appendix [Dl 

The far right plot of Fig. IB.5I shows the effects the SRG has on convergence 
properties of a potential in the oscillator basis. The curves are plotted by taking a 
potential at a large iV max (~ 40) and evolving it in that basis, then cutting it to a 
specified iV max (~ 28) before computing the binding energy 21 !. This procedure was 
repeated at multiple values of ftw for comparisons and we can see several effects. 
First, the SRG evolves the initial potential to a new and different Hamiltonian and 
therefore the optimal fkv for that Hamiltonian's expansion will be different. In fact, 
it will be lower due to the way in which the Hamiltonian has been altered; high 
momentum matrix elements have been reduced and simplified so that the Hamiltonian 
needs less basis size to achieve the same convergence as before. Therefore the basis 
can accommodate information with a smaller hu and gets the benefit of increased 
resolution due to the smaller fkv. In other words, a lower Auv is sufficient so a lower 
Air is affordable. Thus the optimal hw shifts downward. Furthermore, each step 
in the evolution increases the convergence at a given basis size and so the overall 
dependence on tvjj decreases and the curves flatten out as A decreases. 



21 This calculation is analogous to the situation in Ref. [35] where the NN potential was evolved 
in momentum space (N max — ► oo) and then converted to some finite oscillator basis. 
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APPENDIX C 



SYMMETRIC JACOBI OSCILLATOR BASIS IN ONE 
AND THREE DIMENSIONS 



C.l Jacobi Coordinates 

The initial (i.e., unevolved) one- dimensional Hamiltonian for A bosons of equal 
mass m with a local two-body potential has the first-quantized form (in units with 



where the are single-particle coordinates and the ki are single-particle momenta. To 
connect to the nuclear problem of interest that uses potentials given in a momentum 
basis (e.g., chiral effective field theory potentials), we calculate matrix elements using 
harmonic oscillator basis states in Jacobi momentum coordinates. This representation 
also provides a clean visual interpretation of the SRG evolution of potentials. 

With equal-mass particles, a convenient set of relative momentum Jacobi coordi- 
nates is defined by (for j = 1 to A — 1) 



where the ki are the single-particle momenta of the A particles. This is one particular 
choice of normalization for the Jacobi coordinates. This choice is especially useful here 



h=\) 






(C.2) 
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because it sets the reduced masses for each coordinate equal to each other, simplifying 
rotation operations performed later. We define the Fourier transform V(x\ — x^) to 
V(pi,pi) as 

V{pi,p' x ) = J V^^e-^-^dh , (C.3) 

where l\ = (x± — X2)/V2 is the coordinate conjugate to the Jacobi momentum p\. 
We introduce a set of harmonic oscillator states \rij) corresponding to each of the 
coordinates of Eq. (]C.2j) . so a general product basis state has the form 



A-l 

n (c.4) 

with rij = 0, 1, 2, ... , N max for each j. In the next section we discuss how to build 
linear combinations of these states that have the appropriate symmetry. 

C.2 Transformation Brackets 

Transformation brackets are the expansion coefficients in the oscillator basis of 
one system of coordinates in terms of another [901 El] - We apply them to relate 
two different choices of Jacobi coordinates. Here, we show the relevant transforma- 
tion using the three-particle harmonic oscillator states defined in Eq. ( 1C.4I) and then 
generalize at the end. 

The single particle momenta are k\, k-i, and k^. The unprimed Jacobi momenta 



[see Eq.( IC.2j) ] are 



1 

3 V 



P2 = J^((h + k 2 )/2-k 3 ) , (C.5) 
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and the primed coordinates are obtained from exchanging fc 2 and 



Pi 




(C.6) 



After some algebra, the transformation that exchanges the last two particles (i.e., ki 
and kj) can be written as 



which enables us to express the primed oscillator states in terms of the unprimed 
ones. 

We denote the three-particle oscillator basis by \n\n2) = 77^773 j 0) where we have 
set Tiw = 1 for simplicity in this appendix. Note that the transformation brackets 
will not depend on the value of huj. The transformation that exchanges the last two 
single-particle coordinates can again be written as 




(C.7) 
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The derivation of the harmonic oscillator transformation bracket follows directly as 
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The second line is obtained from operating the transformation on the creation opera- 
tors 7]\. The third line is the application of the binomial theorem. The fourth balances 
the oscillator creation and annihilation, and the fifth is just some simplification. 

In the general A-particle system the transformation to exchange the last two 
particles, k^-i and kA, can be written as 



VA-2 

Va-i 



i 

d+l 

d 
d+l 




(C.IO) 



where d — (A — l) 2 — 1 is the number of generators of the rotation group, U(A — 1), or 
the group U(A) with the center of mass coordinate held fixed. An expression for the 
bracket (nA-2 n A-i\ n A-2 n A-i) a(A-2), which appears in Eq. (IC.14I) . is obtained from 
Eq. (IC.9P by substituting the general coordinate transformation Eq. fIC.lOp for the 
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three-particle transformation Eq. (10.81) . or y/\/{d + 1) and \fd]Jl + d) for 1/2 and 
V3/2. 

C.3 Symmetrization 

We carry out the SRG evolution for each A-particle subsystem in a complete basis 
of properly symmetrized states, which will be linear combinations of the basis states 
of Eq. (IC.4I) . The symmetrization procedure is adapted from the procedure developed 
for NCSM calculations [TT1 IT2~| IT3] . This entails symmetrizing first the two-particle 
system and then using a recursive procedure to go from the (A — l)-particle basis 
to an A-particle basis. At each stage we keep only symmetric states, identified as 
eigenstates of the symmetrizer with eigenvalue unity. 

The two-particle system is specified by the oscillator number ri\. The symmetrizer 
is (1 + P 12 )/2, where is the exchange operator between particles i and j. Because 
-P^l^i) = ( — ~L) ni \ n i), the symmetrizer in the two-particle case has eigenvalue one 
acting on states with ri\ even and zero when acting on states with n\ odd. Thus the 
symmetric basis states have n\ even and we simply omit the odd states. Following 
conventions from Ref. [IT], we label these eigenstates as |iV 2 2 2 ), where N 2 is the total 
oscillator number of the symmetric state and i 2 is an arbitrary label that distinguishes 
states degenerate in N 2 - In the two-particle case the notation is trivial, with N 2 = ni 
even and i 2 = 1. We write eigenstate projection coefficients as (iV 2 i 2 ||ni) = 5N 2 ,n 1 (^ + 
(— l) ni )/2. These are referred to in the literature as the "coefficients of fractional 
parentage" . 
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A three-particle basis is specified by product states of the two-body symmetric 
eigenstates, \N 2 i 2 ), an d single-particle states with the oscillator number correspond- 
ing to the third particle, \n 2 ). The symmetrizer for this system is governed by the 
permutation group, S3, which can be defined by just two of its generators. Here we 
choose the permutation operators P12 and P23. The symmetrization operator can be 
written as 

S = hi + P12 + P23 + P2P23 + P23P12 + P2P23P2) ■ (C.ll) 
o 

We build this symmetrizer in the basis \N 2 i 2 ; ^2) = \N2i2) ^2) where the states \N2i2) 
are already eigenstates of P i2 with eigenvalue one, so Eq. (IC.lip reduces to S = 
(l + 2P 23 )/3. 

In this basis, the matrix elements of P23 can be expressed as 

(N 2 i 2 ; n' 2 \P23\N 2 i2\ n 2 ) = Sn'M^^W™*)* > ( c - 12 ) 

where = N2+U2 = N 2 +n 2 = N' and (n' 1 n' 2 \nin2)3 is the one-dimensional harmonic 
oscillator transformation bracket for particles with mass ratio 3 [92]. We construct 
these transformation brackets and generalize to mass ratio d in section IC.2I By diag- 
onalizing this symmetrizer we identify the symmetric eigenstates of the system as the 
ones with eigenvalue unity. We keep only those states and discard the others. This 
set of eigenvectors gives us the coefficients of fractional parentage, {N 2 i 2l n^H-A^Za), 
of the three-boson symmetric eigenstates, l-A^), in terms of the original partially 
symmetrized three-particle space, (iV^; ^2)- Note that i% is not trivial like i 2 , be- 
cause in the three-body system there are eigenstates degenerate in the total oscillator 
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number, N 2 + n2- The label 23 keeps track of those degeneracies. We find in the one- 
dimensional system of bosons that the fraction of symmetric basis states for A = 3 is 
about one-fifth. For A = 4 the reduction in number of states is above 90%. 

The symmetrizer can now be coded directly in the three-particle oscillator space. 
A picture of the resulting matrix is shown in Fig. IC.ll This matrix has the three- 
particle oscillator basis states organized in a block diagonal form because the total 
oscillator number, N = n\ + n 2 , is a conserved quantum number proportional to the 
total energy of a state. 

Remembering that we need only the physical eigenstates of the symmetrizer, we 
diagonalize the matrix and obtain the eigenvectors. We keep only those vectors 
which have nonzero eigenvalue (the physical states), thus resulting in a symmetrizing 
matrix about one-fifth the size of the original oscillator basis. A picture of the sorted 
eigenvectors is shown in Fig. IC.21 These are all the symmetric states as expressed 
in the three-particle oscillator basis. The symmetrizer in this basis is, of course, 
diagonal, and all the eigenvalues are unity. 

The color scheme in Figs. IC.ll and IC.2I is displayed to the right of each graph; 
green is zero, red is positive, and blue is negative. Notice that the eigenvectors are 
normalized to one. The easiest non-trivial case is the second vector where each of the 
two components is l/y/2 ~ .7, which is shown in a bright red. 

To construct the basis states for higher A, we generalize this procedure. To go 
from A — 1 to A we need only to symmetrize between the last two particles, so we 
construct the symmetrizer 




(C.13) 
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10 20 30 40 50 60 



Figure C.l: The full symmetrizer for the \N2i2] ^2) space, with both physical and 
spurious states. A small, iV max = 10, basis is shown for clarity. Note also that 
Matlab has chopped off the last column when making the figure; this is a symptom 
of my very simple plotting routine. 
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2 4 6 8 10 12 14 

Figure C.2: The physical eigenvectors of the symmetrizer. Note vector components 
corresponding to iV max blocks in Fig. IC.ll The same size basis as in Fig. \C.1\ iV max = 
10, is shown. Again, Matlab has chopped off the last column. 
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in the space of (A — l)-particle symmetric eigenstates and the additional Jacobi state, 
n-A-i- We label the basis of this space as {Na-x^a-x] ua-i)- The matrix element of 
the exchange operator in this space is 

{Na^'a-i, "a-iI-^a-i^I-^Va-iU-i; n A -i) 

= ^_ 1 +n' A _ 1! AT 4 _ 1 +n A _ 1 y^(-^A-A-l 11-^4-2^ A-2i^A-2) 

X (A^ A „ 2 2 A „ 2 ;ra A „ 2 ||^A-l«A-l)(«A-2 n A-ll n A-2«A-l)A(A-2)(C ; -14) 

where the sum is over Na-2, 2a-2> ^a-2 and n' A _ 2 . The only significant difference 
from the three-particle case is that we must sum over the components of the A — 1 
subcluster symmetric states to get all the contributions to the exchange of the last 
two bosons, riA-2 and ua-x- The parameter d = A(A — 2) can be derived by taking 
the last two Jacobi coordinates pa-2 and pa-i, as defined in Eq. (1C.2I) . and finding 
the transformation that exchanges particles labeled by &a-i and &a- This procedure 
was explained in section ICT21 

For fermions, we need a complete basis of fully anti-symmetrized states. If we 
consider the one-flavor case, the procedure for our one-dimensional model is a trivial 
modification of Eqs. (1C.llj) and (1C.13I) . namely all odd permutations come with a 



minus sign. Thus, for A = 3 the anti-symmetrizer can be written 

A = ±(1-2P 23 ), (C.15) 

where P 2 3 acts on the flavor space as well. If there are more flavors than particles 
and the interaction is flavor independent, the spatial wavefunction for the ground 
state will be symmetric and correspond to our boson ground state wavefunctions. 
For realistic three-dimensional nuclei, the required construction of an anti-symmetric 
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Jacobi basis with full angular momentum coupling has been worked out for the NCSM 
by Navratil et al. [HJ [121 H3] and is briefly discussed in section IC.51 

C.4 Hamiltonian Matrix Elements 

To obtain the Hamiltonian in the symmetric eigenbasis of the general A-particle 
system, we employ a recursive embedding procedure that utilizes the partially sym- 
metric bases developed for the symmetrization operator. First we treat the kinetic 
energy and then the potential. 

The relative kinetic energy in the three-particle system is the total minus the 
center-of-mass kinetic energies: 

-^rel = -^tot ^cm 

2 u2 l2 ru i U i U ^2 



k\ k\ k\ {k x + k 2 + k 



2m 2m 2m 2(3m) 



pl+p 2 2 



(C.16) 



2m 

where the p^s are defined in Eq. (]C.2|) . Momentum basis states are organized by 
increasing kinetic energy. We can project T re i directly onto the three-particle oscillator 
basis by using the ladder operator definitions of the Jacobi momenta. The projection 
of T re i into the \n\n2) basis is 

III \rri I \ / / I \ V\ J T V\\ \ 

2m 

= ^^K^l^l - Vif + (vt - ^)Vm 2 ) , (C.17) 
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where the rji and 772 operators act on the n\ and spaces, respectively. Continuing, 
we get 

{n x ri 2 \% cX \nxn 2 ) = ^ [(n[ | (r)\ - r/i) 2 |m}^2^ + (n^Kvt ~ ^fM^yJ 



—to 



(V (ni + l)(rai + 2) 5 n ' lini+2 + v^ifai - l)*r»i,m- 
- (2ni + l)5 n ' i>ni ) <5 n ^„ 2 



+ a/ n 2 (n 2 — 1) 5 n ^ n2 _ 2 
- (2n 2 + l)5 n>2 ) 5 re>1 ] . (C.18) 

As noted, we keep only the rii-even states using the projector (iV^ll^i), and we 
can symmetrize the full three-particle system with the symmetric eigenstates, \N3i3) 
whose components are given by (N 3 i 3 \\N 2 i2',n2). 

To derive the A-body kinetic energy in the symmetrized basis, (T A ) sym , we use a 
recursive procedure on the (A — l)-body result to find the A-particle space operator 
matrix elements: 

A-x 



(T A ) sym = (N' A i' A \T A \N A i A ) = {N' A i' A \Y^vVMN A iA) 

i=i 

= (N' A \\N' A ^n' A ^) [(N'^n'^KT^)^ 
+ p A _ 1 /2m\N A _ 1 n A -i)](N A _ 1 n A _ 1 \\N A ) 

= (N'jN'^n'^) [(N^T^^N^S^^ 
+ S N ' A _ ljNA _ 1 {n' A _ 1 \p A _j2m\n A - 1 )]{N A - 1 n A - 1 \\N A ) 

= (N'jN'^n'^) [ (T A _ 1 ) gym i )nA _ j 



1 ^_ 1 ,v A _ 1 (^_ 1 ,n A _ 1 (2n^_i + 1) - K< A _^n A _ x \jn\_ x - n A _ x 



— ( 

4 



^_ 1 -2,n A _V(«A-i + l)(n A -i + 2))] (N A ^n A ^\\N A ) , (C.19) 
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where we have suppressed the i^'s and i^-i's for simplicity after the first line. Inter- 
mediate summations over Na~i, n^-i, Ia-i, ^a-ij n> A-n an d are implicit. 

In the same manner as the kinetic energy, we can recursively embed the potential 
in the A-particle space, starting with the two-body interaction between the first two 
particles. Because we are working in fully symmetrized few-body spaces we do not 
need to consider all pair-wise interactions, but only one such pair and scale by the 
number of interactions. For instance, in the three-particle system the full two-body 
interaction is = V12 + V23+V13 = 3Vi2- In a general A-particle space, this becomes 
= ( 2 )Vi2- The matrix element of a two-body potential, Vu, in the relative coor- 
dinate harmonic oscillator basis, \rii), is (n^V^l^i) = J(n' 1 \p[)(p'\Vi2\p)(p\ni)dpdp' 
where the matrix elements of (p'\Vi2\p) are given by Eq. flC.3p . 

Once in the oscillator basis, embedding in a larger particle space is a straightfor- 
ward process. Starting with the two-body interaction, V12, the two-body oscillator 
symmetric states are isolated using the projector (iV^l^i) which picks out just the n\- 
even states. Embedding this interaction in the three-particle space involves adding 
a new Jacobi coordinate, |ri 2 ), to the existing system. With respect to the two- 
body interaction, V12, this additional coordinate is associated with a delta function, 
5 n2tTl > . Finally we obtain the symmetric three-particle states by using the projector, 
(N 3 i 3 \ \N2i2', n-i). Multiplying by („) = 3 gives us the full strength of the two-body 
interaction. 
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In general we can write this procedure as an expansion of the final A-particle 
symmetric space matrix elements of V\2- 



(V\% ym = (N'J A \V^\N A i A )=()(N' A i' A \Vi2\N A i A ) 

sym&n' A _-,,riA—i 



where again we have dropped the i A 's after the first line for simplicity and intermediate 
sums are implicit. We remind the reader that n A ~i can only take values from to 
N — Na-i, where N is the total oscillator quantum number used to organize the 
states. We start with the two-particle space and work our way up to the A-body 
space, embedding the interactions successively in each sector using Eq. (1C20j) . 

When symmetrizing Va we must embed the symmetrized V A -\ with the appropri- 
ate combinatoric factor included as explained above. This factor derives from the fact 
that we had embedded a 2-body force in the A — 1 space that is now to be extracted 
and embedded in the A-particle space. Thus we must remove the old factor ( A ~ x ) and 
multiply by the new (^) factor, which has the net effect of multiplying by A/ (A — 2). 

Any initial three-body force (discussed below) is embedded in the same manner as 
above except that it originates in the three-particle space. The initial three-body force 
is a function of two Jacobi momenta, which we transform directly into the partially 
symmetrized three-particle oscillator space and then use all of the same embedding 
procedures developed above. Note that two- and three-body forces must be embedded 
in higher spaces with different symmetry factors, (^) and (^) respectively. 
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In previous formulations of this recursive approach [IT] , subsequent potential em- 
beddings are achieved by making a change of coordinates for the last two Jacobi 
momenta. For systems with A > 5, the three-body force requires a similar change of 
coordinates for the last three Jacobi momenta. Such a scheme is unnecessary here. 

C.5 Three-Dimensional NCSM 

In this thesis, only the one- dimensional NCSM was built for SRG studies. How- 
ever, the three-dimensional codes were modified and used for the realistic calculations 
in chapter [5j In three dimensions the Jacobi harmonic oscillator basis differs in sev- 
eral significant ways that complicate the expressions involved considerably, though in 
a straightforward way. However, the basic form of the basis remains the same as in 
the one-dimensional case. 

The three-dimensional NCSM basis is, like the one-dimensional case, organized 
in iV max blocks, so scaling up the basis size or making cuts to study decoupling is 
a simple matter. It is variational in iV max and fko as noted in App. [B] It is built 
recursively from one A-body sector to the next, and the Hamiltonian is embedded 
in subsequent sectors with combinatoric factors, taking advantage of the symmetry 
properties. 

The three-dimensional NCSM was built to calculate systems of nucleons, which 
are fermions, while the one-dimensional model was built based on a system of bosons 
with no spin or angular momentum. So in three- dimensions, we have to build an 
antisymmetrizer to find fully antisymmetric states in which fermions would reside. 
This still employs the same strategy of building on the A — 1 subcluster. So a 
new Jacobi coordinate is added to the A — 1 space and the matrix elements of the 
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ant isy mmet r izer , 



X=^(1-(A-1)P (A _ 1M ), 



(C.21) 



are computed. As before, P(a-i)a is the permutation operator between the Ath and 
(A — l)st nucleons. Here the only difference is the minus sign. And the fact that P 
now also transforms angular momenta. 

Now in three dimensions, we use the three-dimensional isotropic harmonic oscil- 
lator wavefunctions as displayed in Sec. IB.41 The quantum numbers of the nuclear 
system are now n,£,s,j, and t, the radial, orbital angular momentum, spin, total 
angular momentum, and isospin of the state. In the three-dimensional oscillator, the 
conserved energy is N = 2n + £ which will impact embedding procedures from one 
iV^ sized basis to another. Now the anti-symmetrized states are labeled as 



where the (Na-iIa-iJa-iTa-i] ^a-i^a-iJa-iII^a^a^a^a) ar e the three-dimensional 
coefficients of fractional parentage. These are given by the eigenstates of the antisym- 
metrizer in the non-antisymmetrized basis just as their one-dimensional counterparts 
were obtained from the eigenstates of the symmetrizer. 



N A IaJaTa) = y^(iVA-lU-l^A-l?A-i; TlA-liA-ljA-l\ \ N A lA JaTa) 



x \Na-iIa-iJa-\Ta-\] n A -\^A-\3A-\) , 
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To build the antisymmetrizer, Xi the matrix elements of P(a-i)a i n the non- 
antisymmetrized basis can be written as 



where the quantities in curly brackets are the Wigner 6j and 9j symbols used to 
compute the couplings between the angular momenta and isospin quantum numbers. 
The last line is the three-dimensional oscillator transformation bracket as derived in 
various sources with various methods [901 EH E2] • 

Given the scaling properties of this basis discussed in appendix O, several tech- 
niques have been developed to improve the efficiency of the three-dimensional tech- 
nique. In Ref. [93] a method is introduced that allows computation of only some 
antisymmetrizer matrix elements instead of requiring the full antisymmetrizer. An- 
other technique deals with the cumbersome increase of angular momentum channels 
with A. Some authors [M] have developed techniques to better organize these chan- 
nels and make the coding more manageable. Another approach to this problem is 
to convert to the single-particle, or Slater determinant, basis where the all the an- 
gular momentum embedding procedures are more straightforward. However, those 



Sn>,N ^2,{N A - 2 iA-2JA-2TA-2] ^-2^4-2^-211^-1*' 



Ja~i^a-i) 



x{N A _ 2 iA-2JA-2TA-2]n A _2^A-23A-2\\NA-liA-lJA-lTA-l) 
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calculations require large clusters and sophisticated memory handling techniques to 
calculate even the lightest nuclei at small iV max . The details of the single-particle 
basis scaling are discussed in appendix [D] A one-dimensional single-particle basis 
was constructed and is presented in appendix [F] 



186 



APPENDIX D 



SCALING PROPERTIES 



In this appendix we document the computational scaling features of the various 
calculations in this thesis. An introduction to the scaling issues involved in schemes 
such as the NCSM is necessary to appreciate the computational improvements brought 
about by the SRG. 

D.l Momentum Space NN calculations 

The first problems addressed in chapter [3] involved applying the SRG to interac- 
tions projected onto a momentum representation. The mesh used here is most often 
gaussian quadrature with the grid dimension at about 100. Figure HXT1 shows the time 
as a function of the evolution parameter, s to evolve two different potentials, on the 
left a xEFT and on the right AV18, in several partial waves. We expect the evolution 
to be linear in s and indeed find such behavior. On the left plot, as s increases, 
we find no hint of stiffening of the differential equations. The dependence is strictly 
linear all the way to s — 1 (A = 1). 

The right, log- log, plot of Argonne partial waves shows something a little more 
interesting. At smaller s the evolution is also linear, but as s increases, the dependence 
in some partial waves stiffens somewhat. Some investigations suggest this may happen 
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Evolution Parameter, s Evolution Parameter, s 



Figure D.l: SRG timings vs the parameter s show no stiffening for large s (small A). 
The left shows timing results for N 3 LO (500 MeV) potentials in several partial waves 
and the right shows the same for Argonne vi S but in a log- log plot. 

when the potential has been driven to the desired decoupling and due to the larger 
number of matrix elements the flow is working hard to maintain a degree of accuracy 
that is unnecessary. It seems that freezing the high energy parts of the potential, 
when they have reached a sufficient evolution, works extremely well to improve the 
time scaling. 

D.2 One Dimensional Jacobi Model 

In the one-dimensional version of the no-core shell model, we must build a basis 
of harmonic oscillator states (see Chap. [U]). This basis provides for a variational cal- 
culation with the accuracy directly connected to the the basis size, -/V max . Therefore, 
we need to understand how the basis scales and the resources that are consumed 
to build larger bases. Figure ID. 21 has on the left the size of matrices for different 
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Figure D.2: Plots showing the scaling of One-D NCSM basis size and time to build 
those bases as a function of iV max for A = 2, 3, 4, 5, and 6 bosons. 

A-body spaces as a function of N XDax and in the center a log-log plot showing the 
time to build the list of states for those bases. Note the log-log inset on the left plot 
shows the same dependence in N max as the timing plot indicating the build time is 
linear with the number of matrix elements. Here we can see dramatic increase in size 
and time as A increases from 2 to 6. The far right plot recasts this data as basis 
dimension versus A for several choices of iV max . Again we can see the importance of 
bringing down the iV max required for convergence in larger A systems. At small iV max 
the size dependence on A is roughly linear. As iV max increases we get back into the 
exponential region. 

Again we must check on the scaling behavior of the SRG evolution with increasing 
basis size. The left panel of Fig. ID. 31 shows the evolution time in seconds as a function 
of the linear flow parameter, s, for several sample iV max values. All are linear in s, 
again indicating no stiffening of the problem. On the right panel we display the same 
quantities in a log-log plot for better visibility. The slope of a straight line in this type 
of plot indicates the exponent of a power-law dependence, here the slope is unity and 
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Figure D.3: SRG timings vs the parameter s show no stiffening for large s (small A). 
The log-log plot on the right gives a better view of the low iV max samples. 

the limits are set to make this apparent. The vertical shift indicates the difference in 
the slopes on the linear scale. 

D.3 Three Dimensional Jacobi NCSM 

In the left plot in Fig. ID. 41 we can see a similar scaling of basis size with iV max as 
in the one- dimensional version. This is yet another area where the one- dimensional 
model mirrors the behavior of realistic calculations and lends credence to its predictive 
power. The right panel of Fig. ID. 41 shows the SRG evolution is again robustly linear 
in s despite the different and complicated basis and embedding procedures of the 
realistic NCSM. This data had to be obtained from file write times so the curves are 
shifted vertically by an indeterminate amount of disk writing time. 

Practical coding considerations have so far prevented realistic calculations in the 
relative Jacobi NCSM for nuclei larger than A = 4. As covered in appendix [C] the 
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Figure D.4: The left plot shows scaling of NCSM basis size with iV max for the two- 
, three-, and four-body systems. The right plot shows evolution timings vs. s in 
realistic NCSM calculations for selected A = 3 partial waves. The legend indicates 
the dimension of the matrices being evolved. 



calculation of larger systems is accomplished by converting the Jacobi matrix ele- 
ments into the single-particle, or m-scheme, basis. This scaling is discussed in the 
next section with Fig. ID. 51 These problems have nothing to do with SRG implemen- 
tation, but they have an incidental impact on how quickly progress can be made. 
Implementing the SRG independently in the m-scheme basis is a different body of 
work and will be an important check though the return may not be high for evolving 
extremely large sparse matrices. The initial potential would have to be in a basis close 
to convergence; previous work using the SRG within m-scheme calculations [35| 152] 
did the evolutions in the momentum basis, which is effectively iV max = oo. 
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Figure D.5: Size of m-scheme basis as a function of iV max for several light nuclei. 
Note the logarithmic scale in basis size. 
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D.4 Lab Frame in One and Three Dimensions 

In the decoupling studies of chapter [3] we used evolved NN forces embedded in 
the single particle NCSM, referred to as m-scheme, to compute light nuclei. It is im- 
portant to note how much benefit the SRG provides to these particular calculations. 
Figure ID. 51 shows a plot on a log scale of the basis size verses N max for several lighter 
nuclei up to A = 7. The boxes indicate the extent of computational abilities for 
various systems. The box marked "VJowk Cluster - with Disk storage of Hamilto- 
nian" indicates the limits of our abilities with our local cluster of 32 processors with 
the limitation of writing the large files to disk, which is a considerable slowdown. 
Also the basis dimension is not quite as daunting as it seems since these matrices 
are fairly sparse; the number of actual nonzero elements is in the range of the basis 
dimension itself. Still this is a very large number of matrix elements to compute and 
store. Clearly we would like to have converged results for the N max values within the 
larger box marked "In-Core VJowk" so that we can perform the calculations without 
expensive supercomputer resources. As shown in section 13.51 the decoupling afforded 
by the SRG enables these calculations to achieve just such a convergence. 

Figure ID. 61 shows the scaling of basis sizes in the one-dimensional lab-frame basis. 
Note the dramatic increase in basis size over the one- dimensional Jacobi basis. This 
calculation has not been developed further due to its poor scaling properties as com- 
pared to the Jacobi version. However, it is an important check on the accuracy of 
calculations performed in the Jacobi basis, and further developments for this purpose 
will be straightforward as discussed in section IF. 41 

In one dimension, the lab-frame basis is not as efficient as the Jacobi basis in 
achieving convergence with basis size and therefore the Jacobi basis is preferred. 
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Figure D.6: A plot showing the scaling of One-D NCSM lab-frame basis or "m- 
scheme" size as a function of iV max for A = 2, 3, 4 bosons. Compare to Jacobi plots 
of Jacobi basis size in Fig. ID. 21 



However, in three-dimensions, angular momentum considerations change the situa- 
tion. In the Jacobi basis for a given A-body calculation, one must include multiple 
angular momentum channels, or partial waves, to obtain accurate results. The cod- 
ing involved in correctly accounting for all the necessary channels is apparently very 
involved due to the iterative procedure that builds on successive symmetrized A-body 
clusters. In contrast, the m-scheme strategy builds an A-body basis and symmetrizer 
from scratch and embeds the two- and three-body interactions directly into that space. 
Therefore the balance is tipped in favor of using large supercomputing to handle the 
large bases of the m-scheme, instead of continuing to develop code in the Jacobi basis. 
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APPENDIX E 



SPURIOUS STATES FROM G s = H ho 



In this appendix we will consider the spurious bound states which appear during 
the evolution with certain choices of the SRG generator, G s . In the one dimensional 
model, we explored various features of these spurious states qualitatively through 
observing the evolution of potential matrix elements in the ID basis. Watching the 
flow of individual elements can be very instructive when trying to understand the 
qualitative behavior of the SRG. Here we plot various potentials in the oscillator 
basis at N max = 12, which is convenient for looking directly at the potentials and 
seeing the interplay among matrix elements. 

Figure IE. II shows a film strip of evolving potentials in the three-body basis of the 
one-dimensional model. The top shows evolution using the choice G s = T re i and the 
bottom uses G s = H^ Q . The advantage of if ho over T rc \ with respect to the amount 
of diagonalization achieved in this basis is obvious as noted in chapter |H However, 
problems arise when we try to isolate the evolving two- and three-body forces. 

Figure IE. 21 shows the result of evolving the two-body potential using G s — Hho in 
that basis before embedding in the A = 3 space and computing the two-body-only 
contribution to the A = 3 bound state. The resulting evolution curve is analogous 
to that of Fig. 14.51 where we chose G s = T vc \ and the result was a nice smooth 
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Figure E.l: A series of plots of the oscillator basis Hamiltonian in a reduced N mS/X 
basis as it is evolved using G s = T rc i (top) and G s = H ho (bottom) to A = 7, 5, 4, 3, 2. 
The size of the basis shown is only iV max = 12 for visibility. Larger bases look 
qualitatively similar. 




Figure E.2: An example of the discontinuous curves of two-body-only calculations 
characteristic of the G s = H ho SRG. The two-body-only curves computed by evolving 
the A = 2 potential using G s = and embedding the result in the A = 3 basis 
to obtain the three-boson binding energy. The straight line labeled "2+3-body" is a 
check on the unitary evolution in the A = 3 space. 
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Table E.l: A sampling of evolving A = 3 spectra under the choice G s = if ho- The 
table uses values at N max = 20, which is large enough to display the effect, but small 
enough to display the entire spectrum. 



state 


A — 1UU 


A — O 


A — 4 


A — o 


A — Z 


1 


-2.547 


-2.772 


-2.920 


-10.156 


-7.597 


2 


0.049 


-0.295 


-0.456 


-10.090 


-7.497 


3 


1.713 


1.192 


0.370 


-3.377 


-5.129 


4 


3.103 


2.479 


0.760 


-2.994 


-4.199 


5 


5.709 


4.847 


2.126 


-0.895 


-3.649 


6 


6.615 


5.595 


2.701 


-0.571 


-3.544 


7 


7.525 


6.830 


3.771 


0.672 


-1.171 


8 


10.828 


9.221 


5.090 


1.113 


-0.878 


9 


12.276 


10.627 


7.192 


1.634 


0.114 


10 


13.401 


11.291 


7.670 


2.614 


0.190 



curve indicating the effect of induced three-body forces, in the lone bound state, as a 
function of A. Now, using G s = ifho, the curves are discontinuous 22 ! and the magnitude 
of induced three-body forces is much larger than expected. This plot turns out to 
be contaminated by the appearance of multiple three-body bound states due to the 
evolution of the two-body matrices with this G s . 

Table IE. 1 1 shows a sample of the first few states in the spectrum as the initial 
Hamiltonian is evolved with G s = H^. The evolution proceeds as expected up to a 
A value where new states begin to appear. As A is increased further all the bound 
states deepen dramatically, though the precise behavior is dependent on the choice 
of G s . 

To further explore the origins of these states, we looked at the NN-only A = 3 

matrices as they are evolved and embedded. Figure IE. 31 compares the choices of T re i 

22 Notc that the lines in Fig. IE. 21 are drawn to guide the eye from one point to the next. More 
sample A's may reveal this drop to be simply very steep. 
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Figure E.3: A series of plots of the A = 3 oscillator basis Hamiltonian in a reduced 
iV max basis as it is evolved using G s = if ho- 




Figure E.4: A series of plots of the A = 2 oscillator basis Hamiltonian in a reduced 
-^Vmax basis as it is evolved using G s = if ho- 



and ifho with film strips of these. Notice that while the T ie \ version converges on 
a shape similar to T rc i itself, the version moves toward the diagonal and then 
recedes from it, ending up at a shape similar to the T re i component but of opposite 
sign. This suggests that the fully unitary transformations are kept constant by means 
of large three-body forces being induced to compensate for some unknown effect in 
transformation of the two-body forces. 

In figure IE. 41 a film strip is shown of evolving A = 2, matrices corresponding 
to the A = 3 evolutions shown above. Notice how some matrix elements along the 
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diagonal mysteriously decrease so that they are not being driven to a monotonically 
increasing pattern as set by G s = H ho . As seen in chapter [31 the SRG tries to drive 
the matrix toward a shape dictated by the choice of G s . With the matrix elements 
out of order according to the SRG form, it appears to induce strength in off-diagonal 
matrix elements coupling these smaller states. This in turn demands large off-diagonal 
three-body forces to be induced to compensate for this reordering. The exact cause 
of the disrupted ordering of states is not well understood presently but is a subject 
of ongoing investigation. 

Moving up to the next many-body sector, we can look at the effect of the H^ Q 
choice in a four-body system. In Table IE.2I we show the evolving spectra for the 
two-body-only and two-plus-three calculations of the A = 4 system. The former is on 
the top block of the table and the latter is on the bottom. In the top we can see the 
effect is magnified in the larger system. In the bottom the onset of the spurious states 
appears to be earlier, which is consistent with the different shape of two-plus-three 
body calculations of the four-body system. That is, the induced four-body force may 
be compensating for additional spurious contributions occurring in the ^4 = 3 sector. 

In order to study the scale at which these spurious states appear, we also tried 
setting the potential to zero, leaving the initial Hamiltonian to be H s= q = T re i. We 
find that the spurious states appear suddenly at a specific value of the evolution 
parameter, A = \fhuj. This is expected since with V = 0, and most physical constants 
set to one, fko is the only scale left in the problem, having units of energy or [A] 2 . 
By increasing the sampling of A around the value, y/huj, we find that the onset of 
spurious states is indeed discontinuous, at least in the case of zero initial potential. 
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Table E.2: A sampling of evolving spectra for A = 4 under the choice G s = H^. 
Again the table uses values at iV max = 20. The upper block displays the two-body- 
only calculation while the lower block includes the two- and induced three-body forces. 



state 


A = 100 


A = 5 


A = 4 


A = 3.5 


A = 3 


A = 2 


1 


-4.609 


-5.496 


-9.697 


-42.757 


-10.156 


-33.059 


2 


-1.558 


-2.752 


-8.415 


-34.195 


-10.090 


-27.979 


3 


0.812 


-0.792 


-6.021 


-34.112 


-3.377 


-27.549 


4 


2.254 


0.446 


-3.215 


-22.265 


-2.994 


-27.152 


5 


2.733 


0.990 


-2.301 


-21.183 


-0.895 


-26.934 


6 


5.483 


2.999 


-1.974 


-20.878 


-0.571 


-25.851 


7 


6.742 


4.252 


-1.912 


-17.088 


0.672 


-25.692 


8 


7.620 


5.152 


-0.665 


-16.396 


1.113 


-23.746 


9 


8.077 


5.741 


0.147 


-15.059 


1.634 


-22.920 


10 


9.372 


6.465 


0.951 


-14.269 


2.614 


-21.483 


1 


-4.609 


-4.674 


-29.265 


-11.822 


-6.029 


-4.039 


2 


-1.558 


-2.969 


-11.135 


-5.366 


-4.016 


0.745 


3 


0.812 


-0.235 


-8.177 


-5.155 


-3.004 


1.856 


4 


2.254 


0.735 


-7.319 


-4.190 


-1.838 


1.983 


5 


2.733 


2.475 


-4.210 


0.154 


0.398 


5.157 


6 


5.483 


5.361 


0.083 


2.472 


2.378 


8.707 


7 


6.742 


5.529 


1.500 


3.532 


4.328 


11.501 


8 


7.621 


6.095 


2.001 


4.142 


6.327 


11.619 


9 


8.077 


7.414 


6.263 


5.964 


8.342 


13.510 


10 


9.372 


7.473 


6.382 


7.188 


9.963 


16.678 




Figure E.5: A series of plots of the A = 2 oscillator basis Hamiltonian in a reduced 
N max basis as it is evolved using G s = H d osc . 



200 



Table E.3: A sampling of evolving 3 H spectra under the choice G s = in the 
realistic three-dimensional NCSM. The table uses values at iV max = 20 which is large 
enough to display the effect, but small enough to display the entire spectrum. 



state 


A = 100 


A = 5 


A = 2 


A = 1.5 


A = 1.2 


A = 1 


1 


-7.4758 


-7.5540 


-8.7177 


-10.0938 


-145.6345 


-118.5970 


2 


2.1186 


2.0412 


0.9895 


-0.8918 


-130.2285 


-116.3076 


3 


6.5599 


6.4646 


5.1126 


3.1332 


-123.0114 


-109.9408 


4 


6.8592 


6.8040 


5.8956 


4.5028 


-106.4009 


-91.1749 


5 


12.7642 


12.7365 


11.0206 


6.5762 


-105.4125 


-88.3281 


6 


13.5303 


13.3832 


11.7719 


10.0159 


-97.4654 


-87.8343 


7 


15.7025 


15.6187 


13.9796 


10.6716 


-94.8004 


-86.9702 


8 


16.4911 


16.4027 


14.7041 


11.5136 


-88.8669 


-84.5732 



We also tried other choices of G s that are diagonal in the oscillator basis. In 
Fig. IE.5l we show the evolution the A = 2 Hamiltonian using G s = Hd tOSC , where i?d,osc 
is the diagonal part of the evolving Hamiltonian in the oscillator basis. Here we can 
still see the matrix elements evolving out of order and subsequent development of 
large off-diagonal elements in the evolution to push strength to low momentum. Here 
the effect is more dramatic than in the ifho case, since the moving strength is jumping 
several states. This again forces large induced three-body forces to compensate for 
these evolved short range interactions. 

We repeated this experiment in the three-dimensional NCSM code for calculations 
of the triton. The scale of the spurious states, relative to the initial bound states in 
that calculation, was similarly large. The spurious states appeared in the NN-only 
calculations in qualitatively similar ways. Table lR3l recreates the results of table lR2l 
using the three-dimensional NCSM. The first few states in the spectrum are shown 
at selected A's. 
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As can be seen in Fig. IE.5l the SRG evolution on the two-body Hamiltonian seems 
to be driving information from high to low iV max in an off-diagonal way. This pattern 
is indicative of spurious states being developed as the SRG builds large strength 
matrix elements at high-momenta and compensates for that by inducing large three- 
body forces. 

So, while the lone constraint on r] s of anti-hermiticity grants a freedom of choice in 
choosing the form, G s , of evolution, we must be careful of the effects that choice may 
have on the physics during renormalization of different parts of the interaction. Work 
is currently underway to investigate various choices of G s which achieve convergence 
not only for binding energies but for other long-range observables such as the RMS 
radius. [12] 
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APPENDIX F 



SINGLE PARTICLE COORDINATE OSCILLATOR BASIS 



An important check on the Jacobi oscillator basis calculations done in this thesis 
is to repeat them in a different basis. An oscillator basis in single particle, or lab- 
frame coordinates has already been developed in three-dimensions and is used to 
compute nuclei above the practical range of the Jacobi basis. Here we develop a one- 
dimensional analog to the single-particle basis in the spirit of Ref. [95] and compare 
the results with our Jacobi basis calculations. 

The general strategy used here for the lab-frame oscillator basis approach to the 
many-body problem is very similar to the Jacobi basis version. The allowed basis 
states are explicitly constructed and the states' quantum numbers listed. The nuclear 
interaction is embedded in the constructed basis. A symmetrization operator is built 
to isolate the physical states of the system. The symmetric A-body Hamiltonian is 
then diagonalized to obtain the energy spectrum. However there are several important 
differences in each of these steps; we will treat them in turn. 

F.l Embedding 

Here we are working exclusively in single particle coordinates, k\ and k 2 , but we 
have the potential in terms of Jacobi momenta because we already had a function 
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for it from the previous one-dimensional work. We must first convert our interaction 
from the Jacobi basis to the single particle basis. We will use the Talmi-Moshinsky 
transformation brackets that we are already familiar with from the Jacobi basis in 
appendix O Formally, we can express this transformation of the potential as: 

(nin 2 |F|nX) = {n x n 2 \Nn)(Nn\V\ N'n'^N'n^n'^ 

= (nin2\Nn)SN,N'{n\p)V(p,p')(p\n)(N'n'\n' 1 n 2 ) (F.l) 

In our MATLAB implementation, Eq. (IF. If) is accomplished by matrix multiplica- 
tions. The matrix representing the transformation from momenta to oscillators, 
(n\p), is the same set of single-coordinate oscillator wavefunctions as used previ- 
ously. The matrix representation of V(p, p') is read in directly from the expression of 
Eq. (14.31) as before. The delta function, 5n,n>, is embedded so that we are working in 
a full two particle basis organized by the total oscillator number, N to t = N + n (i.e., 
\Nn) = |00), |01), |10), |02), |11), |20), |03>, ...)■ The transformation brackets used here 
can be obtained from Eq. (IC.9I) by using d = 1 instead of d = 3 used there. This 
gives the rotation from the Jacobi center of mass, and first relative coordinate to the 
two single-particle coordinates. 

As in the Jacobi case, the relative kinetic energy can be built from creation and 
annihilation operators: 

/ \rp I / / \ / I (^1 ~~ ^2) 2 \ I l\ 

(nin 2 T re i K n 2/ = ™i n 2 o \ n i n 2> 

2m 

= ———{n 1 n2\{a[-a 1 -a! 2 + a2)\n 1 n 2 ). (F.2) 

Or, we can take the kinetic energy as built in the Jacobi basis and convert using the 
transformation brackets as for the potential above. These two approaches are identical 
and important checks that our construction is consistent. However, for A > 2 the 
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only practical option is the first option if we do not want to build transformation 
brackets for each A-body space conversion. 

To calculate properties of an A-body system, we must be able to embed the two-, 
three-, and higher-body forces into the A-body space. First we build a list of the 
possible states in that space, \nin2ns...nA), and organize them into groups of total 
N = Ylt=i n i U P to the maximum, iV max . For example, for A = 3 the list starts out 
as: 0,0,0; 0,0,1; 0,1,0; 1,0,0; 0,0,2; 0,1,1; 0,2,0; 1,0,1; 1,1,0; 2,0,0; etc up to N = iV max • 

It is a relatively simple matter to look through the lab-frame states and place 
matrix elements of the two-body lab-frame interaction in the appropriate matrix 
elements of the A-body lab-frame basis. For a two-body interaction, we look through 
every pair of incoming and outgoing particles in every matrix element of the A-body 
basis and require a delta function in the initial and final energy of all other particles in 
the system. If these conditions are met, the appropriate two-body interaction matrix 
element is added to the A-body representation of V^ 2 \ A similar routine must be 
coded separately for embedding the three-body forces. This routine would likewise 
consider all triplets of particles and require a delta function in the others. This 
method of embedding is not the most efficient especially with regard to obtaining 
matrix elements between symmetric states as discussed below. 

F.2 Symmetrization 

In lab-frame coordinates, to get a symmetric state of bosons (or antisymmetric 
for fermions), we must simply take the permanent (determinant) of matrices in the 
full oscillator representation. For example, the symmetric state with one particle in 
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an oscillator state rij = 2 and the other with oscillator would be: 

(<02| + (20|)1/ (2) (|02) + |20))/2 = (1/(02; 02) + 1/(02; 20) + 1/(20; 02) + 1/(20; 20))/2; , 

(F.3) 

where the 2 is the product of both normalization factors, 1 / a/2. 

Here we have used existing machinery to accomplish the permanent by building 
a symmetrizer matrix, S, and requiring that it be a projector, S 2 = S. Instead of 
the Jacobi version, where we needed some complicated combination of transformation 
brackets to build S, simple factors of one over the number of permutations of a given 
state ensure that S is a projector for the single particle basis. By keeping the eigen- 
states of this matrix with eigenvalue one (and discarding those with eigenvalue zero), 
we obtain a non-square matrix analogous to the "coefficients of fractional parentage" 
used in the Jacobi case. The Hamiltonian and other operators expressed in the full 
single particle basis can be symmetrized using the physical eigenstate matrix just as 
done in the Jacobi case to obtain an operator in that basis. 

For A = 2, both of these matrices can be seen in Fig. IF. 11 The symmetrizer shown 
on top has a simple structure with factors of 1/2 for those states with 2 permutations 
(i.e.: 1 02) and 1 20) ) and 1 for those which are already symmetric like (i.e.: |11)). Note 
that the eigenvectors of displayed in the bottom panel, have the correct normalization 
factors, 1 / y/2 for a symmetrized two-particle state. 

For A > 2, we can build a symmetrizer operator from the list of total possible 
states discussed above. Now we have to be more careful of the permutations of 
states involved in symmetrizing. The number of permutations of a given state is 
computed and recorded as the list is built. By counting the number of particles 
with a certain oscillator number, the number of permutations of the state is given 
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100 ' ' ' ' 1 1 1 1 1 ' ^"-1 

10 20 30 40 50 60 70 80 90 100 




5 10 15 20 25 30 35 40 45 50 



Figure F.l: A sample of the symmetrizer for the 2-body lab- frame problem on the 
top, and the corresponding non-zero eigenstates of the symmetrizer on the bottom. 
Here N max = 12 as a more realistic iV max would not be visible on the page. 
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by (^4!/ EL max ki) where ki is the number of particles in the oscillator state i. To 
build the symmetrizer, we simply loop through the matrix elements of the full single- 
particle basis and compare the initial and final states. If one state is a permutation 
of the other and their number of permutations are recorded as the same, then one 
over that number is recorded for that matrix element of the symmetrizer. 

Figure IF.2I shows an example of the symmetrizer and the corresponding symmetric 
eigenstates for the A = 3 system up to iV max = 6. Note the smaller colors compared 
to the A = 2 case because with more particles, there are more permutations of a 
given energy content. For instance, the faintest yellow corresponds to a value of 1/6 
for states with 6 permutations (i.e.: 1 021) ) . The darker yellow is 1/3 for states like 
|112). 

The matrices displayed in Figs. lF.ll and lF.2l for the physical symmetric eigenstates 
are the lab-frame analog to the coefficients of fractional parentage used in the Jacobi 
basis. Now, however, the construction of an A-body basis does not rely on an it- 
erative procedure on the A — 1 symmetrized basis. We can build the A-body basis 
from scratch, embed the two- and three- (and higher) body forces in that basis, and 
symmetrize them with the A-body symmetrizer built as just described. 

This approach to building a symmetrized basis for lab-frame momenta is not the 
same as employed in three-dimensions. The overhead of building a large set of states 
and distilling out the symmetric states is very large in three-dimensions. Convergence 
is only feasible in one-dimensional calculations because of an overall reduction in basis 
size do to the omission of angular momentum quantum numbers. However, this is 
a useful check on the general procedure used in the Jacobi case, and can still be a 
useful tool for scaling up in A in one dimension. The method used here of building 
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a relatively simple symmetrizer operator can also be viewed as a first step to check 
future algorithms which may be less straightforward to code. 

F.3 Center of Mass Separation 

A major difference from the Jacobi basis approach is the presence of the center 
of mass energy in the single-particle basis. For each value of the total center of mass 
energy, the A-body system has the same intrinsic spectrum. However, in our basis 
these copies of the spectrum are mixed up together, so we must separate them to 
obtain meaningful results. We can add some quantity which will boost the unwanted 
copies up in energy and isolate them. The harmonic oscillator Hamiltonian, H cm , has 
the straightforward behavior in the oscillator basis that it doesn't mix up the intrinsic 
physics. Therefore, we can subtract off its effects to regain the desired spectrum. 

So, in order to isolate the different center of mass states we will have to add 

P(H cm -(n + l/2)frw) (FA) 

to the total interaction. The constant f3 is a large positive number which will provide 
large separation between the spectra at each of the center-of-mass eigenstates denoted 
by n (n = is the center-of-mass ground state). The value of n determines which 
center-of-mass solution we will isolate exactly, free of factors of /3huj. Spectra at 
center of masses less than n will be kicked downward and those greater than n will 
be kicked upward. 
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Writing out the center-of-mass Hamiltonian using creation and annihilation oper- 
ators is a simple matter in the general A-body space: 

(n 1 n 2 ...n A \H cm \n 1 n 2 ...n A ) 

hw A A 

(n 1 n 2 ...n A \(^al + a*) 2 - (y]a\ - a i ) 2 )\n! 1 n' 2 ...n' A ) 



8mA 

i=l i=l 

(n 1 n 2 ...n A \2 ^(aja^ + a i a ] j )\n' l n 2 ...n' A ) 

A A A 



8mA 

hw 
8mA 

TlhJ 

-{n 1 n 2 ...n A \'2 

i=l 



8mA 

A A 



{nin 2 ...n A \2 s ^(a\a i + did]) ^ (a\a j )\n' 1 n' 2 ...n' A ) 

i=l i=l j=i+l 
A 

(n 1 n 2 ...n A \2 ^(2n- + l)5 nijn >. 



+ 4 E E v^v^+i)^ ' i+ i,n l fin' 3 -i,n J \n 1 n 2 ...n A ) . (F.5) 

i=l j=i+l 

This is written in the same basis as the potential and kinetic energy and can be 
symmetrized by the same process as described in section IF.21 

Table IF. 11 shows an example of the separation of different center of mass spectra 
due to H cm . Columns showing the labels of the center-of-mass and intrinsic states 
are meant to aid the reader's eye when viewing these spectra. The left shows the 
spectra from different center-of-mass energies entangled. For instance, all the ground 
states are collected together as the lowest energies and next the first excited states. 
Separating these by hand is not a trivial matter, especially in larger systems. In 
contrast, the right side shows the same spectrum but now separated by a strong j3 
term. The spectrum from the lowest center of mass appears at the head of the list. 
The next spectrum is multiplied by /3ftw (here (3 = 1000 and hio = 5), the third by 
2j3friu, etc. 
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Table F.l: A sample of the spectrum resulting from the lab-frame basis at a small 
iV max before and after separation by a H cm term. The columns right of each list of 
energies show the center of mass and intrinsic energy level numbers, n and /. The 
last column shows the effective N max for the right-most spectra. The value of (3hjj 
here is 5000. Note the "separated" values should be scaled by 10 4 . 



mixed 


{n 


J) 


separated (xlO 4 ) 


(n 


,0 


eff. N max 


-0.812294429336271 


(0, 


0) 


-0.000081229442933 


(0, 


0) 


6 


-0.487362126673450 


(1, 


0) 


0.000409822038397 


(0. 


1) 


6 


-0.487362126673446 


(2, 


0) 


0.001080597080727 


(0. 


2) 


6 


-U.4ozi f yybbz 1 4ioo 


[6, 


U ) 


U.UUzbo iuZlDoo IDl 


{(), 


6 ) 


/:' 



-0.482179966274182 


(4, 


o) 


0.499951263787332 


(1. 


0) 


5 


4.098220383969666 


(o, 


1) 


0.500489798221359 


(1. 


1) 


5 


4.817772759274547 


(5, 


o) 


0.501990571277867 


(1. 


2) 


5 


4.817772759274551 


(6, 


o) 


0.999951263787332 


(2. 


0) 


4 


4.897982213589284 


(1, 


1) 


1.000489798221359 


(2. 


1) 


4 


4.897982213589287 


(2, 


1) 


1.001990571277867 


(2. 


2) 


4 


10.805970807270517 


(o, 


2) 


1.499951782003372 


(3. 


0) 


3 


12.846759311682062 


(4, 


1) 


1.501284675931167 


(3. 


1) 


3 


12.846759311682076 


(3, 


1) 


1.999951782003373 


(4. 


0) 


2 


19.905712778674406 


(1, 


2) 


2.001284675931168 


(4. 


1) 


2 


19.905712778674410 


(2, 


2) 


2.500481777275925 


(5. 


0) 


1 


26.870215837516941 


(o, 


3) 


3.000481777275927 


(6. 


0) 






Each of the spectra obtained in this way are identical in a complete Nfko space, 
but we must work in a truncated space. The most accurate spectrum is that of n = 
because the other spectra begin at larger and larger center of mass energies and have 
fewer basis states available to contain their excited states. In this way we can see the 
effect of truncation in A^ max as we look up in center-of-mass energy. There are fewer 
and fewer states in the spectrum until the center-of-mass energy is A^ max and there is 
no more "space" for excitations in the basis. Note that this truncation is identical to 
that done explicitly in the Jacobi calculations. On the right side of table IFTTI we show 
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the effective N max for each center-of-mass spectrum. The ground states for each of 
these spectra are identical to the ground states from the Jacobi calculations truncated 
to the iV max shown in the last column. 

F.4 Evolution Results 

Just as in the Jacobi basis, implementation of the SRG's flow equations is a 
straightforward matter in the lab-frame coordinates. The Hamiltonian and generator 
are given to a differential equation solver and the resulting evolved Hamiltonian is 
diagonalized to produce a spectrum of states for that A-bodj system. Also as in 
the Jacobi basis, evolved n-body forces must be isolated and embedded with the 
appropriate symmetry factor into the A-body space. At present, existing algorithms 
are only sufficient for embedding two-body operators into an A-body space. This 
limits us to NN-only and fully unitary calculations of A-body spaces. Additional 
codes are needed before the m-scheme version can be used to check the hierarchy 
of induced many-body forces in systems with A > 3. However, for A = 3 we can 
perform the calculation analogous to that of Fig. 14.51 and find identical results to 
within expected numerical errors. Further coding to check the hierarchy for larger A 
is a straightforward next step for this tool. 



213 



